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Table 12 on page 66 should read 


Mode 

Buckling load parameter 

STARS solution 

Exact 

solution 

1 
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Line 6.1.2 on page 80 should read 


6.1.2 ISTMN, NLVN, GR 
Format (215, E10.4) 
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The first seven lines of ASE input data on page 136 should read 


AERO TEST MDEL - ANTT SYMMETRIC CASE - NCNTRL = 5 
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REVISED RIGID BCDY Y TRANSLATIONS, ROLL, YAW, PLUS AILERON AND RUDDER MODES. 
NACH NO.- 0.90 ALTITUDE: SEA LEVEL 

1 17 3 10 100000 


The k-type flutter input data on page 140-141 should read 
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Table 17 in chapter 7 on page 141 should read 


Table 17. AERO test model: An aeroelastic antisymmetric analysis using a direct interpolation 

for AERO paneling. 


Mode 

Instability 

number 

k-SOLN 

p-k 

ASE 

Velocity, 

keas 

Frequency, 

rad/sec 

Velocity, 

keas 

Frequency, 

rad/sec 

Velocity, 

keas 

Frequency, 

rad/sec 

Fuselage first bending 

FI 

506.15 

77.49 

533.53 

77.54 

591.16 

76.03 

Wing first bending 

F2 

615.75 

7226 

620.85 

73.02 

579.69 

80.09 

Fin first bending 

F3 

— 

— 

856.64 

138.92 

742.81 

128.87 


Analysis notes: 

1 ) F - Flutter point 

2) D - Divergence point 

3) Mach = 0.90 

4) Altitude = Sea level 
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Velocity, V, keas 





Velocity, V, keas 


Figure 31. STARS ATM-pk flutter analysis — damping (g 1 ), frequency (B), velocity (v) plot, 
antisymmetric case, using direct interpolation where g’ = g x 200. 
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SUMMARY 


A unified and highly graphics-oriented analysis capability that includes structures, aerodynamics, 
and control disciplines has been achieved by integrating two new modules in the original STARS 
(S7ructural Analysis Routine^) finite element computer program system (ref. 1). Thus, in addition to 
performing static, buckling, vibration, and dynamic response analyses of undamped and damped 
structures, including those having rotating and prestressed components, the program is also capable 
of performing aeroelastic (flutter and divergence) and aeroservoelastic stability analyses that may 
include digital or analog controller elements. 

The element library in the STARS-SOLIDS module consists of one-dimensional (1-D) line ele- 
ments; two-dimensional (2-D) triangular and quadrilateral shell elements; and three-dimensional 
(3-D) tetrahedral and hexahedral solid elements. Improved composite as well as sandwich shell 
elements are recent additions to this element library. These elements enable the solution of 
structural problems that include truss, beam, space frame, plane, plate, shell, and solid structures, or 
any combination thereof. Associated algebraic equations are solved by exploiting inherent matrix 
sparsity. Zero, finite, and interdependent deflection boundary conditions can be implemented by the 
program. The associated dynamic response analysis capability provides for initial deformation and 
velocity inputs, whereas the transient excitation may be either forces or accelerations. An effective 
in-core or out-of-core solution strategy is automatically employed by the program, depending on the 
size of the problem. Data input may be at random within a data set, and the program offers certain 
automatic data generation features. Input data are formatted as an optimal combination of free and 
fixed formats. 

The STARS-AERO module enables computation of unsteady aerodynamic forces pertaining to 
subsonic and supersonic flow employing the doublet lattice and the constant pressure procedures, 
respectively. Subsequent flutter and divergence analyses are performed by utilizing the standard k 
and p-k methods. 

Associated aeroservoelastic (ASE) analysis is performed by the STARS -ASE module that 
employs a Pade and least squares formulation for curve fitting the unsteady aerodynamic forces, 
followed by a state-space formulation that incorporates the elastic structural and aerodynamic 
effects. Such state-space matrices are further augmented by analog and/or digital controller 
elements. A hybrid formulation is adopted in the case of a digital controller. Frequency response 
characteristics may be calculated to yield frequency, damping, gain, and phase results for both open- 
and closed-loop systems. 

The program, developed in modular form for easy modification, is written in standard FORTRAN 
to run on a variety of computers such as IBM 6000, DEC VAX 3600, Cray Y-MP, and Cray 2, among 
others. Extensive graphics capability exists for convenient model development as well as for 
postprocessing of analysis results, employing any PLOT- 10 compatible terminal such as Tektronix, 
CIT, DEC-VT, and also the E/S PS 390 machines. 
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1. INTRODUCTION 


The general-purpose digital computer program, STARS (STructural Analysis Routine^), has been 
designed as an efficient tool for analyzing practical structures, as well as for supporting relevant 
research and development activities; it has also proved to be an effective teaching aid. All such 
activities are mutually enhancing and interrelated (fig. 1). The current version of the program, 
capable of solving linear elastic aero-structural-control problems (fig. 2), will be continuously 
updated to include other forms of analysis. 

In an effort to optimize the program layout, the subroutines in the various modules have been 
grouped into several links. Interaction between the user and the program is effected through a 
display terminal with or without graphics capabilities; however, a graphics terminal is useful in the 
accurate preparation of data input and in visualizing structural geometry and analysis results. Thus 
with reference to figure 3, which pertains to the typical STARS-SOLIDS module, link 1 relates to the 
input phase of the program. The user may input the data directly, utilizing some limited, automated, 
data-generation capabilities inherent in the program; alternatively, the STARS preprocessor 
program, PREPROC, may be utilized effectively for complete data generation of a structure. Once 
the data have been entered into the system, the user may create an image of the model on the 
terminal display screen. 

Subsequent correction or modification of the model may be easily implemented on an interactive 
basis. Once the model format is verified, the user may proceed to run link 2 of the program, which 
involves major numerical manipulation of input data relative to static, stability, and free vibration 
analysis of the structural model. Nodal displacements caused by static loading and the structural 
mode shapes pertaining to the stability and free vibration problems may then be displayed using the 
graphics terminal. Link 3 of the program, the response link, enables computation of structural 
displacements caused by dynamic loading, as well as element stresses resulting from static and 
dynamic loads input. The STARS postprocessor program, POSTPLOT, is used for extensive plotting 
of results generated by the two links. 

The program can solve static and dynamic problems of nonrotating and rotating structures of 
general configurations with arbitrary displacement boundary conditions. For static problems, a 
multiple set of input data is permissible; for dynamic response problems, a single set of force or 
acceleration data is the usual input. The structural material may be isotropic, orthotropic, or 
anisotropic. Both viscous and structural damping occurring in practice may be included in the 
dynamic analysis. A bandwidth minimization option is available, its utilization being highly desirable 
to ensure economical solution of associated problems. 

The free vibration and dynamic response analyses of structural systems rotating along an 
arbitrary axis are useful features of the STARS program. Such a structure may have a combination of 
nonrotating and rotating parts, and each part may have a different spin rate. Both rigid body and 
elastic modes may be computed by the program, and the dynamic response analysis is formulated 
accordingly. Relevant algebraic computations are performed in single or double precision, using 
either real or complex arithmetic operations. 
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Figure 2. Major modules of STARS, 
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Figure 3. STARS-SOLIDS overview. 


A schematic of the associated aeroelastic and aeroservoelastic (ASH) analyses is depicted in 
figure 4. Thus, once the frequencies and mode shapes of the structure are derived from finite element 
analysis employing the STARS-SOLIDS module, the STARS-AERO module is next utilized to 
compute the unsteady aerodynamic forces on the structure. An alternative option enables input of 
measured modal data in lieu of calculated data. A flutter solution is then achieved using the k and 
p-k methods. The user has to input details of the aerodynamic paneling to achieve the aeroelastic 
analysis. 

Subsequent ASE analysis of the structure may be achieved by first employing the STARS-ASE- 
PADE submodule. The user provides essential data to perform a polynomial curve fitting of 
unsteady aerodynamic forces resulting in the state-space matrices. For an alternative open- loop 
flutter analysis, such data consist of information on polynomial tension coefficients, previously 
calculated generalized masses, and damping and modal characteristics as well as a set of velocity 
values. Additional data, in lieu of velocity values, relating to coordinate transformations from earth - 
to body-centered coordinate systems and also sensor locations, are needed for the subsequent ASE 
analysis for frequency response calculations and also for determination of damping and frequency 
values. This is achieved by the STARS -ASE-FRESP submodule when the primary data input 
relates to analog and/or digital controller blocks connectivity, associated transfer function polynomial 
descriptions, as well as gain inputs, specifications for system outputs and inputs, and also 
connection details between the plant and the blocks. 

Section 2 provides a concise description of the STARS-SOLIDS module of the program, as well 
as highlights of some of its important features, and section 3 depicts the data input procedure. 
Section 4 provides summaries of input data and analysis results for a number of sample test cases 
relevant to the module. Chapter 5 describes the various features of the newly created aeroelastic 
(AERO) and ASE analyses capabilities, whereas chapter 6 provides data input details of various 
related submodules. Finally, a representative, integrated aero-structural-control sample problem is 
worked out in detail in chapter 7. 
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Figure 4. STARS -ASE flowchart. 
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2. STARS-SOLIDS PROGRAM DESCRIPTION 


The structure to be analyzed by STARS may be composed of any suitable combination of one-, 
two-, and three-dimensional (1-, 2-, and 3-D) elements. The general features of the STARS- 
SOLIDS module include the following: 

1. A general-purpose, compact, finite element program. 

2. Elements: bars, rods, beams, 3-D line elements, triangular and quadrilateral plane, plate, 

shells, as well as sandwich panels and composite elements, tetrahedral and hexahedral solids. 

3. Geometry: any relevant structure formed by a suitable combination of the elements in (2). 

4. Material: general, isotropic, orthotropic, and anisotropic material. 

5. Analysis: natural frequencies and mode shapes of usual and rotating structures with or 

without structural damping, viscous damping, or both, including initial load (prestress) effect; 
stability (buckling) analysis; dynamic-response analysis of usual and rotating structures; and static 
analysis for multiple sets of mechanical and thermal loading. 

Special features of the STARS program include the following: 

1. Random data input within a subset. 

2. Matrix bandwidth minimization. 

3. Automatic node and element generation. 

4. General nodal deflection boundary conditions. 

5. Multiple sets of static load input. 

6. Preprocessor and postprocessor. 

7. Plot of initial geometry. 

8. Plots of mode shapes, nodal deformations, and element stresses as a function of time, 
as required. 

Structural geometry is described in terms of the global and/or the local-global coordinate system 
(GCS/LGCS) having a right-handed Cartesian set of X-, Y-, and Z-coordinate axes. Each structural 
node is assumed to have six degrees of freedom (DOF) consisting of three translations, UX, UY, 
UZ, and three rotations, UXR, UYR, UZR, which are the undetermined quantities in the associated 
solution process. Details of some important features of the program are summarized below. 
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2.1 Nodal and Element Data Generation 

The STARS program provides simple linear interpolation schemes that enable automatic genera- 
tion of nodal and element data. Generation of nodal data is dependent on the occurrence of such 
features as nodes lying on straight lines and common nodal displacement boundary conditions, 
whereas generation of element data is possible if the finite element mesh is repetitive in nature with 
elements possessing common basic elemental properties. The program enables input of data 
employing a number of rectangular local-global coordinate systems (LGCS) relevant to various 
substructures. 

A separate preprocessor, PREPROC, has been developed for automated generation of nodal, ele- 
ment, and other associated input data for any continuum. The preprocessor is an interactive graphics 
structures modeling program. It is capable of generating complex structures through duplication, 
mirror imaging, and cross-sectioning of modular representative structures. 


2.2 Matrix Bandwidth Minimization 

This feature enables effective bandwidth minimization of the stiffness, inertia, and all other 
relevant system matrices by reordering input nodal numbers, taking into consideration first-order, as 
well as second-order, nodal connectivity conditions. With reference to figure 5, the existing nodal 
numbering may be modified (ref. 2) to minimize bandwidth of associated matrices. Therefore, any 
node with minimum first-order connectivity may be chosen as the starting node. Accordingly, any 
one of nodes 1, 4, 7, 10, 13, and 16, all of which have a minimum first-order nodal connectivity of two, 
may be selected as the first node to start the nodal numbering scheme. However, nodes 1,4, 10, and 
13 possess a higher second-order connectivity condition than do nodes 7 and 16. For example, 
nodes connected to node 1 (namely nodes 2 and 18) are, in turn, connected to a total of seven nodes, 
whereas such a connectivity number for either node 7 or 16 happens to be only six. As such, either 
node 7 or 16 may be chosen as the starting node for the renumbering scheme. A revised nodal 
numbering that minimizes matrix bandwidth is shown in parentheses in figure 5. The present 
minimization scheme also takes into consideration the presence of nodal interdependent 
displacement boundary conditions. 
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Figure 5. Bandwidth minimization scheme. 


2.3 Deflection Boundary Conditions 

The nodal displacement relationships may be classified as zero, finite, and interdependent 
deflection boundary conditions (ZDBC, FDBC, and IDBC). Details of such a formulation are 
provided in section 3.4. Thus, in addition to prescribed zero and finite displacements, the motion of 
any node in a particular degree of freedom can be related in any desired manner to the motion of the 
same or any other combination of nodes in any set of specified directions. 


2.4 Prescribed Loads 

A structure may be subjected to any suitable combination of mechanical and thermal loadings. 
The loads in the mechanical category may be either concentrated at nodes or distributed. Thus, 
uniform pressure may be applied along the length of line elements acting in the direction of the local 
y- and z-axes. Such uniform surface loads are assumed to act in the direction of the local z-axis of 
the shell and solid elements, acting respectively on the shell and solid base surfaces. 
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The effect of thermal loading can be incorporated by the appropriate input of data pertaining to 
uniform element temperature increases, as well as thermal gradients. 

2.5 Static Analysis 

Static analysis, performed by setting the parameter IPROB = 8 in the input data, is effected by 
solving the set of linear simultaneous equations 

KU = P (1) 

where 

K = system elastic stiffness matrix 

U = nodal displacement vector 

P = external nodal load vector 

IPROB = integer designating problem type (defined in section 3.1) 

A multiple set of load vectors is represented by the matrix P incorporating effects of both mechanical 
and thermal loading. The equations are solved once, initially by Gaussian elimination, and solutions 
pertaining to multiple nodal load cases are obtained by simple back substitution. 

2.6 Elastic Buckling Analysis 

A buckling analysis is performed by solving the eigenvalue problem 

(K E +yK G )U = 0 (2) 

in which K E and K G are elastic stiffness and geometric stiffness matrices, respectively; U represents 
the buckled mode shapes and y is the buckling load. This is achieved by setting IPROB = 9. 

2.7 Free Vibration Analysis 

The matrix equation of free vibration for the general case of a spinning structure with viscous and 
structural damping is expressed (ref. 3) as 

[K e (1 + i*g) + K 0 + K']U + (C c + C D )U + MU = 0 (3) 

in which a dot indicates differentiation with respect to time; the previously undefined terms are de- 
scribed as follows: 

K' = centrifugal force matrix 

C c = Coriolis matrix 

C D = viscous damping matrix 
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M = inertia matrix 

g = structural damping parameter 

i* = imaginary number, V— 1 

Such a structure may have individual nonrotating and also rotating components spinning with 
different spin rates along arbitrary axes. 

Various reduced sets of equations pertaining to specific cases of free vibration are given as 
follows: 

1. Free, undamped vibration of nonrotating structures (IPROB = 1): 

K e U + MU = 0 (4) 

2. Free, undamped vibration of spinning structures (IPROB = 2): 

KU + C c xl + MU = 0 (5) 


with K = K E + K G + K . 

3. Free, damped vibration of spinning structures (IPROB = 4, 5), defined by equation (3). 

4. Free, damped vibration of nonspinning structures (IPROB = 6, 7): 

K e (1 + i*g) U + C D U + MU = 0 (6) 

The eigenvalue problem pertaining to the IPROB = 1 and 9 cases is real in nature, but the rest of the 
above problems involve complex-conjugate roots and vectors. In the special case of a prestressed 
structure, the matrix K G is automatically included in equation (6). 

In addition, STARS solves the quadratic matrix eigenvalue problem (IPROB = 3) associated with 
a dynamic element formulation (ref. 4), 

[k e - xm - x\m 2 - K^iu = o (7) 

which is quadratic in terms of the eigenvalues % = X 2 and where both M 2 and K 4 are the higher 
order dynamic correction matrices, X being the natural frequencies. This option is currently being 
updated to include a number of elements. 

Prestressed structures caused by initial loads may also be analyzed, in which case the relevant 
eigenvalue problem for undamped structures has the form 

(K e + K g - XM)V = 0 (8) 
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in which the geometrical stiffness matrix K G is a function of initial stresses; similar formulations are 
obtained for structures with various forms of damping. 


2.8 Dynamic Response Analysis 

The modal superposition method is employed for the dynamic response analysis following the 
computation of structural frequencies and modes. As an example, for a nonrotating, undamped 
structure, the associated eigenvalue problem of equation (4) is first solved to obtain the first few 
eigenvectors O and also the eigenvalues. The vectors may consist of a set of rigid body modes O r 
and a number of elastic modes d> e which are next mass-orthonormalized so that the matrix product 

O t MO = [I] (9) 

is a unit matrix. A transformation relationship 

U = <Dri (10) 


is substituted in the dynamic equation 


MU + KU = P(t) (11) 

and when premultiplied by yields a set of uncoupled equations 

q r =O r T P(t) (12) 

and 

+ (13) 

incorporating rigid body and elastic mode effects, respectively; P(t) is the externally applied, time- 
dependent forcing function, and Q 2 is a diagonal matrix, with terms 0)i, co* being the natural frequen- 
cies. Solutions of equations (12) and (13) can be expressed in terms of Duhamel's integrals, which, 
in turn, may be evaluated by standard procedures (ref. 5). In the present analysis, the externally 
applied, time-dependent forcing function must be applied to the structure in appropriate small, incre- 
mental steps of rectangular pulses. The forcing function may be either load or acceleration vectors; 
the program also allows application of initial displacement and velocity vectors to the structure. For 
spinning, as well as for damped, structures identified as IPROB = 2, 4, 5, 6, and 7, <D T are replaced by 

their transjugate ® T in the relevant dynamic response formulation. 
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2.9 Shift Synthesis 


The program provides special eigenvalue switching provisions in the analysis to ensure num- 
erical stability. Such a problem may be encountered in the analysis of aerospace structures, which 
are designed to be strong and lightweight. For example, the elements of the mass matrix of equa- 
tion (4) may have numerical values much smaller than those of the stiffness matrix. In such cases, 
the effect of the mass matrix in the (K - X%I)y = 0 formulation may be insignificant. Such a problem 
also occurs in the presence of rigid body modes characterized by "zero" frequencies. An eigenvalue 
shift strategy has been developed to accommodate such situations. 

Thus, the eigenvalue problem pertaining to equation (4) representing the problem defined as 
IPROB = 1 may be written as 


(K-^2M)y = 0 (14) 

in which X is the natural frequency of free vibration, and y is the eigenvector. The stiffness and mass 
matrices must be suitably perturbed to handle rigid body modes and to maintain numerical stability 
by negating effects of rounding error. Thus, equation (14) is rearranged as 

[k + 4 M - (x”+ 4) m] y = 0 (15) 


or 


in which 


(k - Xm) y = 0 

(16) 

K = K + 4M 

✓s 

07) 

M = FM 

(18) 

F 

X = ^- + 4 

(19) 

F 

F - ^ 

(20) 

10 7 

(21) 


where IK^J and IM^J typically denote the norms of the diagonal elements and the number 10 7 relates 
to the computational accuracy of the VAX 1 1 computer. Once the eigenvalue problem defined by 
equation (16) is solved, the natural frequencies are simply obtained as 

X = V(x- 4)f 
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A similar procedure is adopted for the analysis of free vibration problems defined by IPROB = 6 and 
7, as well as for the buckling analysis (IPROB - 9). 

In the case of spinning structures, a somewhat similar strategy is used in perturbing appropriate 
matrices to ensure effective computation of rigid body modes, as well as numerical stability. 


2.10 Formulation for Nodal Centrifugal Forces in Finite Elements 

The STARS program can perform dynamic analyses of structures with nonrotating and rotating 
parts having different spin rates. A general derivation for the in-plane centrifugal forces generated in 
various elements due to the arbitrary spin rate, along with related formulation of the associated 
normal components, is given in detail in reference 6. Reference 7 provides details of a block Lanczos 
algorithm developed for efficient, free vibration analysis of spinning structures. 

Once the nodal centrifugal forces have been derived, as previously mentioned, and stored in array 
P, the element stresses in the structure caused by these forces are simply obtained by solving 
equation (1) (repeated here for convenience), 


KU = P 


The stresses are next utilized to derive the structural geometrical stiffness matrix K G required for 
solving the free vibration problems defined in section 2.7. 


2.11 Material Properties 

The structural material may be general in nature. Thus, the finite element material properties may be 
isotropic, orthotropic, or anisotropic. In the most general case of solid elements having anisotropic 
material properties, defined as material type 3, the stress-strain matrix is expressed as 

5 = Ee (23) 

with Eij being elements of the general material matrix of order 6 by 6, defining the relationship 
between the stress vector 5 and the strain vector e. The elements of the upper symmetric half of the 
E matrix, as well as coefficients of thermal expansion and material density consisting of 28 
coefficients, are the required data input for the pertinent material type. In this connection, it may be 
noted that the material data input is designed in such a way as to be quite general; the user may 
easily incorporate effects of various related features, such as varying material axes orientation, by 
appropriately calculating the elements of the material matrix. If the material is orthotropic, the input 
scheme remains the same for the anisotropic case. 

Material type 2 pertains to thin shell elements displaying anisotropic or orthotropic material 
properties; it requires an input of 13 coefficients. For isotropic material classified as material type 1, 
only four coefficients constitute the required input data. 

The isotropic case for sandwich shell elements is designated as material type 4, whereas type 5 
pertains to the corresponding orthotropic-anisotropic case. 
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2.12 Output of Analysis Results 


A dynamic response analysis, in general, yields an output of nodal deformations and element 
stresses as appropriate functions of time. Additional printouts provide summaries of maximum 
deformations and stresses/loads, as appropriate, as well as principal stresses and relevant angles. 
For line elements, member endloads and moments constitute the usual output of results. In the case 
of thin shell elements, the stresses o xx , o yy and a xy are calculated at the centroid of the element 
and at both its top and bottom surfaces. For solid elements, all six components of stresses (c xx , 

G yy Oz z, cr xy , a yz , a zx ) are computed at the center of the volume of the element. Since free vibration 
analysis constitutes a vital preliminary for the dynamic response analysis, the natural frequencies 
and associated modes are computed by the program and printed out, as desired. Similar results are 
obtained for elastic buckling analysis. For static problems, the nodal displacements and element 
stresses are computed for multiple-load cases. 

Special printout options make possible a selective output of analysis results. Thus, such 
computed data as stiffness and inertia matrices may be printed out, as desired. Initially, the program 
automatically prints out the generated nodal coordinates, element data, and other relevant input 
data. The POSTPROCESSOR program may be effectively used for color graphics depiction of 
solution results. 


2.13 Discussion 

Additional analysis features such as finite, dynamic element discretizations, improved dynamic 
analysis capabilities, and various efficient numerical techniques are continuously being implemented 
in the program. A nonlinear analysis capability is also being developed in parallel. Improved 
preprocessing and postprocessing of data, using E/S PS 390, DEC-VT, CIT, Tektronix, or other 
graphics terminals, are being used to permit efficient modeling and analysis, as well as display, of 
the results pertaining to practical structural problems. 

An automatic data conversion program has also been developed to convert NASTRAN (ref. 8) 
program data into STARS format. 
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3. DATA INPUT PROCEDURE 


(STARS -SOLIDS) 


3.1 Basic Data 


3.1.1 PRIMARY JOB TITLE 
Format (FREE) 

3.1.1. 1 ADDITIONAL JOB DETAILS 
Format (Al, FREE) 

1. Description: Various job-related descriptions, any number of input lines. 

2. Notes: 


First line input is required, and subsequent lines of input must have a C in the first 
column; up to 80 characters per line are accepted. 

3.1.2 NN, NEL, NMAT, NMECN, NEP, NET, NLGCS, NMANGL, NSTACK, MAXLEL 
Format (FREE) 

1. Description: Basic data parameters (structural). 

2. Notes: 


NN 

NEL 

NMAT 


= total number of nodes 

= total number of elements 

= total number of element material types 


NMECN 


= number of material elastic constants, a maximum of numbers, as follows: 
= 4, for isotropic material 

= 13, for orthotropic-anisotropic material for 2-D elements (shell, types 2, 3, 
6, and 7) 

= 10, for isotropic sandwich panel material (shell, types 2 and 3) 

= 25, for orthotropic-anisotropic sandwich panel material (shell, types 2 
and 3) 

= 28, for orthotropic-anisotropic material for 3-D elements (solid, types 4 
and 5) 


NEP 

NET 

NLGCS 


= total number of line element property types (element type 1) 

= total number of shell element thickness types (element types 2 and 3) 
= total number of local-global coordinate systems (LGCS) 


NMANGL = total number of material angle types 


NSTACK = total number of composite shell element stack types 
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MAXLEL = maximum number of layers in a composite shell element 

3.1.3 NTMP, NPR, NSPIN, NC, NBUN, NLSEC, NCNTRL, NOUT, NEXP 

Format (FREE) 

1. Description: Basic data parameters (loads and displacements). 

2. Notes: 

NTMP = total number of element temperature types 

NPR = total number of element uniform pressure types 

NSPIN = total number of different element spin types 

NC = number of sets of nodal loads for IPROB = 8 

= 0, for IPROB = 1 through 7 
= 1, for IPROB = 9 

NBUN = total number of interdependent and finite nodal displacement connectivity 
conditions (includes IDBC and FDBC in section 2.3, being equal to number 
of lines of input) 

NLSEC = total number of line element special end conditions excluding commonly 
occurring cases of purely rigid or hinged ends 

NCNTRL = total number of control surface rigid body modes used for ASE analyses; 
may also be utilized for generating perfect rigid body modes 

NOUT = total number of output nodes where direct modal interpolation is effected; to 
be set to 0 for alternative interpolation scheme effected by GRIDCHG 
program 

NEXP = total number of uniform external in-plane pressures for membranes 

3.1.4 IPROB, IEIG, IDRS, IBAN, IPLUMP, IMLUMP, INMM, IINTP 

Format (FREE) 

1. Description: Data defining nature of required solution. 

2. Notes: 

IPROB = index for problem type, to be set as follows: 

= 1, undamped, free vibration analysis of nonspinning structures 
= 2, undamped, free vibration analysis of spinning structures 
= 3, quadratic matrix eigenproblem option for DEM (dynamic element 
method) analysis 
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IEIG 


= 4, free vibration analysis of spinning structures with diagonal viscous 
damping matrix 

= 5, as for IPROB = 4 with structural damping 

= 6, free vibration analysis of nonspinning structures with general viscous 
damping 

= 7, as for IPROB = 6 with structural damping 

= 8, static analysis of structures with thermal and multiple mechanical load 
cases 

= 9, elastic buckling analysis 

= integer defining eigenproblem solution type 

= 0, for solution based on a modified, combined Sturm sequence and inverse 
iteration method 

= 1, for an alternative solution technique based on a block Lanczos procedure 
(recommended for computation of first few roots and vectors when the 
lower bound PL = 0 for cases IPROB = 1, 2, 3, and 9) 


IDRS 


IBAN 


= index for dynamic response analysis 
= 0, no response analysis required 
= 1, performs response analysis 

= bandwidth minimization option 
= 0, performs minimization 
= 1, minimization not required 
= -1, option to perform minimization only and exit 


IPLUMP = index for nodal external loads 
= 0, no load input 

= 1, concentrated nodal load input for IPROB = 8 and 9, as well as for 
IPROB = 1 through 7 for prestressed structures 


IMLUMP = index for nodal lumped scalar mass 
= 0, no lumped mass 

= 1, lumped nodal mass input (IPROB = 1 through 7) 


INMM = index for nodal 6 by 6 mass matrix 
= 0, no mass matrix 

= 1, nodal mass matrix input (IPROB = 1 through 7) 


IINTP 


= integer defining modal data for direct interpolation 
= 0, no interpolation required 

= 1, performs interpolation on STARS calculated modal data 
= 2, performs interpolation on externally supplied modal data; for example, 
GVS results 


3. Notes: 


A dynamic response analysis is achieved by specifying appropriate values for IPROB and 
IDRS; at end of problem solution, extensive options are available for plotting nodal defor- 
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mations, mode shapes, and element stresses by utilizing the postprocessor program 
POSTPLOT. 

Initial static load (prestress) effect: in the case of dynamic problems, the presence of 
nonzero values of integers IPLUMP, NPR, and/or NTMP activates computation of 
prestressing effect. 

3.1.5 IPREC, IPLOT, XPRINT, INDATA, IERCHK 
Format (FREE) 

1. Description: Additional basic data. 

2. Notes: 

IPREC - specification for solution precision 
= 1, single precision 
= 2, double precision 

IPLOT = index for graphics display 
= 0, no plotting needed 

= 1, performs display of input geometry; if satisfactory, a restart option 
enables continuation of current analysis 

IPRINT = output print option 

= 0, prints final results output only 

= 1, prints global stiffness (K), mass (M), damping or Coriolis (C) matrices, 
as well as detailed output on deformations, stresses, and root convergence 
characteristics 

= 2, prints output as in IPRINT = 1, but omits K, M, and C matrices 
= 3, output as in IPRINT = 0, but omits eigenvector printouts 

INDATA = input data option 

= 0, basic matrices are automatically computed 

= 1, to read upper symmetric banded half of basic matrices K, M, and C from 
user input files, row-wise 

IERCHK = integer defining level of input data-error checks required by user 
= 0, usual level of error checkouts 
= 1, additional extensive data checkouts 

Mass matrix: nodal lumped mass matrix is added to consistent mass matrix to evolve the final 
mass matrix. 

3.1.6 INDEX, NR, INORM, PU, PL, TOL (Required if IPROB * 8) 

Format (FREE) 

1. Description: Data specifications for eigenproblem solution. 
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2. Notes: 


INDEX 


NR 

INORM 


= indicator for number of eigenvalues and vectors to be computed 
= 1, computes NR smallest roots (and vectors) lying within bounds PU, PL 
= 2, computes all roots (and vectors) lying within bounds PU, PL 

= number of roots to be computed (any arbitrary root number input for 
INDEX = 2) 

= index for vector normalization; any desired vector row number 
= 0, normalizes with respect to a scalar of displacement vector Y having 
largest modulus 

= -1, normalizes with respect to a scalar of Y or YD (velocity) vector having 
largest modulus 


PU = upper bound of roots 

PL = lower bound of roots 


TOL = tolerance factor (eq. (21)) 

= 0, defaults to 25. 0E + 08 

~ X, defaults to X (X = 1.0E + 07 may be useful for computation) 


3.1.7 IUV, IDDI, NTTS, NDELT (Required if IDRS = 1) 

Format (FREE) 

1. Description: Data related to dynamic response analysis. 

2. Notes: 


3.1.8 


IUV = index for initial displacement (U) and velocity (V) input 
= 0, no initial data 

= 1, either initial displacement or velocity or both are nonzero vectors 


IDDI 


NTTS 


= index for dynamic data input 
= 1, nodal load input 
= 2, nodal acceleration input 

= total number of sets of load or acceleration data input 


NDELT = number of sets of uniform time increments for response calculation 


G (Required if IPROB = 5 or 7) 

Format (FREE) 


1. Description: Structural damping in formulation [K = K(1 + i*G)]. 
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2. Notes: 


G = structural damping parameter 

i* = imaginary number, f~i 

K = system stiffness matrix 

3.1.9 Mil 

Format (FREE) 

1. Description: Half- bandwidth of K, M, or C. 

3.1.10 ( ( B(I,J), I = 1, N ), J = 1, NC ) 

Format (6E10.4) 

1. Description: Load matrix of order N = NN*6. 

3.1.11 ( ( K(I,J), 1=1, Mil ), J = 1, N ) 

Format (6E10.4) 

1. Description: Stiffness matrix. 


(Required if INDATA - 1) 


(Required if INDATA = 1 and IPROB = 8) 


(Required if INDATA = 1 
and IPROB = 1 through 8) 


3.1.12 ( ( M(I,J), J = 1, Mil ),I = 1,N) 
Format (6E10.4) 

1. Description: Mass matrix. 

3.1.13 ( ( C(I,J), J = 1, Mil ), I = 1, N ) 
Format (6E10.4) 


(Required if INDATA = 1 
and IPROB = 1 through 7) 


(Required if INDATA = 1 
and IPROB = 2 through 5) 


1. Description: Coriolis (IPROB = 2, 4, 5) or dynamic correction (IPROB = 3) matrix. 

3.1.14 ( ( CD(I,J), J = 1, Mil ), I = 1, N ) (Required if INDATA = 1 

Format (6E10.4) and IPROB = 4 through 7) 

1. Description: Viscous damping matrix. 


2. General note: 


Each set of data input in succeeding sections is preceded with a relevant comment 
statement having a dollar sign ($) at the first column, followed by optional descriptive 
words. 


3. Note: 


If INDATA = 1, no further input is required. 
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3.2 Nodal Data 


3.2.1 $ NODAL DATA 

3.2.2 IN, X, Y, Z, UX, UY, UZ, UXR, UYR, UZR, ILGCS, IZDRCS, HNC 
Format (I5,3E10.4,9I5) 

1. Description: NN sets of nodal data input in GCS/LGCS, at random; table 1 provides a 

description of the input data. 

Table 1. Arrangement of nodal data input. 


Node 

Nodal 

Nodal zero displacement 

Local-global 

coordinate 

ZDBC reference 

Increment 

number 

coordinates 

boundary conditions (ZDBC) 

system type 

coordinate system 

(IN) 

(X) (Y) (Z) 

(UX) (UY) (UZ) (UXR) (UYR) (UZR) 

(ILGCS) 

(IZDRCS) 

(UNC) 

* 

* * * 

1 2 3 4 5 6 

* * * * * * 

* 

* 

* 


2. Notes: 

a. A right-handed Cartesian coordinate system (X, Y, Z) is to be chosen to define the 
global coordinate system (GCS). 

b. The asterisk (*) indicates required data input in GCS/LGCS. 

c. Each structural node is assumed to have six degrees of freedom (DOF) consisting of 
three translations, UX, UY, UZ, and three rotations, UXR, UYR, UZR, usually labeled 
as displacement degrees of freedom 1, 2, 3, and 4, 5, 6, respectively. 

d. For nodal zero displacement boundary conditions (ZDBC) defined in coordinate 
system referred to as IZDRCS, set value to 

= 0, for free motion, 

= 1, for constrained motion. 

e. For node generation by increment, set IINC 
ss 0, for no increment, 

— I, to increment node number of previous input by I until current node number is 
attained; coordinates of intermediate nodes are linearly interpolated. 

f. In automatic node generation (note (e)), all relevant data of generated intermediate 
nodes pertain to that of the last data set of the sequence. 

g. Third-point nodes for line elements are assumed to lie on element local x-y plane and 
may be chosen as any existing active node or dummy nodes with UX through UZR set 
to 1. 

h. Final data are automatically formed in increasing sequence of node numbers. 
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3. Notes 


ILGCS = integer specifying local-global coordinate system number (set to 0 if data is 
in GCS), defining nodal data 

IZDRCS = integer defining zero displacement boundary condition reference coordinate 
system (set to 0 for data in GCS or an ILGCS number) 

3.2.3 $ LOCAL-GLOBAL COORDINATE SYSTEM DATA (Required if NLGCS * 0) 

3.2.4 ILGCS, IDMOD 
Format (215) 

3.2.5 XOR, YOR, ZOR, X2, Y2, Z2, 

X3, Y3, Z3 
or 

XOR, YOR, ZOR, Dll, D12, D13, 

D21, D22, D23, D31, D32, D33 
Format (2(6E10.4,/)) 

1. Description: NLGCS sets of local-global coordinate system (LGCS) definition data, 

at random. 

2. Notes: 


(IDMOD = 1) 
(IDMOD - 2) 


IDMOD = integer specifying nature of input data 

= 1, input involves global coordinates of the origin of the LGCS (XOR, YOR, 
ZOR) and two data points (X2 through Z3, pertaining to two points located 
on LGCS X-axis and X-Y plane, respectively) in GCS 
= 2, involves input of origin of LGCS (XOR, YOR, ZOR) and elements of 
direction cosine matrix of the LGCS 

3. Special note: 

If IINTP = 2, no further data input is required until 3.5.7. 
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3.3. Element Data 


General note: Element data input may be at random within each data group. 

3.3.1 $ ELEMENT CONNECTIVITY 

3.3.2 IET, IEN, ND1, ND2, ND3, ND4, ND5, ND6, ND7, ND8, IMPP, IEPP/ITHTH, ITMPP, 
IPRR, 1ST, INC 

Format (1615) 

1. Description: NEL sets of element data input; definition of input data is given in table 2. 


Table 2. Element data layout. 


Element 

type 

(IET) 

Element 

number 

(IEN) 

Node number for vertices 

1 2 3 4 5 6 7 8 

(ND1) (ND2) (ND3) (ND4) (ND5) (ND6) (ND7) (ND8) 

IMPP 

IEPP/ 

ITHTH 

ITMPP 

IPRR 

1ST 

INC 

Line (bars, 
rods, beams, 
3-D lines) 

1 

* 

* 

* 

** 

IEC1 

IEC2 

V 



* 

X 

* 

* 

* 

* 

Shell quadri- 
lateral (plane, 
plate, shear, 
shell - usual and 
sandwich) 

2 

* 

* 

* 

* 

* 




A 

* 

— 1_ 

* 

* 

* 

* 

Shell triangular 
(plane, plate 
shear, shell - usual 
and sandwich) 

3 

* 

* 

* 

* 





A 

* 

t 

* 

* 

* 

* 

Solid 

hexahedron 

4 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 

* 

Solid 

tetrahedron 

5 

* 

* 

* 

* 

* 





* 


* 

* 

* 

* 

Shell quadri- 
lateral 

composite element 
(plane, plate 
shear, and shell) 

6 

* 

* 

* 

* 

* 




□ 

* 


* 

* 

* 

* 

Shell triangular 
composite element 
(plane, plate 
shear, and shell) 

7 

* 

* 

* 

* 





□ 

* 


* 

* 

* 

* 

Rigid element 

(pin-ended bar, 
rigid body) 

8 

* 

* 

* 

** 

IEC1 

IEC2 

o 








* 

Prestressed 

rectangular 

membrane 

1 1 

* 



* 

* 

* 

* 




9 

* 



t 

* 

* 

* 

* 
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* 


** 

IECI 


IMPP 

IEPP(x) 


=s data as defined; under element type 8, individual rigid elements are 
characterized by appropriate data entry 

= third point node for element types 1 and 8 

= integer defining line element end condition pertaining to end I 
= 0, rigid-ended 

= 1, pin-ended in three rotational degrees of freedom 

= J, denoting special end condition number, to be set greater than 1; for scalar 
springs, set IECI to a negative value less than -1 

= integer defining material property type number 

= integer defining line element property type 


ITHTH(t) = integer defining shell element thickness type 


ITMPP = integer defining element temperature type 


IPRR = integer defining element pressure type 


1ST = integer defining element spin type 


INC 


V 

A 




= integer for element generation by increment 

- 0, no increment 

= J, increments' node numbers of previous elements by J until current element 
nodal numbers are reached 

= ILGCS, integer defining LGCS associated with a zero-length scalar spring 
element; defaults to GCS 

= IMANG, integer defining material angle type number, suitable for layered 
elements 

= ISTACK, stack type number, used for integrated composite elements 
(types 6 and 7) 

= dependent degree of freedom at ND1 to be rigidly connected to all six 
degrees of freedom at ND2; rules concerning interdependence of nodes and 
degrees of freedom are defined in section 3.4.4 

= -1, if all six degrees of freedom at ND1 are involved 

= 0, for pin-ended rigid bar elements 


0 


= integer defining prestress type 


Rigid elements may be specified to span any length, including 0. 



Rigid pin-ended bar elements may be simulated by setting EEC1 = EEC2 = 1. 

In automatic element generation (see INC, above), the generated intermediate elements 
acquire properties the same as the last element in current sequence. Also, a special 
option enables repetitive use of an element with an input format (13, 12, 1515); the integer 
IET is then replaced by NELNO and IET, where NELNO is the total number of similar 
elements connecting the specified nodes. 

Sandwich shell elements may be generated by individual inputs of membrane, bending, 
and transverse shear effects. Furthermore, the composite shell elements consisting of 
layered composites can be formed for varying stacks of materials. 

For element type 8, defined by two nodes, if the first node has some ZDBC constraints, 
the latter should also be applied to the second node. 

3. Element description: 

The various elements (fig. 6) and associated degrees of freedom are depicted as follows. 
The global coordinate system (GCS) is represented by X, Y, Z, whereas x, y, z relates to 
local coordinate system (LCS). 


4. Notes: 

a. A right-handed Cartesian coordinate system (x, y, z) is to be chosen to define any 
element local coordinate system (LCS). 

b. Any node may be chosen as the first vertex of an element, the local x-axis being along 
the line connecting vertices 1 and 2. 

c. For line elements, the local x-y plane is defined as the plane contained by vertices 1 , 

2, and the specified third-point node. 

d. The vertices of thin shell elements are usually numbered in a counterclockwise 
sequence when observed from any point along the local positive z-axis, being also 
utilized as plane and plate-bending elements, as appropriate. 

e. For solid elements, the y-axis lies in the plane formed by vertices 1-2-3 and 1 -2-3-4 for 
the tetrahedral and hexahedral elements, respectively; the z-axis is perpendicular to 
the x-y-plane, heading toward the fourth node for the tetrahedron and the plane 
containing the other four nodes for the hexahedral element. 

f. The vertices of the solid elements are also numbered in a counterclockwise sequence 
when viewed from any point on the positive z-axis lying above the plane under 
consideration; the fifth vertex of the hexahedron is to be chosen as the node directly 
above vertex 1. 

g. For layered composite shell element types 6 and 7, the layering sequence starts with 
the layer that has maximum ~z coordinate expressed in element LCS. 
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(a) Line element. (b) Quadrilateral shell element. 



(c) Triangular shell element. (d) Hexahedral solid element. 



2 10135-5 


(e) Tetrahedral solid element. 
Figure 6. STARS-SOLIDS element types. 
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5. Structural modeling: 


Since each node is assumed to possess six displacement degrees of freedom, any 
individual structural form may be simply represented by suppressing appropriate 
displacement terms. The following rules may be adopted: 

Truss structures: to allow only two nodal translational deformations in the plane 
of the structure; to use line elements. 

Plane frame: all three in-plane displacements, namely, two translations and one 
rotation, are retained in the formulation; to use line elements. 

Plane stress/strain: displacement boundary conditions are similar to truss 
structures; to use shell elements. 

Plate bending: only the three out-of-plane displacements consisting of one 
translation and two rotations are considered for the analysis; to use shell 
elements. 

Solid structures: the three translational degrees of freedom are retained in the 
analysis; to use solid elements. 

Shell, space frame: all six degrees of freedom are to be retained in the solution 
process; to use shell and line elements, respectively. 

Suppression of derived nodal motion may be achieved by using zero and interdependent 
displacement boundary conditions (ZBDC, IDBC) defined in sections 3.2 and 3.4, 
respectively. 

3.3.3 $ COMPOSITE SHELL ELEMENT STACK DESCRIPTION DATA (Required for 

composite shell elements (types 6 and 7), and only if NSTACK ^ 0) 

3.3.4 ISTACK, NLAYER 
Format (215) 

3.3.5 (IMATC(I), THCL(I), IMANGC(I), 1=1, NLAYER) 

Format (15, El 0.4, 15) 

1. Description: NSTACK sets of composite shell element data. 

2. Notes: 


ISTACK = 
NLAYER = 
IMATC(I) = 


stack number 

number of layers in the stack 

material type number for the composite layer 
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THCL(I) - thickness of the composite layer 

IMANGC(I) = integer specifying material angle type number (IMANG) 

3 3 6 $ SPECIFICATION FOR MATERIAL AXES ORIENTATION (Required if 

NMANGL * 0) 

3.3.7 IMANG, IMAMD, ILGCS 
Format (315) 

3.3.8 Dll, D12, D13, D21, D22, D23, (IMAMD = 1) 

D31.D32, D33 

Format (2(6E10.4,/)) 
or 

THETA (IMAMD = 2) 

Format (E10.4) 

1. Description: NMANGL sets of material angle definition data. 

2. Notes: 


3.3.9 


IMAMD = integer defining material angle data input mode 

= 1, involves input of elements of direction cosine matrix of material axes 
with respect to LGCS/GCS (set ILGCS = 0 for data in GCS) 

= 2, requires input of material axis angle (THETA) with shell element local 
x-axis 


ILGCS 

THETA 


= integer specifying local-global coordinate system number (set to 0 if data is 
in GCS) 

= material axis angle with respect to shell element local x-axis 


$ LINE ELEMENT BASIC PROPERTIES (Required for line elements only) 


3.3.10 IEPP, A, JX, IY, IZ, SFY, SFZ 
Format (15, 6E10.4) 


1. Description: 

2. Notes: 


NEP sets of line element basic property data in element local coordinate system 
(LCS). 


IEPP = integer denoting line element property type 
A = area of cross section 

JX = torsional moment of inertia about element x-axis 

IY = moment of inertia about element y-axis 
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IZ 


= moment of inertia about element z-axis 


SFY = A/ASY, shear area (ASY) factor along y-axis 

SFZ = A/ASZ, shear area (ASZ) factor along z-axis 

For no shear area effect, SFY and SFZ are to be set at 0.0. 

3.3.1 1 $ LINE ELEMENT SPECIAL END CONDITIONS (Required for line elements only 

if NLSEC * 0) 

3.3.12 ILSEC, (k(I), I = 1, 6) 

Format (15, 6E10.4) 

1. Description: NLSEC sets of line element special end conditions data in LCS. 

2. Notes: 

ILSEC = element end condition type (to be set greater than 1), referring to members 
attached to the nodes at its ends by flexible connections, or members with 
free end degrees of freedom in LCS (corresponds to IEC1 and IEC2) 

= set to a negative value, less than -1, for scalar springs connecting two 
nodes (corresponds to IEC1) 

k(I) = additional spring stiffness along Ith translational (x-, y-, and z-direction) 
degree of freedom and actual rotational Ith spring stiffness (x, y, and z 
rotational constraint) 

= -2, for rigid rotational Ith constraint 

- ~1, for release of corresponding member end degree of freedom, relevant 
also to ILSEC value set greater than 1 

= stiffness values for scalar springs associated with a negative ILSEC value 
less than -2 

Such elements may have 0 or any finite length. 

To simulate only specified end condition, set Young's modulus E = 0 for the corresponding 
material type, IMPP. 

3.3.13 $ SHELL ELEMENT THICKNESS (Required for shell elements 

(types 2 and 3) only) 

3.3.14 ITHTH, TM, TB, TS 
Format (15, E10.4) 

1. Description: NET sets of element thickness data. 

2. Notes: 
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ITHTH = element thickness type 

TM = membrane element thickness 

TB = bending element thickness 

TS = transverse shear element thickness 

Above shell thickness pertains to sandwich elements; in the absence of data for TB and 
TS, the shell element thickness T is taken as TM. 

For consistent mass matrix formulation, shell thickness T is taken as TM. 

33.15 $ ELEMENT MATERIAL PROPERTIES 

3.3.16 IMPP, MT 
Format (215) 

3.3.17 E, MU, ALP, RHO 

Ell, E12, E14, E22, E24, E44, E55, 

E56, E66, ALPX, ALPY, ALPXY, RHO 

Ell, E12, E13, E14, E15, E16, E22, 

E23, E24, E25, E26, E33, E34, E35, 

E36, E44, E45, E46, E55, E56, E66, 

ALP1, ALP2, ALP3, ALP4, ALP5, ALP6, RHO (material type 3); or 

EM, EB, ES, MUM, MUB, MUS, 

ALPM, ALPB, ALPS, RHO (material type 4); or 

E11M, E12M, E14M, E22M, E24M, E44M, 

El IB, E12B, E14B, E22B, E24B, E44B, 

E55S, E56S, E66S, ALPXM, ALPYM, ALPXYM, 

ALPXB, ALPYB, ALPXYB, ALPXS, ALPYS, ALPXYS, RHO (material type 5) 

Format (4(7E10.4,/» 


(material type 1); or 
(material type 2); or 


1. Description: NMAT sets of element material property data; the individual material matrices 

are derived from the 6 by 6 symmetric matrix for general solid material. 


2. Notes: 


IMPP = material number 


MT = material type 

= 1, isotropic 

= 2, orthotropic-ani so tropic, shell elements 
= 3, orthotropic-anisotropic, solid elements 
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E 

E ij 

MU 

ALP 


= 4, isotropic, sandwich shell elements incorporating individual membrane, 
bending, and transverse shear effects 

= 5, orthotropic-anisotropic sandwich shell elements with individual effects, 
as above 

= Young's modulus 

= elements of material stress-strain matrix (i = 1, 6 ; j = 1, 6) 

= Poisson's ratio 

- coefficient of thermal expansion for isotropic material 


ALPX, ALPY, 

ALPXY = coefficients of thermal expansion, shell elements 


ALP1 through 

ALP6 = coefficients of thermal expansion, solid elements 
RHO = mass per unit volume 

For sandwich elements (material types 4 and 5), relevant notations defining such 
properties utilize a postscript of M, B, or S for membrane, bending, or transverse shear 
stiffness, respectively. 

3.3.18 $ ELEMENT TEMPERATURE DATA (Required if NTMP * 0) 

3.3.19 ITMPP, T, DTDY, DTDZ 
Format (2(I5,3E10.4)) 


1. Description: NTMP number of element temperature types; table 3 shows compatible input 

data. 


Table 3. Element temperature data input. 


Element 

type 

T 

DTDY 

DTDZ 

1 

* 

* 

* 

2,3, 6, 7 

* 


* 

4,5 

* 




2. Notes: 


ITMPP = element temperature increase type 
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T 


= uniform temperature increase; relates to all elements 


DTDY = temperature gradient along element local y-axis; relates to line dements 
only 

DTDZ = temperature gradient along element local z-axis; relates to line and shell 
elements 

* = compatible input data 

3.3.20 $ ELEMENT PRESSURE DATA (Required if NPR * 0) 

3.3.21 IPRR, PR 

Format (5(I5,E10.4)) 

1. Description: NPR sets of element pressure data. 

2. Notes: 

IPRR = element pressure type 

PR = uniform pressure 

Pressure directions for line elements: uniform pressure is allowed in local y- and z- 
direction only, and the program calculates as input both end loads and moments; while 
pressure corresponding to a first nodal input pertains to y-direction, a subsequent input for 
the same node signifies pressure acting in the z-direction. 

Pressure directions for shell elements: uniform pressure is allowed in local z-direction 
only; the program computes nodal load input. 

Pressure directions for solid elements: uniform pressure is allowed on base surfaces 
defined by nodes 1-2-3-4 and 1-2-3 for hexahedral and tetrahedral elements, respectively, 
acting in local z-direction; the program computes nodal load input data. 

3.3.22 $ PRESTRESSED RECTANGULAR MEMBRANE ELEMENT DATA 

(Required if NEXP * 0) 

3.3.23 IIEXP, SX, SY 
Format (I5,2E10.4) 

1 . Description: NEXP sets of prestressed membrane stress data. 

2. Notes: 

IIEXP = integer defining stress combination type 

SX, SY = membrane stresses in the element x- and y-directions, respectively 
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3.4 Data in Global or Local-Global Coordinate System 
General note: Data input may be at random within each data group. 

3.4.1 $ ELEMENT SPIN RATE DATA (Required if NS PIN * 0) 

3.4.2 1ST, SPX, SPY, SPZ, ELGCS 
Format (15, 3E10.4, 15) 

1. Description: NSPIN sets of spin data. 

2. Notes: 


1ST = spin type 

SPX, SPY, SPZ = components of element spin rate in global/local- global X-, Y-, and 
Z-directions, respectively 

ILGCS = local-global coordinate system number, as defined in section 3.2.2 

3.4.3 $ DISPLACEMENT BOUNDARY CONDITION DATA (Required if NBUN * 0) 

3.4.4 INI, IDOFJ, INIP, IDOFJP, CONFCT, IDRCS, NDBCON 
Format (415, E10.4, 215) 

1. Description: NBUN sets of nodal interdependent displacement boundary condition (IDBC) 

data. 


2. Notes: 


INI 

IDOFJ 

INIP 

IDOFJP 

CONFCT 

IDRCS 

NDBCON 


node number I 

Jth DOF associated with node I 
node number I' 

J'th DOF associated with node I' 
connectivity factor 

displacement boundary condition reference coordinate system 

integer defining displacement boundary condition increment 
0, no increment 

an integer, to increment IDOFJ and IDOFJP by 1 until IDOFJ reaches 
NDBCON value 
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J and J' vary between 1 and 6. 

3. Additional notes: 

The nodal displacement boundary conditions’ relationship is expressed as 
Ui j = am,n Um,n 

= a i,jUij+ ai' t j'Ui',j , + ... 

The input scheme is shown in table 4. 


Table 4. Data layout for displacement boundary conditions. 


Node 

1 

DOF 

Node 

2 

DOF 

Connectivity 

coefficient 

Terminology 

Reference 
coordinate system 

Incremental 
DOF value 

i 

j 

i' 

j' 

a i\j' 

IDBC 

IDRCS 

NDBCON 

i 

j 

i 

j 

aj 

FDBC 

IDRCS 

NDBCON 

i 

J — 

i 

. j - 

0 

ZDBC 

IDRCS 

NDBCON 


in which 


i, i' 

= node numbers, 

JJ 

= degrees of freedom, 

a i,j> a i',j‘ 

= connectivity coefficients. 


IDBC, FDBC, and ZDBC are, respectively, the interdependent, finite, and zero 
displacement boundary conditions. The ZDBC may also be conveniently implemented by 
following the rules given in table 1, which is generally recommended for such cases. It 
should be noted that the dependent degrees of freedom appearing in columns 1 and 2 may 
not appear subsequently in columns 3 and 4 as independent degrees of freedom. 

However, the independent degrees of freedom may be subsequently related. 

3.4.5 $ NODAL LOAD DATA (Required if IPLUMP * 0) 

3.4.6 IN, IDOF, P, IDOFE 
Format (215, E10.4, 15) 

1. Description: NC sets of nodal force data. 

2. Notes: 

IN = node number 
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IDOF and XDOFE are, respectively, the start and end degrees of freedom assigned with 
the same P value; default value for IDOFE is IDOF. 

P = nodal load 

Each data set is to be terminated by setting a negative value for IN. 

3.4.7 $ NODAL MASS DATA (Required if IMLUMP * 0) 

3.4.8 IN, IDOF, M, IDOFE 
Format (215, E10.4, 15) 

1. Description: Nodal lumped mass data. 

2. Notes: 

M = nodal mass 

Other definitions are as in section 3.4.6. 

3.4.9 $ NODAL MASS MATRIX IN LGCS/GCS (Required if INMM * 0) 

3.4.10 IN, ELGCS 
Format (215) 

3.4.11 (VNMDAT(I), I = 1, 36) 

Format (6(6E10.4,/)) 

1. Description: User input of 6 by 6 nodal mass matrix. 

2. Notes: 

The user may input data for only the upper symmetric elements; numbers in lower half 
may be set to zero as the program automatically symmetrizes the matrix. 

For data in GCS, set LGCS = 0. 

Each data set is to be terminated by setting a negative value for IN. 

3.4. 12 $ NODAL INITIAL DISPLACEMENT AND VELOCITY DATA (Required if IUV = 1 

and IDRS = 1 ) 

3.4.13 IN, IDOF, UI, VI 
Format (215, 2E15.5) 

1. Description: Initial displacements and velocities data. 

2. Notes: 
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IN 


= node number 


ID OF = degree of freedom 

UI = initial displacement value 

VI = initial velocity value 

Data set is terminated if IN is read as -1. 

3.4.14 $ NODAL FORCE ACCELERATION DATA (Required if NTTS * 0 and IDRS = 1) 

3.4.15 TZ 

Format (E15.5) 

3.4.16 IN, IDOF, PZ 
Format (215, E15.5) 

1. Description: NTTS sets of dynamic nodal load (IDDI = 1) or acceleration (IDDI = 2) input 

data. 

2. Notes: 

TZ = time-duration of load application 

PZ = nodal force or acceleration data 

Each data set is terminated by setting IN value to -1; other definitions are as given in 
section 3.4.6. 

3.4.17 $ INCREMENTAL TIME DATA FOR RESPONSE CALCULATION (Required if 

NDELT * 0 and IDRS = 1) 

3.4.18 DELT, IDELT 
Format (E15.5, 15) 

1. Description: NDELT sets of uniform incremental time input data for dynamic response 

calculations. 

2. Notes: 

DELT = uniform incremental time step 

IDELT = total number of uniform time steps in the data set 

3.5 Additional Basic Data 

3.5. 1 $ VISCOUS DAMPING DATA (Required if IPROB = 4 or 5) 
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3.5.2 (C(I,1),I= 1, N ) 

Format (6E10.4) 

1. Description: User input of diagonal viscous damping matrix. 

2. Notes: 

C = diagonal viscous damping matrix 

N = order of matrix 

3.5.3 $ COEFFICIENTS FOR PROPORTIONAL VISCOUS DAMPING (Required if 

IPROB = 6 or 7) 

3.5.4 ALPHA, BETA 
Format (2E10.4) 

1. Description: Proportional viscous damping formulation C = ALPHA*K + BETA*M 

2. Notes: 

ALPHA and BETA are damping parameters. 

K and M are system stiffness and mass matrices. 

3.5.5 $ USER INPUT OPTION FOR VISCOUS DAMPING MATRIX (Required if IPROB = 6 

or 7 and ALPHA and BETA set to 0) 

3.5.6 ( ( C(I,J), J = 1, Mil), 1=1, 6) 

Format (6E10.4) 

1. Description: NN sets of user input of banded viscous damping matrix C(N,M1 1) in blocks of 

six rows of bandwidth Mil, one row at a time (N = 6*NN). 

2. Notes: 

Data file must conform to IDBC, FDBC, and ZDBC, inherent in the problem. 

3.5.7 $ MEASURED MODAL DATA INPUT (Required if IINTP = 2) 

3.5.8 ( INODM(I), (DISPLM(I,J), J = 1,6), I = 1, NN ) 

Format (15, 6E10.4) 

1. Description: Measured modal displacement data input, NR sets of data. 

2. Notes: 


Each data set to be terminated by setting INODM(I) value to -1. 
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3.5.9 $ OUTPUT POINTS SPECIFICATION FOR DIRECT INTERPOLATION OF MODAL 

DATA (Required if NOUT * 0) 

3.5.10 ( IOUTP(I), ( ICONP(IJ), J = 1, 6 ), I = 1, NOUT) 

Format (715) 

1. Description: To read output point and up to six connecting points 

2. Notes: 

IOUTP(I) = output points on AERO interpolation lines 

ICONP(I,J) = STARS finite element nodes whose deflections will be averaged to 
calculate the deflection value at the interpolation point 

3.5. 1 1 $ RIGID CONTROL MODES DATA INPUT (Required if NCNTRL * 0) 

3.5.12 INS, IDOF, DISP, INE, ININC 
Format (215, E10.4, 215) 

1. Description: Modal displacement data for NCNTRL number of modes. 

2. Notes: 

INS, INE = starting and end node numbers; default value for INE is INS 
IDOF = degree of freedom, a value between 1 and 6 
DISP = associated displacement 

ININC = integer defining nodal incremental value; to increment INS by ININC until 
INE is attained 

Each data set is to be terminated by setting INS value to —1. 
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4. SAMPLE PROBLEMS 


(STARS-SOLIDS) 


This section provides the input data, as well as relevant outputs, of 12 typical test cases 
involving static, stability, free vibration, and dynamic response analyses of representative 
structures. The input data are prepared in accordance with the procedures described in section 3 . 
Details of such analyses are in the descriptions that follow in which each structural geometry is 
described in a right-handed, rectangular coordinate system, and the associated input data are defined 
in consistent unit form. 


4.1 Space Truss: Static Analysis 

The static analysis of the space truss depicted in figure 7 (ref. 9) was performed to yield nodal 
deformations and element forces. A load of 300 lb acts at node 7 along the axial direction of the 
member connecting nodes 7 and 9; another load of 500 lb is applied at node 10 in the direction of the 
structural base centerline. Also, the three members in the upper tier of the structure are subjected to 
a uniform temperature increase of 100°. Two rigid elements are, however, introduced between nodes 
5 and 8 and nodes 7 and 9. 


z 



Important data parameters: 

Young's modulus, E 
Poisson's ratio, p 
Coefficient of 

thermal expansion, a 


1.0 x 10 7 
0.3 

12.5 x 10- 6 
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STARS input data; 


SPACE TRUSS - MECHANICAL AND THEFfAL LEADING - REARS FRCM NCDES 5-8, 7-9 

11.21.1.4.1.0. 0.0. 0.0 

1 . 0 . 0 . 1 . 0 . 0 . 0 . 0.0 

8 , 0 , 0 , 1 , 1 , 0 , 0,0 

2 , 0 , 1 , 0,1 

$ NODAL DATA 


1 


0.0 

72.0 


0.0 

1 

1 

1 

0 

0 

0 

2 


0.0 

-72.0 


0.0 

1 

1 

1 

0 

0 

0 

3 

124 

,.68 

0.0 


0.0 

1 

1 

1 

0 

0 

0 

4 

13 

1.86 

48.0 


72.0 







5 

13.86 

-48.0 


72.0 







6 

96.972 

0.0 


72.0 







7 

27. 

708 

24.0 

144.0 







8 

27. 

708 

-24.0 

144.0 







9 

69.276 

0.0 

144.0 







10 

41. 

,568 

0.0 

216.0 







11 

144.0 

36.0 


0.0 

1 

1 

1 

1 

1 

1 

ELEMENT CONNECTIVITY 









1 

1 

1 

4 11 

1 

1 

0 

0 

0 

1 

1 


1 

2 

2 

4 11 

1 

1 




1 

1 


1 

3 

2 

5 11 

1 

1 




1 

1 


1 

4 

3 

5 11 

1 

1 




1 

1 


1 

5 

3 

6 11 

1 

1 




1 

1 


1 

6 

3 

4 11 

1 

1 




1 

1 


1 

7 

4 

5 11 

1 

1 




1 

1 


1 

8 

5 

6 11 

1 

1 




1 

1 


1 

9 

6 

4 11 

1 

1 




1 

1 


1 

10 

4 

7 11 

1 

1 




1 

1 


1 

11 

5 

7 11 

1 

1 




1 

1 


8 

12 

5 

8 11 



-1 






1 

13 

6 

8 11 

1 

1 




1 

1 


1 

14 

6 

9 11 

1 

1 




1 

1 


1 

15 

6 

7 11 

1 

1 




1 

1 


1 

16 

7 

8 11 

1 

1 




1 

1 


1 

17 

8 

9 11 

1 

1 




1 

1 


8 

18 

9 

7 11 

1 

1 







1 

19 

7 

10 11 

1 

1 




1 

1 

1 

1 

20 

8 

10 11 

1 

1 




1 

1 

1 

1 

21 

9 

10 11 

1 

1 




1 

1 

1 


$ LINE ELEMENT BASIC PROPERTIES 
1 0.01389 


$ ELUENT NftTERIAL PROPERTIES 
1 1 

10.0E6 0.3 12.5E-06 

$ ELEMENT TEMPERATURE DATA 
1 100.0 
$ NCDAL LOAD DATA 
10 1 -500.0 

7 1 -259.8 

7 2 150.0 

-1 


STARS analysis results - nodal deformations and element stresses: 


LOAD CASE NO. 1 


NODE 


EXT 

INT 

X-DISPL. 

Y-DISPL. 

Z-DISPL. 

X-ROTN. 

Y-ROIN. 

Z-ROIN. 

1 

1 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

2 

2 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

3 

3 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

4 

4 

-0.302075E+00 

0.233780E+00 

-0.331469E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

5 

5 

-0.294402E+00 

0.233737E+00 

-0.297469E+00 

-0.533925E-01 

-0.976475E-01 

-0.252952E+00 

6 

6 

-0.358100E+00 

0. 34402 9E+00 

0.558126E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

7 

7 

-0.161162E+01 

0.575101E+00 

-0. 385535E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

8 

8 

-0.125416E+01 

0.575116E+00 

-0.226667E+00 

-0.533925E-01 

-0.976475E-01 

-0.252952E+00 

9 

9 

-0.143291E+01 

0.884640E+00 

0.696849E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

10 

10 

-0.460499E+01 

0.916627E+00 

0.131711E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

11 

11 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

12 

12 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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EIEt-ENT STRESSES 


EIENENT ENDl 
NO. 

Ef©2 EM>3 

EM>4 

PX1/PX2 

SXT 

PY1/PY2 

SYT 

PZ1/PZ2 

SXYT 

W1/W2 

SXB 

MY1/MY2 

STB 

M21/MZ2 

SXYB 


ENDS END6 EN77 

END8 

SXX 

SYY 

SZZ 

SXY 

SYZ 

SZX 

1 

1 

4 


0.785575E+03 

-0.785575E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

2 

2 

4 


-0.756511E-02 

0.756511E-02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+CO 

3 

2 

5 


0.464123E+03 

-O.464123E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+CO 

4 

3 

5 


0.807432E-01 

-0.807432E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

5 

3 

6 


-0.116939E+04 

0.116939E+04 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

6 

3 

4 


-0.146366E+03 

0.146366E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

7 

4 

5 


-0.627136E-01 

0.627136E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

8 

5 

6 


0.177002E-02 

-0.177002E-02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

9 

6 

4 


0.150008E+03 

-0.150008E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
0.000 COOE'OG 

10 

4 

7 


0.705240E+03 

-0.705240E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
0. OOOOOOE+OO 

11 

5 

7 


0.452271E-01 

-0.452271E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

12 

5 

8 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

0. OOOOOOE+OO 
O.OOOOOOE+OO 

13 

6 

8 


-0.180786E+00 

0.180786E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

14 

6 

9 


-0.927916E+03 

0.927916E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

15 

6 

7 


-0.321364E+03 

0.321364E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE-OO 

16 

7 

8 


0.424805E-01 

-0.424805E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OC 

O.OOOOOOE+OO 

17 

8 

9 


0.837402E-01 

-0.837402E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

18 

9 

7 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

19 

7 

10 


0.290373E+03 

-0.290373E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

20 

8 

10 


0.290373E+03 

-0.290373E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

21 

9 

10 


-0.110159E+04 

0.110159E+04 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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4.2 Space Frame: Static Analysis 


A space frame with rigid connections, shown in figure 8 (ref. 10), is subjected to nodal forces and 
moments. Results of such analysis are presented below. 



Important data parameters: 


Young's modulus, E 
Poisson's ratio, |i 
Cross-sectional area, A 
Member length, i 


30.24 x 10 6 
0.2273 
25.13 
120 
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STARS input data; 


SPACE FRAME CASE 

6. 4. 1.4.1. 0. 0.0. 0.0 

0 , 0 , 0 , 1 , 0 , 0 , 0 , 0,0 

8 . 0 . 0 . 1 . 1 . 0 . 0.0 
2 , 0 , 1 , 0,1 

S NX1AL DATA 


1 

120.0 

0.0 

0.0 

1 

1 

1 

1 

1 

1 

2 

120.0 

0.0 

120.0 

0 

0 

0 

0 

0 

0 

3 

0.0 

0.0 

120.0 

0 

0 

0 

0 

0 

0 

4 

0.0 

-120.0 

120.0 

1 

1 

1 

1 

1 

1 

5 

0.0 

120.0 

120.0 

1 

1 

1 

1 

1 

1 

6 

10.0 

10.0 

0.0 

1 

1 

1 

1 

1 

1 

ELEMENT CONECTIVITY 








1 

1 1 

2 6 

0 0 




1 

1 


1 

2 2 

3 6 

0 0 




1 

1 


1 

3 3 

4 6 

0 0 




1 

1 


1 

4 3 

5 6 

0 0 




1 

1 



$ LINE ELEMENT BASIC PROPERTIES 

1 25.13 125.7 62.83 62.83 

$ ELEMENT MATERIAL PROPERTIES 

1 1 

30.24E06 0.2273 

S NCDAL WAD DATA 

2 4 100000.0 

2 2 4000.0 

3 1 -4000.0 

3 3 -2000.0 

3 4 -50000.0 

-1 


STARS analysis results: 


LOAD CASE NO. 1 


NODE 

EXT INT X-OISPL. Y-OISPL. Z-DISPL. 


X-ROTN. Y-ROTN. 


Z-RCTO. 


1 

2 

3 

4 

5 

6 


1 

2 

3 

4 

5 

6 


O.OOOOOOE+OO 
-0.125288E+00 
-0.125397E+00 
O.OOOOOOE+OO 
0.00000 OE+OO 
O.OOOOOOE+OO 


O.OOOOOOE+OO 

0.347953E+00 

0.103330E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 
0.196027E-04 
-0.804 94 6E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


O.OOOOOOE+OO 

-0.239969E-02 

-0.580122E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

-0.121545E-02 

-0.283265E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 
0.323397E-O2 
0 . 910380E-03 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


ELEMENT STRESSES 


ELEMENT 

NO. 

ENDl 

EN)2 EN03 

EN04 

PX1/PX2 

SXT 

PY1/PY2 

SYT 

PZ1/PZ2 

SXYT 

MX1/MX2 

SXB 

MY1/MY2 

SYB 

NE1/MZ2 

SXYB 


END5 EM) 6 EfC7 

ENDS 

sxx 

SYY 

SZZ 

SXY 

SYZ 

SZX 

1 

1 

2 


-O.124139E+03 

0.124139E+03 

-0.931688E+03 
0. 931688E+03 

0.261767E+04 

-0.261767E+04 

-0.417341E+05 

0.417341E+05 

-0.193156E+06 

-0.120964E+06 

-0.785065E+05 

-0.332961E+05 

2 

2 

3 


-0.690813E+03 

0.690813E+03 

0.232395E+03 

-0.232395E+03 

-0.129390E+04 

0.129390E+04 

0.234814E+05 
-0. 234814E+05 

0.397457E+05 

0.115522E+06 

0.255969E+05 

0.229051E+04 

3 

3 

4 


-0.654366E+03 
0.654 366E+03 

0.523179E+03 

-0.523179E+03 

-0. 980653E+03 
0.980653E+03 

0.365551E+04 

-0.365551E+04 

0.437120E+05 

0.739664E+05 

0.234344E+05 

0.393472E+05 

4 

3 

5 


0.654 3 66E+03 
-0.654366E+03 

0.131882E+04 

-0.131882E+04 

0.249337E+04 

-0.249337E+04 

-0.365551E+04 

0.365551E+04 

-0. 164730E+06 
-0.134475E+06 

0.870856E+05 

0.711728E+05 



4.3 Plane Stress: Static Analysis 


Figure 9 (ref. 11) depicts the finite element model of the symmetric half of a deep beam subjected 
to a concentrated load. Figure 10 shows the deformed shape of the structure. 


P/2 


1 2 3 4 5 16 
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Figure 9. Deep beam example. 


Important data parameters: 

Young's modulus, E = 30 x 10 6 

Poisson's ratio, jj. = 0.2 

Nodal load, P = 10,000 

Beam side length, i =20 
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STARS input data; 


PLANE STRESS CASE - TEZCAN 

66.50.1.4.0. 1.0.0. 0.0 

0 , 0 , 0 , 1 , 0 , 0 , 0 , 0,0 

8 . 0 . 0 . 1 . 1 . 0 . 0.0 
2, 0,2, 0,1 

$ NCDAL DATA 


1 


0. 


0. 


0. 

0 

0 

1 

1 

1 

0 




55 


18. 


0. 


0. 

0 

0 

1 

1 

1 

0 

0 

0 

6 

61 


20. 


0. 


0. 

1 

1 

1 

1 

1 

1 




2 


0. 


2. 


0. 

0 

0 

1 

1 

1 

0 




62 


20. 


2. 


0. 

0 

0 

1 

1 

1 

0 

0 

0 

6 

3 


0. 


4. 


0. 

0 

0 

1 

1 

1 

0 




63 


20. 


4. 


0. 

0 

0 

1 

1 

1 

0 

0 

0 

6 

4 


0. 


6. 


0. 

0 

0 

1 

1 

1 

0 




64 


20. 


6. 


0. 

0 

0 

1 

1 

1 

0 

0 

0 

6 

5 


0. 


8. 


0. 

0 

0 

1 

1 

1 

0 




65 


20. 


8. 


0. 

0 

0 

1 

1 

1 

0 

0 

0 

6 

6 


0. 


10. 


0. 

0 

1 

1 

1 

1 

1 




66 


20. 


10. 


0. 

0 

1 

1 

1 

1 

1 

0 

0 

6 

$ ELEMENT CONNECTIVITY 












2 

1 

1 

7 

8 

2 

0 

0 

0 

0 

1 

1 





2 

10 

55 

61 

62 

56 





1 

1 

0 

0 

0 

6 

2 

11 

2 

8 

9 

3 





1 

1 





2 

20 

56 

62 

63 

57 





1 

1 




6 

2 

21 

3 

9 

10 

4 





1 

1 





2 

30 

57 

63 

64 

58 





1 

1 




6 

2 

31 

4 

10 

11 

5 





1 

1 





2 

40 

58 

64 

65 

59 





1 

1 




6 

2 

41 

5 

11 

12 

6 





1 

1 





2 

50 

59 

65 

66 

60 





1 

1 




6 


$ SHELL ELUMENT THICKNESS 

1 0.1 0.1 0.1 

$ ELEMENT NATERIAL PROPERTIES 
1 1 

30.0E6 0.2 

S NCDAL LOAD DATA 
6 1 5000. 

-1 


STARS analysis results: 

LOAD CASE NO. 1 


NODE 


EXT 

INI 

X-DISPL. 

Y-DISPL. 

Z-DISPL. 

X-ROTN. 

Y-ROTN. 

Z-ROTN. 

1 

1 

0.643004E-02 

0.784168E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

2 

2 

0.686331E-02 

0.806537E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

3 

3 

0.736951E-02 

0.882648E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

4 

4 

0. 81482 6E-02 

0.991451E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

5 

5 

0.942354E-02 

0.984274E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

6 

6 

0.135581E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

7 

7 

0.645125E-02 

0.339509E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

8 

8 

0. 68583 6E-02 

0.327585E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

9 

9 

0.735772E-02 

0.281 376E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

10 

10 

0.810701E-02 

0.135669E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

11 

11 

0.944763E-02 

-0.202705E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

12 

12 

0.106295E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

13 

13 

0.650795E-02 

-0.353817E-04 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

14 

14 

0.683854E-02 

-0.546589E-04 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

15 

15 

0.7300S9E-02 

-0.111159E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

16 

16 

0.797687E-02 

-0.215918E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

17 

17 

0.877384E-02 

-0.212967E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

18 

18 

0.919536E-02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

19 

19 

0.654701E-02 

-Q.242513E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

20 

20 

0. 67588 3E-02 

-0.242862E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

21 

21 

0.714184E-02 

-0.252697E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

22 

22 

0.764636E-02 

-0.249769E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

23 

23 

0.812575E-02 

-0.171392E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

24 

24 

0.833486E-02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

25 

25 

0. 64974 8E-02 

-0.274906E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

26 

26 

0.659352E-02 

-0.248456E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

27 

27 

0.688814E-02 

-0.219499E-03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 



28 

28 

29 

29 

30 

30 

31 

31 

32 

32 

33 

33 

34 

34 

35 

35 

36 

36 

37 

37 

38 

38 

39 

39 

40 

40 

41 

41 

42 

42 

43 

43 

44 

44 

45 

45 

46 

46 

47 

47 

48 

48 

49 

49 

50 

50 

51 

51 

52 

52 

53 

53 

54 

54 

55 

55 

56 

56 

57 

57 

58 

58 

59 

59 

60 

60 

61 

61 

62 

62 

63 

63 

64 

64 

65 

65 

66 

66 


0.726972E-02 
0.759755E-02 
0.772861E-02 
0.631290E-02 
0. 63324 6E-02 
0.658259E-02 
0.691461E-02 
0.718585E-02 
0.728973E-02 
0. 59621 6E-02 
0.596853E-02 
0.625301E-02 
0.661288E-02 
0.688862E-02 
0. 69900 OE-02 
0.539932E-02 
0. 54894 7E-02 
0.592653E-02 
0.639034E-02 
0.670621E-02 
0.681496E-02 
0.451550E-02 
0.489274E-02 
0.565929E-02 
0.627 52 3E-02 
0.66251 0E-02 
0.673927E-02 
0.300094E-02 
0.427165E-02 
0.557725E-02 
0.624 71 2E-02 
0.660310E-02 
0.671494E-02 
O.OOOOOOE+OO 
0.430391E-02 
0.554128E-02 
0.623057E-02 
0. 65764 5E-02 
0.668556E-02 


-0.181145E-03 
-0.108600E-03 
O.OOOOOOE+OO 
-0.17 681 8E-03 
-0.125004E-03 
-0.84 537 6E-04 
-0 . 587820E-04 
-0.332733E-04 
O.OOOOOOE+OO 
-0. 45301 6E-06 
0.743781E-04 
0.104449E-03 
0.854564E-04 
0.441688E-04 
O.OOOOOOE+OO 
0.206086E-03 
0.303468E-03 
0.294909E-03 
0.206152E-03 
0. 983889E-04 
O.OOOOOOE+OO 
0.394097E-03 
0.491623E-03 
0.394860E-03 
0.214191E-03 
0.879983E-04 
O.OOOOOOE+OO 
0.378987E-03 
0.504907E-03 
0.139993E-03 
-0.457609E-06 
-0.241771E-04 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
-0.765803E-03 
-0.654390E-03 
-0.445990E-03 
-0.226519E-03 
O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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O.OOOOOOE+OO 

O.OOOOOOE+OO 
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O.OOOOOOE+OO 

O.OOOOOOE+OO 


O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0. OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
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O.OOOOOOE+OO 
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O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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O.OOOOOOE+OO 
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O.OOOOOOE+OO 
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ElfMENT STRESSES 


EIENENT EM)1 E1D2 EUD3 END4 
NO. 

EN15 END6 END7 END8 

PX1/PX2 

SXT 

sxx 

PY1/PY2 

SYT 

SYY 

PZ1/PZ2 

SXYT 

SZZ 

NK1/NK2 

SXB 

SXY 

MYl/MSf^ 
SYB 
SC 1 

MZ1/MZ2 

SXYB 

SZX 

1 

1 

7 

8 

2 

0.143381E+03 

0.107015E+03 

-0.260105E+03 

0.143381E+03 

0.107 01 5E+03 

-0.260105E+03 

2 

7 

13 

14 

8 

0.239308E+03 

-0.186143E+03 

-0. 607174E+02 

0.239308E+03 

-0. 18614 3E+03 

-0. 607174E+02 

3 

13 

19 

20 

14 

-0.348179E+03 

-0.216836E+03 

0.459590E+03 

-0.348179E+03 

-0.216836E+03 

0.459590E+03 

4 

19 

25 

26 

20 

-0.163769E+04 

-0.131780E+03 

O.843353E+03 

-0.163769E+04 

-0.131780E+03 

0.843353E+03 

5 

25 

31 

32 

26 

-0.335925E+04 

-O.848665E+02 

0.105358E+04 

-0. 335925E+04 

-0.848665E+02 

0.105358E+04 

6 

31 

37 

38 

32 

-0.5385S0E+04 

-0.127262E+03 

0.125524E+04 

-0.538550E+04 

-0. 127262E+03 

0.125524E+04 

7 

37 

43 

44 

38 

-0.787078E+04 

-0.282554E+03 

0.166298E+04 

-0.787078E+04 

-0.282554E+03 

0. 166298E+04 

8 

43 

49 

50 

44 

-0.112623E+05 

-0.790639E+03 

0.263614E+04 

-0.112623E+05 

-0.790639E+03 

0. 26361 4E*04 

9 

49 

55 

56 

50 

-0.163355E+05 

-0.159126E+04 

0.514415E+04 

-0.163355E+05 

-0.159126E+04 

0. 51441 5Ei 04 

10 

55 

61 

62 

56 

-0.241927E+05 

-0.963766E+04 

0.122654E+05 

-0.241927E+05 

-0.963766E+04 

0.1226548*05 

11 

2 

8 

9 

3 

-0.840131E+02 

0.207460E+03 

-0.233322E+03 

-0.840131E+02 

0.207460E+03 

-0.233322E+03 

12 

8 

14 

15 

9 

-0.761663E+03 

-0.922655E+03 

0.583229E+03 

-0.761663E+03 

-0.922655E+03 

0.583229E+03 

13 

14 

20 

21 

15 

-0.196667E+04 

-0.890847E+03 

0.161039E+04 

-0.196667E+04 

-0.89084 7E+03 

0.161039E+04 

14 

20 

26 

27 

21 

-0.324356E+04 

-0.505300E+03 

0. 220385E+04 

-0.324356E+04 

-0. 50530QE+03 

0.220385E+04 

15 

26 

32 

33 

27 

-0.431822E+04 

-0.342976E+03 

0.250988E+04 

-0.431822E+04 

-0.342976E+03 

0.25O988E+O4 


46 



16 

32 

38 

39 

33 

-0.530784E+04 

-O.532538E+03 

0.288429E+04 

-G.530784E+04 

-0.532538E+03 

0. 288429E+04 

17 

38 

44 

45 

39 

-0.625965E+04 

-0.109058E+04 

0. 356589E+04 

-0.625965E+04 

-0.109058E+04 

0.356589E-D4 

18 

44 

50 

51 

45 

-0.691433E+04 

-0.217278E+04 

0.466161E+04 

-0. 691433E+04 

-0.217278E+04 

0. 466161E+04 

19 

50 

56 

57 

51 

-0.621450E+04 

-0.470548E+04 

0.572053E+04 

-0. 621450E+04 

-0. 47054 8E+04 

0. 572053E+04 

20 

56 

62 

63 

57 

-0.425137E+03 

-0.198628E+04 

0.149334E+04 

-0.425137E+03 

-0.198628E+04 

0. 149334E+04 

21 

3 

9 

10 

4 

-0.472056E+03 

-0. 371184E+03 

0. 221830E+03 

-0.472056E+03 

-0.371184E+03 

0.221830E-*03 

22 

9 

15 

16 

10 

-0.185435E+04 

-0.224936E+04 

0.212951E+04 

-0.185435E+04 

-0.224936E+04 

0.212951E404 

23 

15 

21 

22 

16 

-0.398151E+04 

-0.156003E+04 

0.314191E+04 

-0.3981 51 E+04 

-0. 15600 3E+04 

0. 314 1 91 E + 04 

24 

21 

27 

28 

22 

-0.485999E+04 

-0. 662377E+03 

0. 308726E+04 

-0.485999E+04 

-0.662377E+03 

0. 30872 6E> 04 

25 

27 

33 

34 

28 

-0.506124E+04 

-0.531419E+03 

0.303415E+04 

-0.506124E+04 

-0.531419E+03 

0. 30341 5E+04 

26 

33 

39 

40 

34 

-0.492160E+04 

-0. 933600E+03 

0.320350E+04 

-0.492160E+04 

-0. 933600E+03 

0.320350E+04 

27 

39 

45 

46 

40 

-0.445751E+04 

-0.169963E+04 

0.354638E+04 

-0.445751E+04 

-0. 169963E+04 

0.354638E+04 

28 

45 

51 

52 

46 

-0.340814E+04 

-0.270233E+04 

0.371169E+04 

-0.340814E+04 

-0. 270233E+04 

0.371169E+04 

29 

51 

57 

58 

52 

-O.136223E+04 

-0 . 268085E+04 

0. 255089E+04 

-0.136223E+04 

-0.268085E+04 

0.255089E+04 

30 

57 

63 

64 

58 

-0.304189E+03 

0.448786E+03 

0.372647E+03 

-0.304189E+03 

0.448786E+03 

0. 372647E+03 

31 

4 

10 

11 

5 

-0.673925E+03 

-0.272642E+04 

0.179106E+04 

-0.673925E+03 

-0. 27264 2E+04 

0.179106E+04 

32 

10 

16 

17 

11 

-0.680483E+04 

-0.387664E+04 

0.554921E+04 

-0.680483E+04 

-0.387664E+04 

0.554 92lE-t04 

33 

16 

22 

23 

17 

-0.751821E+04 

-0.893684E+03 

0.401277E+04 

-0.751821E+04 

-0.893684E+03 

0. 401277E+04 

34 

22 

28 

29 

23 

-0.683326E+04 

-0.234739E+03 

0.293327E+04 

-0. 683326E+04 

-0.234739E+03 

0.293327E+04 

35 

28 

34 

35 

29 

-0.583750E+04 

-0.432099E+03 

0.248987E+04 

-0.583750E+04 

-0. 432099E+03 

0.248987E+04 

36 

34 

40 

41 

35 

-0.470407E+04 

-O.105916E+04 

0.240206E+04 

-0.470407E+04 

-0. 10591 6E+04 

0.240206E+04 

37 

40 

46 

47 

41 

-0. 33964 9E+04 

-0.179718E+04 

0.239540E+04 

-0. 33964 9E+04 

-0. 179718E+04 

0.239540E+04 

38 

46 

52 

53 

47 

-0.189864E+04 

-0.213440E+04 

0.207308E+04 

-0. 189864E+04 

-0.213440E+04 

0.207308E+04 

39 

52 

58 

59 

53 

-0.625687E+03 

-0.124948E+04 

0.118445E+04 

-0. 62S687E+03 

-0. 12494 8E+04 

0.118445E+04 

40 

58 

64 

65 

59 

-0.316120E+02 

0.146181E+04 

0.168679E+03 

-0.316120E+02 

0.146181E+04 

0. 16867 9E^03 

41 

5 

11 

12 

6 

-0.239134E+05 

-0.106444E+05 

0.129045E+05 

-0.2391 34 E+05 

-0.106444E+05 

0. 12904 SE-tOS 

42 

11 

17 

18 

12 

-0.158184E+05 

-0.461364E+02 

0.497835E+04 

-0.158184E+05 

-0.461364E+02 

0.497835E+04 

43 

17 

23 

24 

18 

-0.111853E+05 

0.645631E+03 

0.210061E+04 

-0.111853E+05 

0.645631E+03 

0.210061E+04 

44 

23 

29 

30 

24 

-Q.842533E+04 

0.414872E+03 

0.125923E+04 

-0.842533E+04 

0.414872E+03 

0.125923E+04 

45 

29 

35 

36 

30 

-0.642358E+04 

-0.220669E+03 

0. 969578E+03 

-0.642358E+04 

-0.220669E+03 

0.969578E+03 

46 

35 

41 

42 

36 

-O.468073E+04 

-0.101786E+04 

0.883437E+03 

-0.468073E+04 

-0.101786E+04 

0.883437E+03 

47 

41 

47 

48 

42 

-0.301527E+04 

-0.167224E+04 

0.826096E+03 

-0.301527E+04 

-0.167224E+04 

0.826096E+03 

48 

47 

53 

54 

48 

-0.151628E+04 

-0.170116E+04 

0. 664179E+03 

-0.151628E+04 

-0.170116E+04 

0.664179E+03 

49 

53 

59 

60 

54 

-0.461703E+03 

-0.571000E+03 

0.355763E+03 

-0.461703E+03 

-0.571000E+03 

0.355763E+03 

50 

59 

65 

66 

60 

-0.459945E+02 

0.187102E+04 

0.581728E+02 

-0.459945E+O2 

0.187102E+04 

0.581728E+02 


47 





























4.4 Plate Bending: Vibration Analysis 


A square cantilever plate was analyzed to yield the natural frequencies and associated mode 
shapes. Figure 11 depicts the plate with a 4 by 4 finite element mesh, the bottom edge along the 
axis being clamped. 


Y 



Figure 11. Square cantilever plate. 


Important data parameters: 


Young's modulus, E 
Side length, £ 

Plate thickness, t 
Poisson's ratio, ji 
Mass density, p 


10 x 10 6 
10 
0.1 
0.3 

0.259 x IQ" 3 



STARS input data: 


SQUARE 4 BY 4 PLATE NONSPIM'UNG STRUCTURE 
25, 16, 1, 4, 0, 1, 0, 0, 0, 0 

0 , 0 , 0 , 1, 0 , 0 , 0 , 0 , 0 

1 , 1 , 0, 1, 0, 0, 0, 0 

2 , 0 , 2 , 0 , 1 

1, 6, 0, 6000.0, 0.0, 0, 0.0 
$ NCDAL DATA 


1 -5.0 

0.0 


0.0 

1 

1 

1 

1 

1 

1 

0 

0 

0 

5 5.0 

0.0 


0.0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

6 -5.0 

2.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 5.0 

2.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

11 -5.0 

5.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 5.0 

5.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

16 -5.0 

7.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20 5.0 

7.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

21 -5.0 

10.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

25 5.0 10.0 

ELEMENT CONNECTIVITY 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2 11 

2 7 

6 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 4 4 

5 10 

9 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

2 5 6 

7 12 

11 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 8 9 

10 15 

14 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

2 9 11 

12 17 

16 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 12 14 

15 20 

19 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

2 13 16 

17 22 

21 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 16 19 20 25 

SHELL ELEMENT THICKNESSES 

24 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 


1 0.1 

$ ELEMLMT M\TERIAL PROPERTIES 
1 1 

1.0E+07 0.30 0.0 0.259E-3 


STARS output summary - The output summary is presented in table 5. 


Table 5. Natural frequencies of a square cantilever plate. 


Mode 

Natural frequency 
CD, rad/sec 

Nondimensional parameter 
y = g >2 Vpt/D 

1 

214.02 

3.60 

2 

506.62 

8.52 

3 

1248.40 

20.99 

4 

1538.29 

25.87 

5 

1765.53 

29.69 


Note: D = plate flexural rigidity 
= Et 3 /12(l - p 2 ) 
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4.5 General Shell: Vibration Analysis 

A cantilevered circular cylindrical shell is shown in figure 12 in which quadrilateral shell elements 
are used for structural discretization to perform a free vibration analysis. 



10141 


Figure 12. Finite element model of cylindrical shell. 

Important data parameters: 

Side length, a, b =10 

Radius, r =20 

Thickness, t =0.1 

Young's modulus, E = 29.5 x 10 6 

Poisson's ratio, p = 0.3 

Mass density, p = 0.733 x 10 -3 
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STARS input data: 


SHELL EIEMENT 

8 BY 8 

Cm/ED 

SHELL 

81, 64, 

1, 4, 0, 

1, 0, 0, 

0, 0 


0, 0, 

0, 1, 0, 

0, 0, 0, 

0 


1, 1, 

0, 1, 0, 

0, 0, 0 



2, 0, 

2, 0, 1 




1, 6, 

0, 90000.0, 0.0, ' 

0.0 


$ NHAL DATA 




1 

0.0 

0.0 

0.0 

1 

2 

1.25 

0.0 

0.0 

0 

9 

10.0 

0.0 

0.0 

0 

10 

0.0 

1.25 0.2803754 

1 

11 

1.25 

1.25 0.2803754 

0 

18 

10.0 

1.25 0.2803754 

0 

19 

0.0 

2.5 

0.478218 

1 

20 

1.25 

2.5 

0.478218 

0 

27 

10.0 

2.5 

0.478218 

0 

28 

0.0 

3.75 0.5959826 

1 

29 

1.25 

3.75 0.5959826 

0 

36 

10.0 

3.75 

0.5959826 

0 

37 

0.0 

5.0 

0.6350833 

1 

38 

1.25 

5.0 0.6350833 

0 

45 

10.0 

5.0 

0.6350833 

0 

46 

0.0 

6.25 

0.5959826 

1 

47 

1.25 

6.25 0.5959826 

0 

54 

10.0 

6.25 0.5959826 

0 

55 

0.0 

7.5 

0.478218 

1 

56 

1.25 

7.5 

0.478218 

0 

63 

10.0 

7.5 

0.478218 

0 

64 

0.0 

8.75 0.2803754 

1 

65 

1.25 

8.75 0.2803754 

0 

72 

10.0 

8.75 0.2803754 

0 

73 

0.0 

10.0 

0.0 

1 

74 

1.25 

10.0 

0.0 

0 

81 

10.0 

10.0 

0.0 

0 

$ ELEMENT CONNECTIVITY 



2 

1 1 

2 11 

10 0 

0 

2 

8 8 

9 18 

17 0 

0 

2 

9 10 

11 20 

19 0 

0 

2 

16 17 

18 27 

26 0 

0 

2 

17 19 

20 29 

28 0 

0 

2 

24 26 

27 36 

35 0 

0 

2 

25 28 

29 38 

37 0 

0 

2 

32 35 

36 45 

44 0 

0 

2 

33 37 

38 47 

46 0 

0 

2 

40 44 

45 54 

53 0 

0 

2 

41 46 

47 56 

55 0 

0 

2 

48 53 

54 63 

62 0 

0 

2 

49 55 

56 65 

64 0 

0 

2 

56 62 

63 72 

71 0 

0 

2 

57 64 

65 74 

73 0 

0 

2 

64 71 

72 81 

80 0 

0 


$ SHELL ELEMENT THICKEESSES 
1 0.1 


FREE VIBRATION ANALYSIS 


111110 0 0 
00000000 
00000001 
11111000 
00000000 
00000001 
11111000 
00000000 
00000001 
11111000 
00000000 
00000001 
11111000 
00000000 
00000001 
11111000 
00000000 
00000001 
111110 0 0 
00000000 
00000001 
11111000 
00000000 
00000001 
111110 0 0 
00000000 
00000001 

0 0 11 

00110001 
0 0 11 

00110001 
0 0 11 

00110001 
0 0 11 

0 0 1 1 0 0 0 1 
0 0 11 

00110001 
0 0 11 

00110001 
0 0 11 

00110001 
0 0 11 

00110001 


$ ELEMENT MATERIAL PROPERTIES 
1 1 

0 . 2 950E+080 . 3000E+000 . 0000E+000 . 7332E-03 
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STARS output summary - The output summary is presented in table 6. 


Table 6. Natural frequencies of a cylindrical cantilever shell. 


Mode 

Natural frequency 
co, rad/sec 

Nondimensional parameter 
y = coa^ Vpt/D 

1 

686.0745 

11.30 

2 

1108.5908 

18.26 

3 

1918.0797 

31.60 

4 

2703.7155 

44.54 

5 

2962.7536 

48.81 

6 

3904.4432 

64.32 


4.6 General Solid: Vibration Analysis 

A cube idealized by hexahedral solid elements is shown in figure 13. The nodes lying in the 
X-Y plane are assumed to be fixed. Details of the natural frequency analysis of the cube are 
presented herein. 



Figure 13. Cube discretized by hexahedral elements. 


Important data parameters: 


Side length, i 
Young’s modulus, E 
Poisson’s ratio, |i 
Mass density, p 


= 10 

= 10xl0 6 
= 0.3 

= 2.349 x 1(H 




STARS input data 


HEXAHEDRCN CASE - 2 m 2 

27 . 8 . 1 . 4 . 0 . 0 . 0 . 0 . 0.0 

0 , 0 , 0 , 1 , 0 , 0 , 0 , 0,0 

1 . 1 . 0 . 1 . 0 . 0 . 0.0 
2,0, 2,0,1 

1 , 6 , 0 , 150000 . 0 , 0 . 0 , 0.0 
$ NXAL DATA 


1 

00.0 

00.0 


00.0 

1 

1 

1 

1 

1 

1 

2 

5.0 

00.0 


00.0 

1 

1 

1 

1 

1 

1 

3 

10.0 

00.0 


00.0 

1 

1 

1 

1 

1 

1 

4 

10.0 

5.0 


00.0 

1 

1 

1 

1 

1 

1 

5 

10.0 

10.0 


00.0 

1 

1 

1 

1 

1 

1 

6 

5.0 

10.0 


00.0 

1 

1 

1 

1 

1 

1 

7 

00.0 

10.0 


00.0 

1 

1 

1 

1 

1 

1 

8 

00.0 

5.0 


00.0 

1 

1 

1 

1 

1 

1 

9 

5.0 

5.0 


00.0 

1 

1 

1 

1 

1 

1 

10 

00.0 

00.0 


5.0 




1 

1 

1 

11 

5.0 

00.0 


5.0 




1 

1 

1 

12 

10.0 

00.0 


5.0 




1 

1 

1 

13 

10.0 

5.0 


5.0 




1 

1 

1 

14 

10.0 

10.0 


5.0 




1 

1 

1 

15 

5.0 

10.0 


5.0 




1 

1 

1 

16 

00.0 

10.0 


5.0 




1 

1 

1 

17 

00.0 

5.0 


5.0 




1 

1 

1 

18 

5.0 

5.0 


5.0 




1 

1 

1 

19 

00.0 

00.0 


10.0 




1 

1 

1 

20 

5.0 

00.0 


10.0 




1 

1 

1 

21 

10.0 

00.0 


10.0 




1 

1 

1 

22 

10.0 

5.0 


10.0 




1 

1 

1 

23 

10.0 

10.0 


10.0 




1 

1 

1 

24 

5.0 

10.0 


10.0 




1 

1 

1 

25 

00.0 

10.0 


10.0 




1 

1 

1 

26 

00.0 

5.0 


10.0 




1 

1 

1 

27 

5.0 

5.0 


10.0 




1 

1 

1 

ELEMENT CCNNECTIVITy 









4 

1 1 

2 9 

8 

10 

11 

18 

17 

1 



4 

2 2 

3 4 

9 

11 

12 

13 

18 

1 



4 

3 9 

4 5 

6 

18 

13 

14 

15 

1 



4 

4 8 

9 6 

7 

17 

18 

15 

16 

1 



4 

5 10 

11 18 

17 

19 

20 

27 

26 

1 



4 

6 11 

12 13 

18 

20 

21 

22 

27 

1 



4 

7 18 

13 14 

15 

27 

22 

23 

24 

1 



4 

8 17 

18 15 

16 

26 

27 

24 

25 

1 




$ ELEMENT MATERIAL PROPERTIES 
1 1 

1 . 0 E +7 0.3 0.0 2 . 349 E -4 


STARS output summary - The output summary is presented in table 7. 


Table 7. Natural frequencies of a solid cube. 


Mode 

Natural frequency parameter 
a> = oy V ( E/p) , rad/sec 

Exact solution 

✓s j 

CO 

Mesh size 

2 by 2 

4 by 4 

6 by 6 

1 

0.07975 

0.07195 

0.06958 j 

0.06801 

2 

0.07975 

0.07195 

0.06958 

0.06801 

3 

0.1315 

0.1043 

0.09762 

0.09288 

4 

0.1720 

0.1645 

0.1623 

0.1611 

5 

0.2280 

0.1933 

0.1850 

0.1819 

6 

0.2280 

0.1933 

0.1850 

0.1819 
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4.7 Spinning Cantilever Beam: Vibration Analysis 


A cantilever beam spinning about the Y-axis is shown in figure 14. 



Figure 14. Spinning cantilever beam. 


Important data parameters - The structure is assumed to possess both viscous and structural 
damping. 


Young's modulus, E 

= 30 x 10 6 

Cross-sectional area, A 

- 1.0 

Moment of inertia: 


About Y-axis 

= 1/12 

About Z-axis 

= 1/24 

Element length, £ 

= 6 

Nodal translational mass 

= 1 

Nodal mass moment of inertia 

= 1/35 

Scalar viscous damping 

= 0.628318 

Structural damping coefficient 

= 0.01 

Spin rate, Hz 

= 0.1 
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STARS input data 


SPINNING CANTIIEVER BEAM - 10-EUMENT IDEALIZATION - VISC AND STRUCT DAMPING 

12.10.1.4.1.0. 0.0.0.0 

0 , 0 , 1 , 1 , 0 , 0 , 0 , 0,0 

5. 0. 0. 0.0. 1.0.0 

2 . 0 . 2 . 0.1 

1,6,0,500.0,0.0,0.0 

0.01 

$ NXlAL DATA 


1 

0.0 

0.0 


0.0 

1 

1 

1 

1 

1 

1 

2 

6.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

3 

12.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

4 

18.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

5 

24.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

6 

30.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

7 

36.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

8 

42.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

9 

48.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

10 

54.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

11 

60.0 

0.0 


0.0 

0 

0 

0 

0 

0 

0 

12 

25.0 

15.0 


0.0 

1 

1 

1 

1 

1 

1 

ELEMENT CONNECTIVITY 









1 

1 1 

2 12 

0 

0 

0 

0 

0 

1 

1 

0 

1 

10 10 

11 12 

0 

0 

0 

0 

0 

1 

1 

0 


$ LINE ELEMENT BASIC PROPERTIES 

1 1.0 0.125000.083333330.04166667 

$ ELEMENT NATERIAL PROPERTIES 
1 1 

30.0E+06 0.30 

S ELEMENT SPIN RATE DATA 

1 0.0 0.628318 0.0 

$ NCDAL NASS DATA 

21 1.0 3 

31 1.0 3 

41 1.0 3 

51 1.0 3 

61 1.0 3 

71 1.0 3 

81 1.0 3 

91 1.0 3 

10 1 1.0 3 

11 1 1.0 3 

2 4 0.0285714 6 

3 4 0.0285714 6 

4 4 0.0285714 6 

5 4 0.0285714 6 

6 4 0.0285714 6 

7 4 0.0285714 6 

8 4 0.0285714 6 

9 4 0.0285714 6 

10 4 0.0285714 6 

11 4 0.0285714 6 

-1 

$ VISCOUS DAMPING DATA 

0 . OOOOE+OOO . OOOOE+OOO . OOOOE+OOO . OOOOE+OOO . OOOOE+OOO . 0000E+00 
0.628318000.628318000.628318000. 628318000 . 628318000 . 62831800 
0.628318000.628318000. 628318000. 628318000 . 628318000. 62831800 
0.628318000.628318000.628318000.628318000.628318000.62831800 
0 . 62831 8000 .628318000. 628318000. 62 8318000 .62831 8000 .62831800 
0.628318000.628318000.628318000.628318000.628318000.62831800 
0 . 628318000 . 628318000. 628318000. 628318000.628318000.62831800 
0 . 628318000 . 628318000. 628318000. 628318000 . 62831 8000 . 62831800 
0 . 628318000. 628318000. 628318000. 628318000 . 62831 8000. 62831800 
0 . 628318000.628318000. 628318000. 628318000 . 628318000.62831800 
0. 628318000.628318000. 628318000. 628318000.628318000.62831800 
0 . OOOOE+OOO . OOOOE+OOO. OOOOE+OOO . OOOOE+OOO . OOOOE+OOO . 0000E+00 
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STARS output summary - The output summary is presented in table 8. 


Table 8. Spinning cantilever beam. 


Mode 

Structure without 
damping, 
IPROB = 2 

Structure with 
viscous damping, 
IPROB = 4 

Structure with viscous 
and structural damping, 
IPROB = 5 

1 

2.5263 

-.3107 ± 2.4886i* 

-.3195 ± 2.4820i* 

2 

3.4481 

-.3116 ± 3.4200i* 

-.3255 + 3.41231* 

3 

15.3964 

-.3169 + 15.3865i* 

-.3930 ± 15.383 li* 

4 

21.7055 

-.3166 ± 21.7002i* 

-.4243 ± 2 1 .69 1 2i* 

5 

43.1614 

-.3202 ± 43. 1 398i* 

-.4848 ± 43.06271* 

6 

60.9511 

-.3202 ± 60.949 li* 

-.6246 ± 60.9390i* 


Note: Natural frequencies for various problem types are due to a spin rate Q = 0.1 Hz (0.6283 
rad/sec) (i* = iTT). 


4.8 Spinning Cantilever Plate: Vibration Analysis 

The cantilever plate model described in section 4.4 is chosen for this sample problem. The plate 
is spun along the Z-axis with a uniform spin rate £2% = 0.8*co n , co n being the first natural frequency 

of vibration of the nonrotating plate. Table 9 provides the first few natural frequencies of the plate 
in nondimensional form, co being the natural frequencies. Also presented in the table are the results 
of the free vibration analysis of the plate rotating along an arbitrary axis, the spin rate being 

Or ~ 0.8 x coi, with components Q x ~ &Y = * 0.8 co^/ 
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STARS input data: 


9XRRE 4 BY 4 PIATE SPIM4ING STRUCTURE 

25.16.1.4.0. 1.0.0.0.0 

0 , 0 , 1 , 1 , 0 , 0 , 0 , 0,0 

2 . 0 . 0 . 1 . 0 . 0 . 0.0 
2, 0,2, 0,1 

1,6,0,6000.0,0.0,0.0 
$ NCDAL DATA 


1 


-5.0 


0.0 


0.0 

1 

1 

1 

1 

1 

1 

0 

0 

0 

5 


5.0 


0.0 


0.0 

1 

1 

1 

1 

1 

1 

0 

0 

1 

6 


-5.0 


2.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 


5.0 


2.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

11 


-5.0 


5.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 


5.0 


5.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

16 


-5.0 


7.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20 


5.0 


7.5 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

21 


-5.0 


10.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

25 


5.0 


10.0 


0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

ELEMENT CCNNEX7ITVITY 












2 

1 

1 

2 

7 

6 

0 

0 

0 

0 

1 

1 

0 

0 

1 


2 

4 

4 

5 

10 

9 

0 

0 

0 

0 

1 

1 

0 

0 

1 

1 

2 

5 

6 

7 

12 

11 

0 

0 

0 

0 

1 

1 

0 

0 

1 


2 

8 

9 

10 

15 

14 

0 

0 

0 

0 

1 

1 

0 

0 

1 

1 

2 

9 

11 

12 

17 

16 

0 

0 

0 

0 

1 

1 

0 

0 

1 


2 

12 

14 

15 

20 

19 

0 

0 

0 

0 

1 

1 

0 

0 

1 

1 

2 

13 

16 

17 

22 

21 

0 

0 

0 

0 

1 

1 

0 

0 

1 


2 

16 

19 

20 

25 

24 

0 

0 

0 

0 

1 

1 

0 

0 

1 

1 


$ SHELL ELEMENT THICKNESSES 
1 0.1 


$ ELEMENT MATERIAL PROPERTIES 
1 1 

1.0E+07 0.30 0.0 0.259E-3 

$ ELEMENT SPIN RATE DATA 

1 0. 0. 170.86 


STARS output summary - The output is presented in table 9. 


Table 9. Natural frequency parameters of a 
spinning square cantilever plate. 



Natural frequency parameter „ 2 — — 

7=0^ V pt/D 

Mode 

Qz = 0.8co ^ 

£2r = 170.86 rad/sec, 


= 100.00 rad/sec 

Qx = Qy = Qz 

= 57.735 rad/sec 


CO 

y 

(0 

y 

1 

233.57 

3.9286 

146.36 

2.4617 

2 

534.17 

8.9845 

494.34 

8.3146 

3 

1268.30 

21.3323 

1248.60 

21.0009 

4 

1562.00 

26.2722 

1543.26 

25.9570 

5 

1784.80 

30.0196 

1768.80 

29.7505 

6 

2912.50 

48.9871 

2898.28 

48.7479 
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4.9 Helicopter Structure: Vibration Analysis 

A coupled helicopter rotor-fuselage system is shown in figure 15 (ref. 12) along with relevant 
stiffness and mass distributions, which are suitably approximated for the discrete element modeling 
of the structure. Numerical free vibration analysis was performed for the structure with the rotor 
spinning at 10 rad/sec (fly = 10); such results are presented in table 10, along with the results for 
the corresponding nonspinning case. 


1 2 3 4 5 6 7 8 

I 1 1 1 1 — I — I H 


Hub 


-1 ! — H 



Q Y 

1011 12 1314 151617 

i — i — i — i — i — i — t — i 

V 


H 1 1 1_ 


Rotor 
y— Fuselage 


18 19 20 21 22 23 24 25 26 27 28 2930 


(a) Discrete element model. 


16 

12 


Mass, 

slug/ft 


Point mass (node 9) = 16.0 slug 

/-Fuselage 




Rotor 


0 -20 -10 0 10 20 
Length, ft 


10145 


(b) Structural mass distribution. 


16 x 10 6 



Length, ft 

10146 

(c) Structural stiffness distribution. 

Figure 15. Coupled helicopter rotor-fuselage system. 
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STARS input data: 

HELICOPTER STFUCTORE, JOURNAL CASE, SPIN = 10.0 (FIXED UZ, UXR, UYR) 

31.29.7.4.2.0. 0.0.0.0 

0 , 0 , 1 , 0 , 0 , 0 , 0 , 0,0 

2 . 0 . 0 . 1 . 0 . 1 . 0.0 
2, 0,2, 0,1 

1,12,0,53.15,0.0,0.0 
$ NXftL DATA 

1 -25.00000 1.0 0.0 0 0 1 1 1 0 0 

8 -3.12500 1.0 0.0 0 0 1 1 1 0 0 0 1 

9 0.00000 1.0 0.0 0 0 1 1 1 0 0 

10 3.12500 1.0 0.0 0 0 1 1 1 0 0 

17 25.00000 1.0 0.0 0 0 1 1 1 0 0 0 1 

18 -20.00000 0.0 0.0 0 0 1 1 1 0 0 

23 -3.33333 0.0 0.0 0 0 1 1 1 0 0 0 1 

24 0.00000 0.0 0.0 0 0 1 1 1 0 0 

25 3.33333 0.0 0.0 0 0 1 1 1 0 0 

30 20.00000 0.0 0.0 0 0 1 1 1 0 0 0 1 

31 10.00000 0.5 0.0 1 1 1 1 1 1 0 

$ ELEMENT CONNECTIVITY 


1 

1 

1 

2 

24 

1 

1 

1 

16 

16 

17 

24 

1 

1 

1 

17 

18 

19 

9 

2 

1 

1 

18 

19 

20 

9 

2 

1 

1 

19 

20 

21 

9 

3 

1 

1 

20 

21 

22 

9 

4 

1 

1 

21 

22 

23 

9 

5 

1 

1 

22 

23 

24 

9 

6 

1 

1 

23 

24 

25 

9 

6 

1 

1 

24 

25 

26 

9 

5 

1 

1 

25 

26 

27 

9 

4 

1 

1 

26 

27 

28 

9 

3 

1 

1 

27 

28 

29 

9 

2 

1 

1 

28 

29 

30 

9 

2 

1 

1 

29 

9 

24 

31 

7 

2 


$ LINE ELEMENT BASIC PROPERTIES 

1 1 . 1 . 1 . 1 . 

2 100 . 1 . 1 . 1 . 

$ ELEMENT MATERIAL PROPERTIES 
1 1 


2.0E05 

0.3 

0. 

0.3 

2 1 
1.53E06 

0.3 

0. 

1.23 

3 1 

2.E06 

0.3 

0. 

1.23 

4 1 

4.8E06 

0.3 

0. 

8.6 

5 1 

7.8E06 

0.3 

0. 

16. 

6 1 
11.E06 

0.3 

0. 

16. 

7 1 

1.E08 

0.3 

0. 

0. 


S ELEMENT SPIN RATE DATA 

1 0.0 10.0 0.0 

$ NCXAL NASS DATA 

9 1 16. 3 

-1 
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STARS output summary - The output summary is presented in table 10. 


Table 10. Natural frequencies of a helicopter structure. 


Mode number 

Natural frequencies, spin rates 

Mode shape 

Dy= 0 

n 

♦—* 

o 

1,2,3 

0 

0 

Rigid body 

4 

4.642 

11.789 

Rotor 1st symmetric bending 

5 

5.041 

11.793 

Rotor 1st antisymmetric bending 

6 

22.138 

22.229 

Fuselage 1st bending 

7 

27.892 

36.199 

Rotor 2nd antisymmetric bending 

8 

28.278 

37.800 

Rotor 2nd symmetric bending 

9 

37.176 

38.478 

Rotor 3rd antisymmetric bending 


4.10 Rocket Structure: Dynamic Response Analysis 

A rocket is simply idealized by four line elements, as shown in figure 16 (ref. 4), which is 
subjected to a pulse loading function at the base. Results of the dynamic response analysis are 
shown in figures 17 and 18. 


t x 


4 

3 

2 

1 

P 1 = P# 


10147 

(a) Rocket structure. 


p^t) 


10148 

(b) Pulse loading. 


Figure 16. Rocket subjected to dynamic loading. 

Important data parameters * Arbitrary element and material properties data are assumed for the 
analysis to correlate results with available ones expressed in parametric form. 


Young's modulus, E 

- 100 

Poisson's ratio, p 

= 0.3 

Cross-sectional area, A 

= 1.0 

Mass density, p 

= 1.0 
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Length of an element, i 

= 2.5 

Pulse load intensity, Po 

= 10.0 

Duration of load, sec 

= 1.0 

Total time period for 


response evaluation 

= 2.0 


STARS input data: 


DYNAMIC RESPONSE CASE - PRZEMIENIECKI 

6. 4. 1.4. 1.0. 0.0. 0.0 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0,0 

1 . 0 . 1 . 1 . 0 . 0 . 0.0 
2, 0,2, 0,1 

1,3,0,20.0,0.0,0.0 
0, 1,1,2 
$ NXlAL DATA 


i 

0.0 

0.0 

0.0 

0 

1 

1 

1 

1 

1 




5 

10.0 

0.0 

0.0 

0 

1 

1 

1 

1 

1 

0 

0 

1 

6 

5.0 

5.0 

0.0 

1 

1 

1 

1 

1 

1 





$ ELEMENT CONNECTIVITY 

1112600 11 

1445600 110001 

$ LINE ELEMENT BASIC PROPERTIES 

1 1.0 0.0 0.0 0.0 
$ ELEMENT MATERIAL PROPERTIES 
1 1 

100.0 0.3 0.0 1.0 

$ DYNAMIC NDDAL FORCE DATA 

1.0 

1 1 10.0 

-1 

$ INCREMENIAL TIME DATA FOR DYNAMIC RESPONSE ANALYSIS 
0.10 10 
0.20 5 


STARS analysis results at a typical time step: 


DYNAMIC RESPONSE AT TIME - 0.7000E+00 


NODE 

EXT INT X-DISPL. Y-DISPL. Z-DISPL. X-ROTN. YHRDTN. Z-ROIN. 

1 1 0.646322E+O0 O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

2 2 0.490999E+00 O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

3 3 0.191562E+00 O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

4 4 -0.102987E-02 O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

5 5 -0. 49508 0E-01 O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

6 6 O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

EIEMENT STRESSES 


EUMENT END1 EM)2 END3 END4 
NO. 

END5 END 6 END? END8 

PX1/PX2 

SXT 

sxx 

PY1/PY2 

SYT 

SYY 

PZ1/PZ2 

SXYT 

szz 

MX1/MK2 

SXB 

SXY 

MY1/MY2 

SYB 

SYZ 

MZ1/MZ2 

SXYB 

SZX 

1 

1 

2 

-0.1964 00E+02 
0.196400E+02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

2 

2 

3 

0.119775E+02 

-0.119775E+02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

3 

3 

4 

0.770366E+01 

-0.770366E+01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

4 

4 

5 

0.193913E+01 

-0.193913E+01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 
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4.11 Plate, Beam, and Truss Structures: Buckling Analysis 


A buckling analysis was performed for a simply supported square plate model, described in sec- 
tion 4.4, subjected to a uniform unit stress acting along the two edges parallel to the y-axis; relevant 
input data and analysis results are as follows. 


STARS input data: 


SQUARE 4 W 4 PLATE, BC W=0, midline X=0,Y(13}=0, BUCKLING ANALYSIS 

25.16.1.4.0. 1.0.0.0.0 

0 , 0 , 0 , 1 , 0 , 0 , 0 , 0,0 

9.1. 0. 1. 1. 0.0.0 

2 . 0 . 2 . 0.1 

1 , 1 , 0 , 20000 . 0 , 0 . 0 , 0.0 

S MTlAL DATA 


1 

-5.00 

0.0 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

2 

-2.50 

0.0 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

3 

0.0 

0.0 


0.0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

4 

2.50 

0.0 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

5 

5.00 

0.0 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

6 

-5.00 

2.50 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

7 

-2.50 

2.50 


0.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

8 

0.0 

2.5 


0.0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

9 

2.50 

2.50 


0.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

10 

5.00 

2.50 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

11 

-5.00 

5.00 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

12 

-2.50 

5.00 


0.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

13 

0.0 

5.00 


0.0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

14 

2.50 

5.00 


0.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

15 

5.00 

5.00 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

16 

-5.00 

7.50 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

17 

-2.50 

7.50 


0.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

18 

0.0 

7.5 


0.0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

19 

2.50 

7.50 


0.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

20 

5.00 

7.50 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

21 

-5.00 

10.00 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

22 

-2.50 

10.00 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

23 

0.0 

10.0 


0.0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

24 

2.50 

10.00 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

25 

5.00 

10.00 


0.0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

ELEMENT CONNECTIVITY 












2 

1 1 

2 7 

6 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 

4 4 

5 10 

9 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

2 

5 6 

7 12 

11 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 

8 9 

10 15 

14 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

2 

9 11 

12 17 

16 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 

12 14 

15 20 

19 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

2 

13 16 

17 22 

21 

0 

0 

0 

0 

1 

1 

0 

0 

0 


2 

16 19 

20 25 

24 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 


$ SHELL ELEMENT THICKNESSES 

1 0.1 

$ ELEMENT MATERIAL PROPERTIES 
1 1 

1.0 E +07 0.30 

$ NCDAL LOAD DATA 
1 1 .125 

6 1 .250 

11 1 .250 

16 1 .250 

21 1 .125 

5 1 -.125 

10 1 -.250 

15 1 -.250 

20 1 -.250 

25 1 -.125 

-1 
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STARS analytical results - The analytical results pertaining to the buckling load are presented in 
table 11. 


Table 1 1. Critical load of a simply supported square plate. 


Buckling load parameter for Mode 1 

STARS solution 

Exact 

solution 

4 bv 4 8 bv 8 14 bv 14 

3552.892 3562.814 3575.978 

3615.240 j 


The cantilever beam defined in section 4.7 is the subject of the buckling analysis; the relevant 
details are given below. 


STARS input data: 


CANTILEVER BEAM - 10-EIENENT IDEALIZATION - BUCKLING ANALYSIS 
C 

C TEMPERATURE LOADING ADDED 
C 

12.10.1.4.1.0. 0.0.0.0 

1 . 0 . 0 . 1 . 0 . 0 . 0 . 0.0 

9. 1.0. 1.1. 0.0.0 

2 . 0 . 2 . 0.1 

1 , 1 , 0 , 12000 . 0 , 0 . 0 , 0.0 
$ MCDAL DATA 


i 

0.0 

0.0 


0.0 

1 

1 

1 

1 

1 

1 

2 

6.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

3 

12.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

4 

18.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

5 

24.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

6 

30.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

7 

36.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

8 

42.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

9 

48.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

10 

54.0 

0.0 


0.0 

0 

0 

1 

1 

1 

0 

11 

60.0 

0.0 


0.0 

0 

1 

1 

1 

1 

0 

12 

25.0 

15.0 


0.0 

1 

1 

1 

1 

1 

1 

ELEMEJff CCNNECnVITY 









1 

1 1 

2 12 

0 

0 

0 

0 

0 

1 

1 

1 

1 

10 10 

11 12 

0 

0 

0 

0 

0 

1 

1 

1 


S LINE ELEMENT BASIC PROPERTIES 

1 1.0 0.125 0.083333 0.041667 

$ ELEMENT MATERIAL PROPERTIES 

1 1 

30.0E+06 0.30 6.6E-06 

$ ELEMENT TEMPERATURE DATA 

1 - 1.0 

$ NCDAL IQAD DATA 
11 1 - 1.0 
-1 


STARS analytical results - The analytical results are presented in table 12. 
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Table 12. Critical load of a cantilever beam. 


Mode 

Buckling load parameter 

STARS solution 

Exact 

solution 

1 

35.2320 

35.2282 


The simple truss of figure 19 (ref. 4) is also analyzed to determine the critical loads. The 
associated input data and analytical results are given below. 


p 



10151 

Figure 19. Truss structure. 


STARS input data: 


PRZ - ITUSS BLCKLINS Af^LYSIS 

4. 2. 1.4. 1.0. 0.0. 0.0 

0 , 0 , 0 , 1 , 0 , 0 , 0 , 0,0 

9.0. 0. 1.1. 0.0.0 

2 . 0 . 2 . 0.1 

1 , 2 , 0 , 20000 . 0 , 0 . 0 , 0.0 

$ MEAL DATA 


1 100.0 

100.0 

0.0 

0 

0 

1 

1 

1 

1 

2 100.0 

0.0 

0.0 

1 

1 

1 

1 

1 

1 

3 0.0 

0.0 

0.0 

1 

1 

1 

1 

1 

1 

4 0.0 

50.0 

0.0 

1 

1 

1 

1 

1 

1 

ELMOT CCNNBCTTVm 








113 

1 4 

1 1 

0 

0 

0 

1 

1 


12 2 

1 4 

1 1 

0 

0 

0 

1 

1 



$ LINE ELEMENT BASIC PROPERTIES 
1 0.1 


$ ELEMttT MATERIAL PROPERTIES 
1 1 

10.0E03 0.2 

S MEAL LOAD DATA 
1 2 - 1.0 
-1 


STARS analytical results - The analytical results are presented in table 13. 
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Table 13. Critical load of a simple truss. 


| Mode 

Buckling load parameter 

STARS solution 

Exact 

solution 

1 

261.20388 

261.20387 


4.12 Composite Plate Bending: Vibration Analysis 

To illustrate the multiple material angle (as in layered elements), the diverse coordinate system 
capabilities, and thermal effects, a composite plate problem (fig. 20) similar to that in section 4.4 is 
considered for vibration analysis. 



Figure 20. Square composite plate. 


Important data parameters: 

Side length, t =10 

Plate thickness, t = 0.063 

Mass density, p = 0.259 x 10*3 

Material properties - anisotropic, as shown in input data. 
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STARS input data 


TRIANGULAR 4 BY 4 PIATE / COMPOSITE IAYERS / IUGCS TEST CASE / 1MANGL TEST CASE 
25, 32, 1, 13, 0, 0, 2, 8, 10, 2 

1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 

1 , 1 , 0 , 1 , 0 , 0 , 0 , 0 

2 , 0 , 2 , 0 , 1 

1, 10, 0, 6000.0, 0.0, 0.0 
$ NCDAL DATA 


1 

-5.0 

0.0 

0.0 

1 

1 

1 

1 

1 

1 

1 

0 

0 

5 

5.0 

0.0 

0.0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

6 

-2.5 

-10.0 

0.0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

10 

-2.5 

0.0 

0.0 

0 

0 

0 

0 

0 

0 

2 

0 

1 

11 

-5.0 

5.0 

0.0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

15 

5.0 

5.0 

0.0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

16 

-2.5 

0.0 

0.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20 

-2.5 

-10.0 

0.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

21 

-5.0 

10.0 

0.0 

1 

1 

1 

1 

1 

1 

1 

0 

0 

25 

5.0 

10.0 

0.0 

1 

1 

1 

1 

1 

1 

1 

0 

1 


$ IOCAL-GLOBAL COORDINATE SYSTEM DATA 
1 2 



-10.0 


-5.0 


0.0 


0.0 


-1.0 


0.0 






1.0 


0.0 


0.0 


0.0 


0.0 


1.0 





2 

1 
















-10.0 

- 

-10.0 


0.0 

- 

-15.0 

- 

-10.0 


0.0 





- 

-15.0 

- 

-15.0 


0.0 











ELEMENT CONNECTIVITY 













7 

1 

1 

2 

7 

0 

0 

0 

0 

1 

1 

0 

1 

0 

0 


7 

4 

4 

5 

10 

0 

0 

0 

0 

1 

1 

0 

1 

0 

0 

1 

7 

5 

6 

7 

12 

0 

0 

0 

0 

2 

1 

0 

1 

0 

0 


7 

8 

9 

10 

15 

0 

0 

0 

0 

2 

1 

0 

1 

0 

0 

1 

7 

9 

11 

12 

17 

0 

0 

0 

0 

3 

1 

0 

1 

0 

0 


7 

10 

12 

13 

18 

0 

0 

0 

0 

3 

1 

0 

1 

0 

0 


7 

11 

13 

14 

19 

0 

0 

0 

0 

4 

1 

0 

1 

0 

0 


7 

12 

14 

15 

20 

0 

0 

0 

0 

4 

1 

0 

1 

0 

0 


7 

13 

16 

17 

22 

0 

0 

0 

0 

5 

1 

0 

1 

0 

0 


7 

16 

19 

20 

25 

0 

0 

0 

0 

5 

1 

0 

1 

0 

0 

1 

7 

17 

7 

6 

1 

0 

0 

0 

0 

6 

1 

0 

1 

0 

0 


7 

20 

10 

9 

4 

0 

0 

0 

0 

6 

1 

0 

1 

0 

0 

1 

7 

21 

12 

11 

6 

0 

0 

0 

0 

7 

1 

0 

1 

0 

0 


7 

24 

15 

14 

9 

0 

0 

0 

0 

7 

1 

0 

1 

0 

0 

1 

7 

25 

17 

16 

11 

0 

0 

0 

0 

8 

1 

0 

1 

0 

0 


7 

26 

18 

17 

12 

0 

0 

0 

0 

8 

1 

0 

1 

0 

0 


7 

27 

19 

18 

13 

0 

0 

0 

0 

9 

1 

0 

1 

0 

0 


7 

28 

20 

19 

14 

0 

0 

0 

0 

9 

1 

0 

1 

0 

0 


7 

29 

22 

21 

16 

0 

0 

0 

0 

10 

1 

0 

1 

0 

0 


7 

32 

25 

24 

19 

0 

0 

0 

0 

10 

1 

0 

1 

0 

0 

1 


$ COMPOSITE SHELL ELEMENT STACK DESCRIPTION DATA 
1 2 

1 .0315 3 

1 .0315 1 

2 2 

1 .0315 3 

1 .0315 3 

3 2 

1 .0315 4 

1 .0315 4 

4 2 

1 .0315 2 

1 .0315 2 

5 2 

1 .0315 1 

1 .0315 1 

6 2 

1 .0315 6 

1 .0315 5 

7 2 

1 .0315 6 

1 .0315 6 

8 2 

1 .0315 7 

1 .0315 7 

9 2 

1 .0315 8 

1 .0315 8 

10 2 

1 .0315 5 

1 .0315 5 

$ SPECIFICATION FCR MATERIAL AXIS ORIENTAIICN 
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0 


1 2 

0.0 

2 2 0 

0.7854 

3 2 0 

1.5708 

4 2 0 

2.3562 

5 2 0 

3.1416 

6 2 9 

4.7120 

7 2 0 

5.4980 

8 2 0 

3.9260 

$ ELB-ENT JMERIAL PROPERTIES 
1 2 

.7966E+07 . 6638E+06 0.0 .2566E+07 0.0 0.125E+7 .1042E+7 

0.0 .1042E+07 .35e-5 .114e-4 0.0 0.259E-3 

$ ELEMENT TEMPERATURE DATA 

1 10 . 


STARS output summary - The results are printed in table 14. 


Table 14. Natural frequencies of a 
composite square cantilever plate. 


Mode 

Natural frequency 
o), rad/sec 

T = 0 

o 

t— i 

n 

1 

505.41 

373.38 

2 

611.98 

486.79 

3 

967.78 

851.04 

4 

1434.66 

1275.96 

5 

1523.71 

1361.70 

6 

1765.22 

1643.97 




5. STARS-AERO AND ASE PROGRAM DESCRIPTION 


The aeroelastic and aeroservoelastic modules (fig. 2) are recent additions to the original STARS pro- 
gram (ref. 1) that are capable of predicting related stability of such structures as aircraft and spacecraft. 
Thus, once the vibration analysis is performed utilizing the STARS-SOLEDS module, the program con- 
tinues to determine flutter and divergence characteristics as well as open- and closed-loop stability anal- 
yses, as desired. In this connection, a typical feedback control system is shown in figure 21, whereas 
figure 22 presents a flowchart depicting the ASE analysis methodology. References 13 and 14 provide 
some details of the current analysis techniques. 

Detailed numerical formulation in connection with the present aero-structural-control analysis is 
given in section 5.1. The unsteady aerodynamic forces for supersonic flow are computed by a constant 
pressure method (CPM) (ref. 15), whereas the doublet lattice method (DLM) (refs. 16,17) is utilized for 
the subsonic case. Both k and p-k stability (flutter and divergence) solution procedures are available to 
the user. 

For the ASE analysis, the aerostructural problem is recast in the Laplace domain when the general- 
ized aerodynamic forces are curve-fitted using Pade and least squares approximations, thereby yielding 
the state-space matrices. Such matrices can then be augmented by analog elements such as actuators, 
sensors, prefilters, and notch filters, and also the analog controller. Associated equivalent open-loop 
(loop-gain) or open-loop transfer function is obtained by standard procedure, whereas the closed-loop 
formulation is derived similarly by appropriately taking into account the feedback equation. The system 
frequency responses are simply obtained from the appropriate transfer matrices. Associated modal 
damping and frequency values may also be derived by solving the eigenvalue problem of the augmented 
state-space plant dynamics matrix. 

In the case of a digital controller, a hybrid equivalent open-loop or closed-loop transfer function is 
achieved by suitably combining the controller, the open-loop transfer function of the original analog 
system of the plant, and other analog elements; frequency responses are then obtained in a routine man- 
ner. The modal damping and frequency values are obtained by first transferring the augmented analog 
state-space plant dynamics matrix from its usual Laplace (s) to the digital z-plane, adding the same to the 
corresponding matrix for the controller, and finally solving the associated eigenvalue problem. 

Furthermore, the open-loop stability analyses (flutter and divergence) may also be effected with or 
without the controller (analog or digital). This is achieved by solving eigenvalue problems of the appro- 
priately augmented and transformed, as the case may be, plant dynamics matrix for a number of reduced 
frequency values and noting the change in sign of the real part of the eigenvalues. Such a solution with- 
out a controller can be compared with the aeroelastic analysis using the k and p-k methods, whereas the 
relevant solution in the presence of a controller proves to be useful for comparing relevant flight test re- 
sults of modem, high-performance, unstable aircraft. 
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Summing junction outputs 
u = r - Gy 
X = AX + Bu 
y = CX + Du 
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Figure 21. Feedback control system. 


5. 1 Numerical Formulation for Aeroelastic and Aeroservoelastic Analysis 

In the numerical formulation presented here, structural discretization is based on the finite 
element method, whereas the panel methods are adopted for computation of unsteady aerodynamic 
forces. The more specialized matrix equation of motion of such structures relevant to the current 
analysis has the form 


Mq + Cq + Kq + #A e (k)q = P(t) ( 24 ) 

in which relevant terms are defined as follows: 

M inertia matrix, 

C damping matrix, 

K elastic stiffness matrix, 

q dynamic pressure l/2pV 2 , p and V being the air density and true airspeed, 
respectively, 

k reduced frequency cob/V, co and b being the natural frequency and wing semichord 
length, respectively, 

A e (k) aerodynamic influence coefficient matrix for a given Mach number Moo and set 
of ki values, 

q displacement vector, 

P(t) external forcing function, and 

s Laplace variable (= i*co, i* being V— 1). 

A solution (ref. 1) of the related free vibration problem 


Mq + kq = 0 


( 25 ) 


yields the desired roots co and vectors C>. Next, applying a transformation 

q = On (26) 
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to equation (24) and premultiplying both sides by d> T , the generalized equation of motion is derived as 


Mrj + Cq + ICq + <7 Q(k)q = P(t) 


(27) 


in which M = O t MO, etc., the modal matrix <J> = [O x O e $ 5 ], and the generalized coordinate 
q = [q r q e q§] incorporates rigid body, elastic, and control surface motions, respectively. 

Expressing the generalized aerodynamic force matrix Q(k) as Pade polynomials (ref. 4) in 
i*k (=i*cob/V = sb/V) results in 


i*k i*k 

Q(k) = A 0 + i*kAi + (i*k) 2 A 2 + i * k+ p i A 3 + A 4 
where Pj are the aerodynamic lag terms (assuming j = 1 , 2 ), and 


+ . 


i*k 


k 2 


i*k+pj k 2 +pj 2 + k 2 +pj 2 
Further, separating the real and imaginary parts in equation (28) yields 

QrOO = (QrOO - Ao) 

-k 2 I k2 -I — — — i 
k 2 + (J 2 k 2 + p 2 

= S R (k) A 

Ql(k)= QiW/k-Aj 


i*kPj 


1 

A2 


a 3 

- 

la 4 J 


0 




-I 




k z + pf k z + P^ J 


*1 

A 2 


a 3 


-A 4 _ 


= S^k) A 

in which for a small value of k = kj, the coefficients assume the following form: 

Ao = Qrfki) 

Qi(ki) A 3 A 4 


Ai = 


ki Pi P2 


(28) 


(28a) 


(29) 


(29a) 


(30) 

(30a) 


Substituting equation (30a) in equation (29a), the unknown coefficients A 3 and A 4 can be determined. 
However, the resulting solution is sensitive to the choice of pj. On the other hand, if the elements of the 
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Ai matrix are replaced by measured damping coefficients without any lag terms, then the solution is 
insensitive to the pj values. 

Equations ( 29 ) and ( 29 a), computed for an NF number of values of reduced frequencies k;, may be 
combined as 


Q R (k 2 ) 


SR(k 2 ) 

Qte) 


site) 

QrI^nf- 1) 

. Ql(kNF-i). 


Sr^nF- 1) 

_ Si(kNF-i)_ 


A 2 

A3 

LA 4 


( 31 ) 


or 


and a least square solution 


Q = SA 


A = [S T S] l S T Q 


( 32 ) 


( 33 ) 


yields the required coefficients A 2 , A3, and A4. This procedure may be easily extended for a larger 
number of lag terms, if desired. Equation ( 27 ) may be rewritten as 


Mri + Cri + Kii + q A 0 Tl + Ai (^-|q + A 2 + A3X1 + A 4 X 2 + • • • 

and collecting like terms gives 

(k + <7 AqVi + C + ^^jAiJn + M + ^(yf A 2 ri +^A 3 X! +^a 4 x 2 + . 


= 0 


( 34 ) 


• - 0 


( 35 ) 


or 


/\ 


Also 


^+01+ Mil + q A3X1 + q A 4 X 2 + . • • = 0 


X- = 

Aj r 


( 36 ) 


s + 


KM 


from which 




( 37 ) 


( 38 ) 
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Equations (36), (37), and (38) can be rewritten as one set of matrix equations 


T 


0 10 0 


I 




✓N. /N 

/-X 



M 



__ 

-K — C -47A3 -q^A 



I 


x. 


0 I -^8,1 0 


Xi 

I . 


-X 2 - 


b 


_X 2 _ 



0 I o -^p 2 i 

b J 



m'x' = k'x* 

from which 

X’ - (MyiK'X' 

= RX' 

Also, the state-space vector X' may be rearranged as 

x"=[(rir Tie Tlr Xi X 2 ) (n 5 T1 »)] 

= [x u] 


( 39 ) 


( 40 ) 


( 41 ) 


( 42 ) 


and equation (41) partitioned as 


A 

X 


Ru 

Rm 


X 

- u - 


Rn,i 

Rn.n_ 


- U - 


where the first set of matrix equations denotes the plant dynamics, and the second set represents the 
dynamics of control modes. In the case of plant dynamics, the state-space equations become 


( 43 ) 


X = AX + Bu ( 44 ) 

the relevant matrices and vectors being defined as 

✓s 

A plant dynamics matrix 

B control surface influence matrix 
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X generalized coordinates in inertial frame 
u control surface motion input into plant 

and in which the terms AX and Bu represent for an aircraft, for example, the airplane dynamics and 
forcing function on airplane due to control surface motion, respectively. 

Coordinate Transformation 

To incorporate control laws and feedback, it is necessary to transform equation (44) from the earth- 
fixed (inertial) to the body-fixed coordinate system. Since no transformations are applied to elastic and 
aerodynamic lag state vectors, a transformation of the form 

X = T 2 ~HaTi - T 3 ) X +T 2 _ 1 Bu 
= AX + Bu (45) 


in which 


and so forth, Ti being the 12 by 12 coordinate transformation matrix, yields the required state-space 
equation in the body coordinate system. 

Determination of Sensor Outputs 

The structural nodal displacements are related to the generalized coordinates by equation (26), and 
the related sensor motion can be expressed as 


qs = T s a>ri 


= C 0 X (46) 

where Co = [T s O 0 0 0], and in which T s is an interpolation matrix. Similar relations may be derived 
for sensor velocities and acceleration as 


where 


V 


T s d>ii 

Ms. 


-T s a>ri_ 


= CiX 


T s O 0 0 

0 T s O 0 


0 

0 


( 47 ) 
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Equation (45) is next premultiplied by Ci to yield 


CiX = CiAX + CiBu 


= C2X + D211 


and adjoining equations (46) and (48), the following expression is obtained 


y = 


q s 

q s 

As. 


Co 

C2J 


x + 


0 

LD2J 


u 


( 48 ) 


or 


y = CX + Du (49) 

which is the required sensor output relationship, the matrices C and D signifying output at sensor due to 
body and control surface motions, respectively. 

Augmentation of Analog Elements and Controller 

The complete state-space formulation for an aircraft incorporating structural and aeroelastic effects 
is represented by equations (45) and (49), derived earlier. To conduct an aeroservoelastic analysis, it is 
essential to augment such a formulation with associated analog elements like actuators, sensors, notch 
filters, and prefilters along with the controller. Thus denoting the state-space equations of one such 
element as 


x«> = A®X® + 

(50) 

y(i) = C('>X (i > + D®U® 

(51) 


these can be augmented to such original equations (45) and (49) as appropriate; thus, typically, for the 
case of a connection from plant output to the external input, the relevant formulation is as follows: 


or 


■ 

X 

_[ A 0 

X 

+ 

B 

Lx (i, J 

Ib®C A (i) . 

lx®J 


B (i) D. 


[u] 


(52) 


X(i) = A(i)X(i) + B(i)u 


(53) 


noting that uO) = y. Also 
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or 


y 


c 0 

X 

-f- 

D 

y® J 


.D®C C® - 

. X® . 


-D®D. 


y(i)= C (i) X (i) + D (i) u (M) 

becomes the new sensor output expression. 

Any analog element, including a controller, can be augmented in a similar manner. Figure 21 shows 
a typical feedback control system. For such a system, the three sets of relevant matrix equations are 


X = AX + Bu 

(55) 

y = CX + Du 

(55a) 

u = r ~ Gy 

(55b) 

where equation (55b) is the feedback equation. By applying Laplace transformations to equa- 
tion (55), the following relationship is obtained: 

sX(s) = AX(s) + Bu(s) 

(56) 

y(s) = CX(s) + Du(s) 

(56a) 

u(s) = r(s) - G(s)y(s) 

(56b) 

Further, from equation (56) 

X(s) = [si - A]-'Bu(s) 

(57) 


and substitution of equation (57) into equation (56a) yields the required open-loop frequency response 
relationship 

y(s)= [c(sI-Ar'B +D]u(s) 

= H(s)u(s) (58) 

H(s) being the equivalent open-loop (loop-gain) transfer function with the analog controller or the open- 
loop transfer function without the controller. To obtain the closed-loop frequency response relationship, 
equation (58) is first substituted in equation (56b), resulting in 

u(s) = r(s) - G(s)H(s)u(s) (59) 

or 

u(s) = [1 + G(s)H(s)] _1 r(s) (59a) 
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and again substitution of equation (58) yields 


y(s) = (H(s)[I + G(s)H(s)]-‘)r(s) (59b) 

= H(s) r(s) (60) 

in which H(s) is the desired closed-loop transfer function. The frequency responses plots can be simply 

obtained from the transfer matrices H(s) orH(s), as the case may be. Associated damping and 
frequency values for the system, for the loop-gain or open-loop case, may also be calculated by solving 
the eigenvalue problem of the relevant A matrix for various ki values and observing the changes in sign 
of the real part of an eigenvalue. 

In the presence of a digital controller, a hybrid approach (ref. 6) is adopted for the frequency 
response solution. Thus, if A', B\ C\ and D* are the state-space matrices associated with the controller, 
the related transfer function is simply given by 

G(z) = C'[zl - A’]-^' + D' (61) 

and the frequency response relationship for the hybrid analog/digital system can be written as 

(621 

= H*(s)u(s) (62a) 


in which 

H(s) is the open-loop transfer function for the plant and other analog elements 

[ZOH] is the zero order hold complex expression (= e~ ST 

and where H*(s) is now the equivalent open-loop (loop-gain) transfer function of the hybrid system. 

The closed-loop frequency response relationship may be obtained as before by using equations (62a) and 
(56b) 


y(s) = {H(s) [I + G(s) H(s) 1 } r(s) 

= H*(s)r(s) (63) 

To compute the damping and frequencies, the analog plant dynamics matrix A is first transformed into 
the z-plane by the standard discretization procedure which is next augmented to the A* matrix. The 
appropriate eigenproblem solution of the final matrix yields the required results, as before. 

The STARS program has been extended to include capabilities representative of formulations 
presented in this section. 
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6. DATA INPUT PROCEDURE 


(STARS-AERO AND ASE) 


Figure 22 depicts the data input strategy for the entire ASE analysis procedure; such input for the 
solids module is described in section 3. In the following, the data pertaining to the other related analyses 
are given in the appropriate order, in which AERO module data input is compatible with the program 
described in reference 16. 



10154 


Figure 22. ASE analysis data input scheme. 


79 






6 . 1.1 


6.1 GENMASS Data 


(STARS-AERO-GENMASS) 


$ JOB DESCRIPTION 
Format (FREE) 

6.1.2 ISTMN, NLVN, GR 
Format (315) 

1. Description: Generalized mass matrix and modal interpretation generation data. 

2. Notes: 

ISTMN = integer specifying starting mode number 
NLVN = number of laterally vibrating nodes 

GR = gravitational constant 

6.1.3 $ LATERALLY MOVING NODAL NUMBERS DATA (Required if NLVN > 0) 

Format (FREE) 

6.1. 3.1 ( LN(I), 1 = 1, NLVN ) 

Format (15) 

1. Description: NLVN number of nodes input data. 

2. Notes: 

Input of a GR value is needed to convert generalized mass data into generalized weight 
acceptable to AERO module. 

If GRIDCHG is used, then the LN refers to STARS nodes (that is, the input vector, as 
defined in section 3.2.2 of the STARS manual). 

If direct STARS interpolation is used, then the LN refers to nodes as defined in STARS (that 
is, the output vector, as defined in section 3.5.10 of the STARS manual). 

3. Note: 

The data is to be stored in the file GENMASS. DAT. 
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6.2 GRIDCHG Data (STARS-AERO-GRIDCHG) 

6.2. 1.1 $ JOB TITLE 
Format (FREE) 

6.2. 1.2 NELN, NLINES, NOSURF 
Format (315) 

6.2. 1.3 IDELE, NMOD 
Format (215) 

6.2. 1.4 NBLOCK, IRPEAT 
Format (215) 

1 . Description: General input data. 

2. Notes: 

NELN = number of eliminated nodes from input vector 

NLINES = number of output vector interpolation lines 
0 < NLINES < 20 

NOSURF = number of sets of input vector coordinates to be translated 

IDELE = flag for deletion of interpolation elements 

= 0, for no elimination of interpolation element(s) 

= 1, to eliminate interpolation element(s) 

NMOD = number of output points whose values are to be changed to a user- specified 
value (for all modes) 

NBLOCK = number of blocks of added deflections 

IRPEAT = flag for reuse of deflections for different modes 
= 0, for user to input all blocks for all output modes 
= 1, to repeat first subset of block data for all subsequent modes 

6.2. 2.1 $ ELIMINATED INPUT NODES (Required if NELN * 0) 

Format (FREE) 

6.2.2.2 ( NODEL(I), I = 1, NELN ) 

Format (15) 

1. Description: NELN indices of nodes in input vector whose deflections are not used in 

interpolation. 
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2. Notes: 


NODEL(I) = node in input data that is not to be used in interpolation. 

6.2.3. 1 $ NUMBER OF POINTS ON OUTPUT VECTOR LINES 
Format (FREE) 

6.2.3.2 ( NGP(I), 1=1, NLINES ) 

Format (15) 

1. Description: NLINES sets of numbers of output points. Each set makes up part of the output 

vector. 

2. Notes: 

NGP(I) = number of points that will be interpolated to on each line 
0 < NGP(I) < 12 

6.2.4. 1 $ ENDPOINTS OF OUTPUT VECTOR LINES 
Format (FREE) 

6.2.4.2 ( (XTERMl(I), YTERMl(I), XTERM2(I), YTERM2(I), XT(I), YT(I)), 1=1, NLINES ) 
Format (6E10.4) 

1. Description: NLINES sets of endpoints for output vector interpolation lines and optional 

translations. 

2. Notes: 

XTERM 1 (I) , YTERM 1 (I) = inboard coordinates of line I 

XTERM2(I), YTERM2(I) = outboard coordinates of line I 

XT(I), YT(I) = optional translations to be applied to line data in X* and 

Y-directions 

6.2.5. 1 $ SPANWISE COORDINATES OF POINTS ON OUTPUT VECTOR LINES 
Format (FREE) 

6.2.5.2 ( ( YGP(J,I), J = 1, NGP(I) ), I = 1, NLINES ) 

Format (7E10.4) 

1. Description: NLINES sets of spanwise coordinate of desired point in output vector. 
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2. Notes: 


YGP(J,I) = spanwise coordinate of a point desired on an interpolation line, before any 

translation; translation, as defined in section 4.1, will automatically be applied 

6.2.6.1 $ TRANSLATION DATA FOR INPUT VECTOR POINTS (Required if NOSURF * 0) 

Format (FREE) 

6.2. 6.2 XTRAN, YTRAN, ZTRAN 
Format (3E10.4) 

6.2.63 NODNUM 
Format (15) 

1. Description: NOSURF subsets of input vector nodal data. 

2. Notes: 


XTRAN = value to be added to X-coordinate of input vector in set 

YTRAN = value to be added to Y-coordinate of input vector in set 

ZTRAN = value to be added to Z-coordinate of input vector in set 

NODNUM = index of node to be translated 

A data set is terminated if NODNUM is read as -1; a node should not be referenced more 
than once. 


6.2.7. 1 $ INTERPOLATION ELEMENT DATA 
Format (FREE) 

6 . 2.12 NXPT 

Format (15) 

6.2.7.2.1 ( XMESH(I), 1=1, NXPT ) 

Format (7E10.4) 

6.2.73 NYPT 

Format (15) 

6.2.73.1 ( YMESH(I), 1=1, NYPT ) 

Format (7E10.4) 

1. Description: Streamwise and spanwise finite element interpolation boundaries. 
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2. Notes: 


NXPT 


number of stations in X-direction for interpolation grid 
2 < NXPT < 20 


XMESH(I) = actual X-coordinates of streamwise stations in interpolation grid, in 
ascending order 

NYPT = number of stations in Y-direction for interpolation grid 
2 < NYPT < 20 

YMESH(I) = actual Y-coordinates of spanwise stations in interpolation grid, in 
ascending order 

6.2. 8.1 $ INTERPOLATION ELEMENT DELETION DATA (Required if I DELE > 0) 

Format (FREE) 

6.2.8.2 IOPT 
Format (15) 

6.2. 8.3 NCOL, NROW 
NCOL 
NROW 
Format (215) 

1. Description: Data for elimination of finite element interpolation elements. 

2. Notes: 

IOPT = type of elimination 

= 0, to proceed to next set to be eliminated 
= 1, to eliminate following element(s) 

= 2, to eliminate following row of elements 
= 3, to eliminate following column of elements 
= 4, to quit all eliminations 

NCOL = column of interpolation element(s) 

NROW = row of interpolation element(s) 

6.2.9. 1 $ OUTPUT VECTOR MODIFICATION DATA (Required if NMOD > 0) 

Format (FREE) 

6.2.9.2 ( NODE(I), DEFL(NODE(I)), 1=1, NMOD ) 

Format (15, E10.4) 

1 . Description: Sets a deflection to a user input value (for all output modes), where number of 

output modes NTOTAL = NR + NCNTRL - ISTMN + 1. 


(Required if IOPT = 1 ); or 
(Required if IOPT = 2); or 
(Required if IOPT = 3) 
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2. Notes: 

NODE® = output point index 

DEFL(NODE(I)) = new deflection value 

NR = analytically calculated number of roots (from section 3. 1 .6 of the 

STARS manual) 

NCNTRL = number of rigid body control modes (from section 3, 1 . 3 of the 

STARS manual) 

ISTMN = integer specifying starting mode numbers 

6.2.10. 1 $ BLOCK SPECIFICATION OF ADDITIONAL DEFLECTION DATA (Required if 

Format (FREE) NBLOCK > 0) 

6.2.10.2 ( (IBLOCK(I), NADD(I), IBFORE(I)), I = 1, NBLOCK ) 

Format (315) 

1. Description: NBLOCK sets of description of additional deflections to be added to output vector. 

2. Notes: 

IBLOCK(I) - user’s identification number of an added output block of output points 
NADD(I) = number of points in block 

IBFORE(I) = index of existing point in front of which block is to be inserted 
Succeeding values of IBFORE should be greater than the previous ones. 

6.2. 1 1.1 $ DEFLECTION DATA SPECIFICATION FOR BLOCKS (Required if NBLOCK > 0) 

Format (FREE) 

6.2.11.2 ( (NNODE(J), DADD(J)), J = 1, NADD(I) ) 

Format (15, E10.4) 

1 . Description: Added deflection data for each block. 

2. Notes: 

NNODE(J) = index of added point in set 
DADD(J) = deflection of added point in set 
This is repeated for NTOTAL modes. 

If IRPEAT = 1, the same deflections are reused for all modes. 
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6.2. 12. 1 $ EIGENVALUE SPECIFICATION FOR CONTROL MODES (Required if N CNTRL > 0) 
Format (FREE) 

6.2.12.2 ( EIGADD(I), I = 1, N CNTRL ) 

Format (El 0.4) 

1. Description: Eigenvalues for rigid body control modes from STARS. 

2. Notes: 

EIGADD(I) = user-input eigenvalues, in rad/sec, for rigid body control modes 


6.2.13 NOTES ON PROGRAM USAGE 

GRIDCHG is a versatile interpolation program that may be used as an alternative to the preferred di- 
rect interpolation option defined in section 3.1.4 of the STARS-SOLIDS module. It is utilized to inter- 
polate deflections, obtained by a finite element code or ground vibration survey, into the straight line in- 
put points required by the aerodynamic module. Options for separate interpolation of different surfaces 
and for modification by the user of both the input and output vectors exist. 

Input vector: The input vector is a calculated or measured vector with six degrees of freedom read 
from the file FOR096 (if bandwidth minimization is used) or from FOR048 (if bandwidth minimization 
is not used). Both files are STARS binary files. GRIDCHG normally uses only the Z component of the 
vector for the interpolation. However, GRIDCHG does read the input file for the GENMASS program 
as part of its input, and if the variable NLVN is nonzero in that file, then it reads from the file those 
nodes of the input vector for which the Y deflection is to be used (that is, a vertical surface). The 
GENMASS. DAT file must always be present for GRIDCHG to run, even if NLVN is zero. 

Discrete element interpolation: The user defines a set of rectangular elements used for the interpola- 
tion. Each element uses the deflections within its boundaries for a surface fit, with the added stipulation 
that adjacent elements have identical displacements and slopes at edges. The achievable quality of in- 
terpolation is a function of number and distribution of input nodes. The output vector is obtained using a 
surface fit within a particular element. Separate surfaces need to be individually interpolated, and this is 
accomplished by letting the value of a row or column of interpolation elements or columns between the 
surfaces to be set to zero. If the projection of the surfaces in the X-Y plane overlap, the user has the op- 
tion of temporarily modifying the coordinates of input and output vectors to separate them, thereby al- 
lowing their individual interpolation. 

Output vector: The output vector occasionally needs to be modified and/or have data added to it. 
This can be implemented as required. 

Eigenvalues: The STARS-SOLIDS module contains an option which allows additions of user input 
eigenvectors to those analytically calculated. The eigenvalues for those modes are added here. 
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6.3 AERO Data 


(STARS-AERO) 


6.3. 1.1 JOB TITLE - 1:6 (six lines of title cards) 
Format (FREE) 

6. 3.2.1 (LC(I),I = 1,40) 

Format (1015) 


1. Description: 

2. Notes: 


LC(1) 


LC(2) 

LC(3) 


LC(4) 


LC(5) 


LC(6) 


LC(7) 

LC(8) 


Basic data parameters. 


= integer defining flutter and divergence solution algorithm 
= -1, p-k type of solution 
= 0, pressure calculations only 
= 1, k and state-space solutions 
= 2, divergence analysis 

= maximum number of vibration modes to be used in analysis 
0 < LC(2) < 20 

- number of lifting surfaces 

0 < LC(3) < 30, for doublet lattice method (DLM) or constant pressure method 
(CPM) 

= number of reduced velocities, VBO, used in analysis 
If LC(1) = -1, set LC(4) = 6 
If LC(1) = 0 or 1, set 1 < LC(4) < 30 
If LC(1) = 2, set LC(4) = 1 

LC(4) and LC(13) apply to the reduced velocities described in section 6. 3. 4. 2 
and section 6.3.4.4 

= number of air densities at which flutter and divergence solutions are to be 
found 

0 < LC(5) < 10 

If LC(1) = 0, set LC(5) = 0 

= print option for tested aerodynamic forces used to check aerodynamic force 
interpolation 
= 1, print 
= 0, no print 

= print option for aerodynamic pressures 
= 1, print data 
= 0, no print 

= print option for lift and moment coefficients 
= 1, print data 
= 0, no print 
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LC(9) 

LC(10) 

LC(ll) 

LC(12) 

LC(13) 

LC(14) 

LC(15) 

LC(16) 

LC(17) 

LC(18) 


= input frequency-independent additions to the aerodynamic matrix QBAR 
= 1, make additions 
= 0, no additions 

= print option for full set of interpolated generalized forces when used in k 
solutions 
= 1, print data 
= 0, no print 

= index of mode whose frequency is to be used in normalizing flutter 
determinant 

Frequency chosen must be nonzero 
Suggested index is 1 

= index defining flutter determinant formulation 
= 1, for nonzero frequencies [D = K -1 (M + Ae)] 

= 0, in presence of zero frequencies [D = (M + Ae) -1 K] 

K = generalized stiffness matrix 

M = generalized mass matrix 

Ae = aerodynamic force matrix 

IfLC(l) =0, set LC(12) =0 

= index defining interpolation of aerodynamic forces 
= 0, no interpolation, to compute at each input VBO 
= 1 , to compute directly at only 6 VBOs, interpolate to others 
If LC(1) = -1, set LC(13) = 1 
If LC(1) = 0 or 2, set LC(13) = 0 
If LC(1) = 1, set LC(13) = 0 or 1, as desired 

= not used. Set = 0 

= index defining velocity scale in flutter solution output 
= 1, use true airspeed, TAS 
= 0, use equivalent airspeed, EAS 

= index defining addition of structural damping to complex stiffness matrix 
= 1, add a single damping value to all modes 
= -1, add an individual damping value to each mode 
= 0, no damping added 

= print option to display number of iterations required to find each root in a 
p-k solution 
= 1, print 
= 0, no print 

= option for root extrapolation in a p-k solution 

= 1, use root values at two previous velocities for initial estimation of a root 
as 0, use root value at previous velocity as root estimate 



If LC(1) *—1, set LC(18) =0 


LC(19) 

LC(20) 

LC(21) 

LC(22) 

LC(23) 

LC(24) 

LC(25) 

LC(26) 

LC(27) 

LC(28) 


= option for ordering of roots after a p-k solution 
= 1 , to perform ordering 
= 0, no ordering required 
If LC(1) 1, set LC(19) = 0 

= print option for iterated roots in p-k analysis or intermediate results in 
k analysis 
= 1, print 
= 0, no print 

= index for aerodynamics 

= 1, use doublet lattice method or constant pressure method (subsonic and 
supersonic Mach numbers, respectively) 

= index defining generation and storage of aerodynamic influence coefficients 
matrix 

= 0, compute and save 
= 1, read precomputed values from a file 

= print option for input modal vector 
= 1, print 
= 0, no print 

= print option for interpolated deflections and slopes of aerodynamic elements 
- 1, print 
= 0, no print 

= number of modal elimination cycles 
0 < LC(25) < 25 

= index defining additional flutter analysis 
= 0, no additional cycles 

> 0, perform additional flutter analysis cycles with stiffness variations applied to 
a mode 

0 < LC(26) < 20 

= index of mode whose frequency and stiffness is to be varied for the 
LC(26) cycles 

If LC(26) = 0, set LC(27) = 0 

= print option for modal eigenvectors 
= 1, print 
= 0, no print 

If LC(1) = -1, the eigenvectors for the critical flutter root in a user-chosen 
velocity interval are displayed 
If LC(1) = 0 or 2, set LC(28) = 0 


89 



LC(29) 

LC(30) 


LC(31) 

LC(32) 

LC(33) 


LC(34) 

LC(35) 

LC(36) 

LC(37) 


LC(38) 

LC(39) 

LC(40) 

6.3.3. 1 INV 

Format (15) 


If LC(1) = 1, the eigenvectors for all roots between user chosen reduced 
velocities, VBO, and real frequencies are displayed 

= print option for physical vectors corresponding to modal eigenvectors 
= 1, print 
= 0, no print 

= print option for k solution flutter determinant matrix analysis 
= 1, print 
= 0, no print 

If LC(1) = -1 or 0, set LC(30) = 0 

= index defining revisions to generalized mass matrix and modal frequencies 
= 1, revise 
= 0, no change 

= index defining revisions to generalized stiffness matrix 
= 1, revise 
= 0, no change 

= index defining type of aerodynamics 
= 1, steady state 
= 0, oscillatory 

If LC(1) = 2, set LC(33) = 1 

- not used. Set = 0 

= not used. Set = 0 

= not used. Set = 0 

= print option for aerodynamic element geometric data associated with doublet 
lattice and constant pressure methods 
= 1, print 
= 0, no print 

If LC(21) * 1, set LC(37) = 0 

= tape unit for ASCII printout of generalized forces and associated information. 
Suggest LC(38) = 99 

= not used. Set = 0 

= not used. Set = 0 


1. Description: Input vibration data location flag. 
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2. Notes: 


INV = integer defining location of input vectors, modal frequencies, and generalized 

masses 

= 1, STARS binary file 
= 2, this input file 

6.3. 3. 1.1 NMDOF (Required if INV = 2) 

Format (15) 

1. Description: Input vector degrees of freedom. 

2. Notes: 

NMDOF = total number of modal degrees of freedom used to define an input mode shape 
0< NMDOF <200 

6. 3. 3. 1.2 ( QZ(I), 1=1, NMDOF ) (Required if INV = 2) 

Format (7E10.0) 

1. Description: LC(2) sets of NMDOF input deflections. 

2. Notes: 

QZ(I) = principal out-of-plane deflection at point I of input vector 

6.3.3. 1.3 NCARD (Required if INV = 2) 

Format (15) 

1 . Description: Mass matrix specifications. 

2. Notes: 

NCARD = Number of nonzero generalized mass matrix elements 

6.3.3. 1.4 I, J, WW(I,J) (Required if INV = 2) 

Format (215, E10.0) 

1. Description: NCARD sets of data specifying nonzero generalized mass matrix elements. 

2. Notes: 

I = row index of generalized mass matrix 

J = column index of generalized mass matrix 

WW(I,J) = generalized mass (weight) matrix value, lbw 

6.3.3. 1.5 ( OMG(I), I = LC(2) ) (Required if INV = 2) 

Format (7E10.0) 
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1. Description: LC(2) modal frequencies. 

2. Notes: 

OMG(I) = modal frequency in proper order, Hz 

6.3.4.1 BR, FMACH 

Format (2E10.4) 

1. Description: Reference values for aerodynamics. 

2. Notes: 


BR = reference semichord, in. 

FMACH = reference freestream Mach number 

If FMACH < 1.0, doublet lattice method is used 
If FMACH > 1.0, constant pressure method is used 

6.3.4.2 ( VBO(I), I = 1, LC(4) ) (Required if LC(1) - 1) 

Format (7F10.4) 

1 . Description: LC(4) reduced velocities. 

2. Notes: 

VBO(I) = reduced velocity (V/bco) for flutter-divergence analysis 

If aerodynamic interpolation is chosen, then aerodynamic forces will be interpolated at each 
of these VBO(I) values, using the values for RVBO input in section 6.3.4.4; if direct calcula- 
tion is used, the aerodynamic forces will be calculated at each of these reduced velocities. 

0 < LC(4) < 30 

6.3.4.3 NV, VI, DV (Required if LC(1) = -1) 

Format (15, 2F10.0) 

1. Description: Airspeed velocity specification for p-k analysis. 

2. Notes: 


NV = number of velocities used in initial analysis, knots 

1 < N V < 20 

VI = lowest velocity from which to start analysis, knots 

VI > 200, suggested 

DV - velocity increment to be summed to VI during initial analysis, knots 

DV > 250, suggested 
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6.3. 4.4 TOLI, ( RVBO(I), I = 1, 6 ) (Required if LC(1) =-l or LC(13) = 1) 

Format (7E10.0) 

1. Description: Aerodynamic forces interpolation data. 

2. Notes: 

TOLI = tolerance value used for testing the interpolation fit; a nominal value of 
1.0E-03 is recommended 

RVBO(I) = reduced velocity at which aerodynamic forces will be computed, to be used as 
part of the basis in interpolating forces at other reduced velocities 

If aerodynamic interpolation is used, the RVBOs should span the entire range of VBOs of 
section 6.3.4.2. 

For LC(1) = -1, use the following approximations: 

1. RVBO(l) < 1.69 * 12.0 * VI / (BR * WMAX), where 
WMAX = maximum modal frequency, rad/sec. 

2. RVBO(6) > 1.69 * 12.0 * VMAX / (BR * WMIN), where 
VMAX = VI + (NV - 1)*DV, and 

WMIN = minimum modal frequency, rad/sec. 

6.3.5.1 MADD, LADD, MSYM (Required if LC(31) - 1) 

Format (315) 

1. Description: Specifications for changes to mass matrix and modal frequencies. 

2. Notes: 

MADD = number of changes to mass matrix 

LADD = number of changes to modal frequencies 

MSYM = integer specifying symmetry of mass matrix modifications 
= 0, changes are symmetric 
= 1, changes are nonsymmetric 

6.3. 5. 1.1 I, J, WW(I, J) (Required if MADD > 0) 

Format (215, F10.0) 

1. Description: MADD changes to the mass matrix. 

2. Notes: 

I = row index of mass matrix element 
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J 


= column index of mass matrix element 


WW(I,J) = value to be substituted for existing element in mass matrix, lb 
If MS YM = 0, specify only changes to upper triangular elements. 

6.3.5. 1.2 I, OMG(I) (Required if IADD > 0) 

Format (15, F10.0) 

1. Description: IADD changes to modal frequencies. 

2. Notes: 

I = index of mode to be changed 

OMG(I) = new frequency to be substituted for old, Hz 

6.3. 5.2 GDD (Required if LC( 1 6) = 1 ) 

Format (El 0.4) 

1. Description: General structural damping factor. 

2. Notes: 

GDD = A single value for hysteretic damping to be applied to all modes; the 
imaginary term on the diagonal of the complex stiffness matrix will be 
multiplied by the term GDD 

6.3. 5.3 NCD (Required if LC(16) = -1) 

Format (15) 

1. Description: Integer specifying individual structural damping. 

2. Notes: 

NCD = number of individual modes for which hysteretic damping will be specified 

6.3. 5.3.1 ( I, GDP(I) ) (Required if LC(16) = -1 and NCD * 0) 

Format (15, El 0.0) 

1. Description: NCD individual structural damping values. 

2. Notes: 

I = mode index 

GDP(I) = hysteretic damping applied to mode I, as in section 6.3. 5. 3 
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6.3.6. 1 GMAX, GMIN, VMAX, FMAX (Required if LC(1) * 2) 

Format (4F10.0) 

1. Description: Maximum and minimum scales for V-g, V-f print plots. 

2. Notes: 

GMAX = maximum value of damping scales for V-g plots 

GMIN = minimum value of damping scale for V-g plots 

VMAX = maximum value of velocity scale for V-g and V-f plots, knots 
FMAX = minimum value of frequency scale for V-f plots, Hz 

6.3.7. 1 ( RHOR(I), I = 1, LC(5) ) (Required if LC(1) * 0) 

Format (7F10.0) 

1. Description: LC(5) values of air density ratios. 

2. Notes: 

RHOR(I) = density ratio with respect to sea level 
0 < RHOR(I) < 10 

A separate flutter and/or divergence analysis is performed at each density ratio in which the 
aerodynamic force matrix is multiplied by the square root of the density ratio. 

6.3.8. 1 NADDF, NSYM (Required if LC(9) = 1 ) 

Format (215) 

1. Description: Specifications for frequency-independent additions to aerodynamic matrix. 

2. Notes: 

NADDF = number of following additions to the flutter-determinant aerodynamic matrix 

NSYM = index defining symmetry of additions 

= 0, additions are symmetric. Input only upper triangular elements 
- 1, additions are not symmetric 

6.3. 8.1.1 I, J, DETAD(I, I) (Required if LC(9) = 1) 

Format (215, 2E10.0) 

1. Description: NADDF frequency-independent additions to aerodynamic matrix. 

2. Notes: 

I = row index of additions 
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J 


= column index of additions 


DETAD(I, J) = value of addition. DETAD(I,J) is a complex value 
Additions to the aerodynamic matrix QBAR are done in the following manner: 

DETADreaL DETADjmAG 

QBAR = QBAR + k 2 ~ + 1 * k » 

where k is the reduced frequency. 

6.3. 8.2 NADDS, NSYM (Required if LC(32) - 1) 

Format (215) 

1. Description: Specifications for changes to generalized stiffness matrix. 

2. Notes: 

NADDS = number of following changes to the stiffness matrix 

NSYM = index specifying symmetry of changes 

= 0, changes are symmetric (B(IJ) = B(J,I)) 

= 1, changes are not symmetric 

6.3. 8.2.1 I, J, B(I, J) (Required if LC(32) = 1) 

Format (215, 2E10.0) 

1 . Description: NADDS changes to stiffness matrix. 

2. Notes: 

I = row index of changes 

J = column index of changes 

B(I, J) = new value of complex stiffness matrix element 

If NSYM = 0, only the upper triangular elements are input. 

6-3. 8.3 RATOM(I) (Required if LC(26) > 0) 

Format (7E10.0) 

1. Description: LC(26) values of stiffness variations for an input mode. 

2. Notes: 

RATOM(I) = ratio of modal frequency with respect to the original input value, OMG(I) 
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(Required if LC(25) * 0) 


6.3.83. 1 NOTIR, ( NINZ(J), J=l, NOTIR ) 

Format (1015) 

1. Description: LC(25) sets of modal elimination specification for flutter and divergence analysis. 

2. Notes: 

NOTIR = number of deleted modes in a given modal elimination cycle 

NINZ = index of individual deleted mode for a given cycle 

It should be noted that the aero module always does an initial analysis without modal 
deletions before doing any modal elimination analyses as defined in this section. 

6.3.9. 1 VA, VB (Required if LC(28) = 1 and LC(1) * 2) 

Format (2E10.0) 

1. Description: Eigenvector calculation range. 

2. Notes: 

VA = lower bound of the range over which the eigenvectors are to be calculated 

VB = upper bound of the range over which the eigenvectors are to be calculated 

If LC(1) = -1, the range is over velocity, V, knots 

V 

If LC(1) = 1, the range is over reduced velocity, 

63.9.2 FLO, FHI (Required if LC(28) = 1 and LC( 1 ) = 1 ) 

Format (2E10.0) 

1. Description: Eigenvector display range. 

2. Notes: 

FLO = lower bound of the frequency range over which the eigenvectors are to be 
displayed, Hz 

FHI = upper bound of the frequency range over which the eigenvectors are to be 

displayed, Hz 

63.10.1 FL, ACAP 
Format (2F10.0) 

1. Description: Reference length and area. 
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2. Notes: 


FL = reference chord of model, in. (2.0 * BR, normally) 

ACAP = reference area of the model, in^ 

6.3.10.2 NDELT, NP, NB, NCORE, N3, N4, N7 
Format (715) 

1 . Description: Doublet lattice and constant pressure methods geometrical paneling data. 

2. Notes: 


NDELT 


= index defining aerodynamic symmetry 
- 1, aerodynamics are symmetrical about Y = 0 
= -1, aerodynamics are antisymmetrical about Y = 0 
= 0, no symmetry about Y = 0 (single surface only) 


NP 


= total number of "panels" on all lifting surfaces 


NB 


= body identification flag 
= 0, no bodies of any kind 

> 0, number of slender bodies used for doublet lattice analysis 
= -1, constant pressure method body elements exist 
0 < NB < 20 for doublet lattice method 


NCORE = problem size, N*M, where 

N = total number of aerodynamic elements, and 
M = number of modes 


N3 


N4 


N7 


= print option for pressure influence coefficients 
= 1, print 
= 0, no print 

= print option for influence coefficients relating downwash on lifting surfaces to 
body element pressures 
= 1, print 
= 0, no print 

= index specifying calculation of pressures and generalized forces 
= 1, calculate 

= 0, cease computations after influence coefficients are determined 
If LC(1) = — 1 or 1, set N7 = 1 


6.3.11.1 IBOD1, IBOD2 
Format (215) 


(Required if NB = -1) 


1. Description: Aerodynamic elements defining contiguous panels which describe a supersonic 

body for the constant pressure method. 
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2. Notes: 


IBOD1 = first aerodynamic element on first panel (lowest index) 

IBOD2 = last aerodynamic element on last panel (highest index) 

6.3.12.1 6.3.12.1.1 to 6.3.12.1.5 are repeated for NP sets of surface paneling data. 

6.3.12.1.1 XO, YO, ZO, GGMAS 
Format (4F10.0) 

6.3.12.1.2 XI, X2, X3, X4, Yl, Y2 
Format (6F10.0) 

6.3.12.1.3 Zl, Z2, NEBS, NEBC, COEFF 
Format (2F10.0, IX, 213, 3X, F10.0) 

6.3.12.1.4 (THa), 1=1, NEBC) 

Format (6F10.0) 

6.3.12.1.5 (TAU(I), 1=1, NEBS) 

Format (6F10.0) 


1. Description: NP sets of data defining aerodynamic panels and their component aerodynamic 

elements. Section 6.3.12.1.1 translates and rotates panels. Such coordinates are in 
the global (aircraft) system indicating position of the origin of the LCS for each 
panel. Section 6.3.12.1.2 contains coordinates of points defining an aero- 
dynamic panel, while section 6.3.12.1.3 defines boundaries of "aerodynamic ele- 
ments" in the panel. The panel is divided into a number of smaller trapezoids, 
called "aerodynamic elements," by lines of constant percent panel chord and of 
constant percent panel span. Section 6.3.12.1.4 defines chordwise panel stations, 
and 6.3.12.1.5 defines spanwise panel stations. 

2. Notes: 


XO = translational value to be applied to x-coordinates, in. 

YO = translational value to be applied to y-coordinates, in. 

ZO = translational value to be applied to z-coordinates, in. 


GGMAS = panel dihedral or rotation, deg, about global x-axis 

GGMAS is in a right-handed coordinate system; an upright panel would require a positive 
rotation of 90°. 


XI = x-coordinate of panel inboard leading edge, in. 
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X2 


- x-coordinate of panel inboard trailing edge, in. 

X3 = x-coordinate of panel outboard leading edge, in. 

X4 = x-coordinate of panel outboard trailing edge, in. 

Y1 = y-coordinate of panel inboard edge, in. 

Y2 = y-coordinate of panel outboard edge, in. 

Z1 = z-coordinate of panel inboard edge, in. 

Z2 = z-coordinate of panel outboard edge, in. 

Coordinates are in the local coordinate system. 

NEBS = number of element boundaries in the spanwise direction 
2 < NEBS < 50 

NEBS must be set = 2 for each body interference panel 

NEBC = number of element boundaries in the chordwise direction 
2 < NEBC < 50 

COEFF = entered as zero 

TH(I) = chordwise element boundaries for the panel in fraction of chord 

0.0<TH< 1.0 

( TH(1) = 0.0, TH(NEBC) = 1.0 ) 

TAU(I) = spanwise element boundaries for the panel in fraction of span 
0.0 < TAU <1.0 

( TAU(l) = 0.0, TAU(NEBS) = 1.0 ) 

The data is to be repeated NP times in the following sequence: 

1 . Vertical panels or plane of symmetry (y = 0). 

2. Panels on other surfaces. 

3. Body interference panels. These panels must be one element wide (that is, 

NEBS = 2). 

There are (NEBS - 1) * (NEBC - 1) aerodynamic elements on a primary or control surface. 

Indices for aerodynamic elements start at the inboard leading edge element, increase while 
traveling aft down a strip, then outward strip by strip, ending at the outboard trailing edge 
element. 

The sum of all elements for all panels and bodies must not exceed 400. 

The total number of strips must not exceed 150. 
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6.3. 13. 1 The following is repeated for NB sets of slender body data. (Required if NB > 0) 

6.3.13.1.1 XBO, YBO, ZBO 
Format (3F10.0) 

6.3.13.1.2 ZSC, YSC, NF, NZ, NY, COEFF, MRK1, MRK2 
Format (2F10.0, IX, 312, 3X, 1F10.0, 213) 

6.3.13.1.3 (F(I), 1 = 1, NF ) 

Format (6F10.0) 

6.3.13.1.4 ( RAD(I), 1=1, NF ) 

Format (6F10.0) 

1. Description: NB sets of data defining subsonic slender bodies and their component elements. 

Section 6.3.13.1.1 defines X, Y, and Z global reference coordinates, and section 
6.3.13.1.2 defines slender body origin, elements, and any related interference 
panels. Section 6.3.13.1.3 defines slender body element stations, while section 

6.3.13.1.4 defines slender body radii. 

2. Notes: 

XBO = translational value to be added to X-coordinate, in. 

YBO = translational value to be added to Y-coordinate, in. 

ZB O = translational value to be added to Z-coordinate, in. 

ZSC = local z-coordinate of the body axis, in. 

YSC = local y-coordinate of the body axis, in. 

NF = number of slender body element boundaries along its axis 

2 < NF < 50 

NZ = flag for body vibration in z-direction 

= 1, body vibrating 
= 0, body not vibrating 

NY = flag for body vibration in y-direction 

= 1, body vibrating 
= 0, body not vibrating 

COEFF = entered as 0.0 

MRK1 = index of the first aerodynamic element on the first interference panel 
associated with this slender body 
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MRK2 = index of the last aerodynamic element on the first interference panel 
associated with this slender body 

F(I) = x-coordinate of body station defining a slender body element in local 

coordinates, in starting with body nose and proceeding aft 

RAD(I) = radii of body elements at the stations F(J), in. 

NZ must never equal NY. 

Vertically vibrating bodies should be input before laterally vibrating ones; if both vertical and 
lateral body vibrations are desired in a single body, two bodies are input at the same location 
with corresponding NZ and NY. 

A slender body, as defined here, is a frustum of a right angle cone; there are (NF - 1 ) slender 
body elements. 

6.3.14.1 NSTRIP, NPR1, JSPECS, NSV, NBV, NYAW 
Format (615) 

1. Description: General aerodynamics data. 

2. Notes: 


NSTRIP = number of chordwise strips of panel elements on all panels. 

For LC(8) = 0, set NSTRIP = 1 

Printouts of lift and moment coefficients for the strips occur for NSTRIP > 1 
Never set NSTRIP = 0 


NPR1 


JSPECS 


NSV 

NBV 

NYAW 


= print option for pressures in subroutines QUAS or FUTSOL. Use only for 
debugging 
= 1, print 
= 0, no print 

= index describing plane's aerodynamic symmetry about Z = 0 
= 1, antisymmetrical aerodynamics about Z = 0 (biplane or jet effect) 

= -1, symmetrical about Z = 0 (ground effect) 

= 0, no symmetry about plane Z = 0 

- number of strips lying on all vertical panels on the symmetric plane Y = 0 
= number of elements on all vertical panels lying on the plane Y = 0 
= symmetry flag 

= 0, if NDELT = 1 (symmetric about Y = 0) 

= 1, if NDELT = -1 (antisymmetric about Y = 0) 

= 0 or 1, if NDELT = 0 (asymmetric about Y = 0) 


102 



The total number of strips, summed over all panels, may not exceed 150 regardless of the 
value of NSTRIP. 

6.3.14.1.1 ( LIM(I,1), LIM(I,2), LIM(I,3), 1 = 1, NSTRIP ) 

Format (313) 

1. Description: NSTRIP sets of data defining chordwise strips for aerodynamic coefficient 

calculations. 

2. Notes: 

LIM(I,1) = index of first element on each chordwise strip 

LIM(I,2) = index of last element on each chordwise strip 

LIM(I,3) = 0 

For NSTRIP = 1, a blank card is used. 

6.3.15.1 6.3.15.1.1 to 6.3.15.1.2 are repeated for LC(3) sets of primary surface data. 

6.3.15.1.1 KSURF, NBOXS, NCS 
Format (1L5, 215) 

6.3.15.1.2 NLINES, NELAXS, NICH, NISP 
Format (415) 

6.3.15.2 6.3.15.2.1 and 6.3.15.2.2 are repeated for NLINES subsets of data. 

6.3.15.2.1 NGP, XTERM1, YTERM1, XTERM2, YTERM2 
Format (15, 4E10.0) 

6.3.15.2.2 ( YGP(I),I= 1,NGP) 

Format (8E10.0) 

6.3.15.3.1 DIST (Required if NELAXS = 1) 

Format (El 0.0) 

6.3.15.3.2 ( X1(I), Y1(I), X2(I), Y2(I), 1=1, NCS ) (Required if KSURF = T) 

Format (4E10.0) 

6.3.15.3.3 NLINES, NELAXS, NICH, NISP (Required if KSURF = T) 

Format (415) 

6.3. 15.4 6.3.15.4.1 and 6.3.15.4.2 are repeated for NLINES subsets of data. (Required if KSURF = T) 

6.3. 15.4. 1 NGP, XTERM1, YTERM1, XTERM2, YTERM2 
Format (15, 4E10.0) 
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6.3.15.4.2 ( YGP(I), I = 1, NGP ) 

Format (8E10.0) 

6.3.15.5 DIST (Required if NELAXS = 1 and KSURF = T) 

Format (El 0.0) 

1. Description: LC(3) sets of input modal vector data to be applied to interpolation of deflections 

for primary and control surface aerodynamic elements. 

2. Notes: 


KSURF 


NBOXS 

NCS 

NLINES 


= flag indicating control surfaces on a primary surface 
= T, this surface has one or more control surfaces with forward hinge lines 
= F, this surface has no control surfaces 

= number of elements on this surface, including those on control surfaces 

= number of control surfaces on primary surface 

0 < NCS < 5 

= number of lines along which input modal vector data are prescribed 

1 < NLINES < 20 


NELAXS = index defining input vector components 

= 1, translation and pitch rotation are prescribed at each input point 
= 0, only translation is prescribed 


NICH 


= index defining chordwise interpolation-extrapolation from input vector to 
aerodynamic elements 
- 0, linear 
= 1, quadratic 
= 2, cubic 


NISP 


= index defining spanwise interpolation-extrapolation from input vector to 
aerodynamic elements 
5 = 0, linear 
= 1, quadratic 
= 2, cubic 


NGP = number of points on an input vector line 
2 < NGP <12 


XTERM1 


X-coordinate specifying the inboard end of an input vector line in the local 
coordinate system 


YTERM1 


Y-coordinate specifying the inboard end of an input vector line in the local 
coordinate system 
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XTERM2 = X-coordinate specifying the outboard end of an input vector line in the local 
coordinate system 

YTERM2 = Y-coordinate specifying the outboard end of an input vector line in the local 
coordinate system 

YGP(I) - spanwise coordinate of a point along an input vector line, going inboard to 
outboard in the local coordinate system 

XI (I) = X-coordinate of the inboard terminus of the Ith control surface leading edge 

in LCS 

Y1(I) - Y-coordinate of the inboard terminus of the Ith control surface leading edge 

in LCS 

X2(I) = X-coordinate of the outboard terminus of the Ith control surface leading edge 

in LCS 

Y2(I) = Y-coordinate of the outboard terminus of the Ith control surface leading edge 

in LCS 

DIST = displacement reference distance 

6.3.16.1 The following sets of data are repeated NB times. (Required if N B > 0) 

6.3.16.1.1 NGP, NSTRBP, IPANEL 

Format (315) 

6.3.16.1.2 (XGP(I), I = 1, NGP) 

Format (6F10.0) 

1. Description: NB sets of data describing input modal vector to be applied to slender body 

aerodynamic elements deflection. 

2. Notes: 

NGP = number of points on a slender body axis at which input vector data are 
prescribed 
2 < NGP < 20 

NSTRIP as number of interference panels (or strips) associated with a slender body 

IPANEL = index of the first such interference panel associated with a slender body 

XGP(I) = streamwise coordinate of each point at which input modal data are prescribed, 

in LCS 

This data is not to be input for a constant pressure method model 
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6.3.17.1 KLUGLB 
Format (15) 

1. Description: Print option for global geometry. 

2. Notes: 

KLUGLB = print option for aerodynamic elements in global coordinate system 
= 1, print 
= 0, no print 

6.3. 18 NOTES ON PROGRAM USAGE 
Aerodynamic Modules 

The STARS aerodynamic module consists of two unsteady, linear, inviscid, aerodynamic codes: 
the doublet lattice method (DLM) for subsonic analyses, and the constant pressure method (CPM) 

(ref. 14) for supersonic analyses. Flutter and divergence solutions may be obtained by k, p-k, or state- 
space methods. 

Aerodynamic Modeling 

The aerodynamic elements on lifting and interfering surfaces consist of trapezoidal elements parallel 
to the free stream. The aspect ratio of an element should ideally be on the order of unity or less. 

The number of elements required for accurate analysis varies with the model and the reduced 
frequency values. Increasing the number of elements will increase the computational time. Higher 
reduced frequencies require smaller and, therefore, more elements. A guide for element size in the 
streamwise direction is 


A xk < 0.04, 

where k is reduced frequency, and Ax is element length. 

Elements should be concentrated near wing tips, leading and trailing edges, control surface hinges, 
and so forth. As a guide, a cosine distribution of elements over the wing's chord and full span may 
be adopted. 

The total number of aerodynamic elements, including any slender body elements, may not exceed 
400. Similarly, the total number of strips is limited to 150, whereas the number of modes in the analysis 
may not exceed 20. 

The surface element may be thought of as having an unsteady horseshoe vortex bound along the 
quarter chord of the element and trailing aft to infinity. The downwash from the unsteady vortices are 
calculated at a control point located at the three-quarter chord of an element's centerline. Since the 
induced downwash at the center of a vortex is infinite, no control point should ever lie on any vortex 
line, such as along the extension of any element edge, either upstream or downstream. 
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6.4 CONVERT Data 


(STARS-ASE-CONVERT) 


Purpose: Prepare CONVERT.DAT data file; selection of desired modes. 

Description: Enables selection of desired modes. 

6.4.1 $ JOB TITLE 
Format (FREE) 

6.4.2 NM 
Format (15) 

1. Description: General data. 

2. Notes: 

NM = total number of desired modes to form reduced generalized matrices 

6.4.3 $ MODAL SELECTION AND ORDERING 
Format (FREE) 

6.4.3. 1 IOLD, INEW 
Format (215) 

1. Description: Orders the modes to be used for ASE analysis, NM sets of data. 

2. Notes: 

IOLD = old modal number 
INEW = new modal number 

3. Notes: 

Output is the reduced generalized force matrix and is stored in PD.DAT file for subsequent 
input into the ASE module. 


107 



6.5 ASE PADE Data (STARS-ASE-PADE) 

Purpose: Prepare PADE.DAT data file. 

Description: Performs Pade curve fitting of unsteady aerodynamic forces and state-space matrix 

formulation. 

6.5.1 $ JOB TITLE 
Format (FREE) 

6.5.2 NRM, NEM, NCM, NG, NS, NK, NA, RHOR, VEL, CREF, IWNDT, NQD 
Format (FREE) 

1. Description: General input data. 

2. Notes: 

NRM = number of rigid body modes 
NEM = number of elastic modes 

NCM = number of control modes 

NG = number of gusts 

NS = number of sensors 

NK = number of sets of input data at discrete reduced frequencies 

NA = order of Pade equation 

0<NA<4 

RHOR = relative aerodynamic density with respect to sea level 

VEL = true airspeed, ft/sec 

CREF = reference chord, ft 

IWNDT = wind tunnel correction index 

= 0, uses formulation as in reference 13 

= 1, uses wind tunnel data to modify aerodynamic generalized force matrix as in 
reference 13 

NQD = number of velocities for flutter and divergence analysis, to be set to 0 for ASE 
analysis as in reference 13 

6.5.3 $ TENSION COEFFICIENTS 
Format (FREE) 
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6.5.3. 1 (BETA(I), 1 = 1, NA) 

Format (FREE) 

1. Description: Pade approximate's data. 

2. Notes: 

BETA(I) = tension coefficients 

6.5.4 $ GENERALIZED MASSES 
Format (FREE) 

6.5.4.1 ((GMASS(I, J), J = I, NM), I = 1, NM) 

Format (FREE) 

1 . Description: Generalized mass data, upper symmetric half, starting with diagonal element. 

2. Notes: 

NM = total number of modes 

= NRM + NEM + NCM 

GMASS(I) = generalized mass of mode I, slugs 

6.5.5 $ GENERALIZED DAMPING 
Format (FREE) 

6.5.5. 1 (DAMP(I), 1 = 1, NM) 

Format (FREE) 

1. Description: Generalized damping data. 

2. Notes: 

DAMP(I) = generalized damping applied to mode I 

6.5.6 $ NATURAL FREQUENCIES 
Format (FREE) 

6.5. 6.1 (OMEGA(I), I = 1, NM) 

Format (FREE) 

1. Description: Modal frequency data. 

2. Notes: 

OMEGA(I) = natural frequency of mode I, rad/sec 



6.5.7 $ VELOCITIES FOR FLUTTER AND DIVERGENCE ANALYSIS (Required if NQD > 1 ) 
Format (FREE) 

6.5.7. 1 (VEL(I), I = 1, NQD) (Required if NQD > 0) 

Format (FREE) 

1. Description: True airspeed data for flutter and divergence analysis, ft/sec. 

2. Notes: 

VEL(I) = airspeed at which to solve a matrix for frequency and damping 

6.5.8 $ AIRCRAFT ANGLES, DEGREES OF FREEDOM (Required if NQD = 0) 

Format (FREE) 

6.5.8. 1 PHI, THETA, PSI, US, VS, WS, PS, QS, RS, PHID, IHAD, PSID, NDOF 

Format (FREE) (Required if NQD = 0) 

1. Description: Data for transformation of earth- to body-centered coordinate systems. 

2. Notes: 

PHI 

THETA 
PSI 

US, VS,WS 
PS, QS, RS 
PHID, IHAD, PSID 
NDOF 

6.5.9 $ SENSOR DATA 
Format (FREE) 

6.5.9.1 IFLSI 
Format (FREE) 

1. Description: Flag for identification of sensor interpolation points in presence of GVS 

data only. 

2. Notes: 

IFLSI = 1, for antisymmetric case 


= roll angle, deg 
= pitch angle, deg 
= yaw angle, deg 
= body axes velocities 
= angular rates 
= Euler angle rates 

= number of aircraft degrees of freedom; a negative sign indicates 
antisymmetric case 

(Required if NQD = 0 and NS > 0) 
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= -1, for symmetric case 
= 0, for non-GVS case 


6.5. 9.2 XS, YS, ZS 

Format (FREE) 


(Required if NQD = 0 and NS > 0) 


LX, MY, NZ, THX, THY, THZ 
Format (FREE) 

1. Description: Sensor location and orientation; NS sets of data. 

2. Notes: 


XS = X-coordinate of sensor 

YS = Y-coordinate of sensor 

ZS = Z-coordinate of sensor 

LX = direction cosine for accelerometer normal in X 

MY = direction cosine for accelerometer normal in Y 

NZ = direction cosine for accelerometer normal in Z 

THX = direction cosine for pitch axis about X 

THY = direction cosine for pitch axis about Y 

THZ = direction cosine for pitch axis about Z 


3. Notes: 


For the case IFLSI * 0, the user must modify file VEC_AND_COORDS.DAT by defining 
appropriate sensor location. This is done by setting the fourth column of the relevant nodes 
in the nodal coordinates section of the data file to the appropriate value of IFLSI. 
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6.6 ASE FRESP Data 


(STARS-ASE-FRESP) 


Purpose: Prepare frequency response analysis data file. 

6.6.1 $ JOB TITLE 
Format (FREE) 

6.6.2 NX, NY, NU, NV, NXC, DELTAT, TDELAY, MAXBC, MAXPO 
Format (FREE) 

1. Description: System parameters. 

2. Notes: 

NX = number of states in the plant 

s= [2*(NRM + NEM) + NA*(NRM + NEM + NCM)] (Refer to section 6.5) 

NY = number of outputs from the plant 

= (number of rows of C matrix) 

= (2*NS*3) 

NU = number of inputs to the plant 

= (2*NCM) 

NV = number of external inputs to the system 

NXC = total number of continuous states (plant plus analog elements) 

DELTAT = sample time for digital elements 
TDELAY = system time delay 
MAXBC = maximum number of block connectivity 
MAXPO = maximum polynomial order plus one 

6.6.3 NB, NYBTUV, IADDRA, IADDCB, IADDRC, NLST, NDRESP, IRP, ITRP 
Format (FREE) 

1. Notes: 

NB = number of analog and digital elements in the system including the summing 

elements and excluding the plant 

NYBTUV = NYTOV + NBTOU (See section 6.6.10.2 for definitions.) 

IADDRA = additional rows of A due augmentation of control elements. Appropriate 

summation of orders of polynomial of all analog elements for open- as well as 
closed-loop solutions to be derived from block connectivity input. 
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IADDCB = additional columns of B due augmentation of control elements 

IADDRC = additional rows of C due augmentation of control elements, equal to the 
number of connecting links into block/connecting junction 

NLST = total number of frequency range specifications for frequency response 
computations 

NDRESP = number of times the loops are broken for open-loop response evaluation 
IRP = frequency response problem number to be evaluated 

ITRP = total number of frequency response cases 

6.6.4. 1 $ BLOCK CONNECTIVITY 
Format (FREE) 

1. Note: 

Analog blocks to precede digital blocks 

6.6.4.2 IBN, ICN1, ICN2, ICN3, IEXI, ISLPCL, IELPCL 
Format (715) 

1. Notes: 

IBN = integer defining block number 

ICN1, ICN2, 

ICN3 = connecting block numbers, up to 3 

IEXI = integer defining external input number 

ISLPCL = integer defining starting block of the closed-loop system 

IELPCL = integer defining closing block of the closed-loop system 

2. Notes: 

A symbolic gain block indicating closing of loop is identified by presence of starting and 
closing blocks. 

6.6.5. 1 $ TRANSFER FUNCTION DESCRIPTION, AS ORDER OF POLYNOMIALS, FOR EACH 
BLOCK 
Format (FREE) 
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1. Notes: 


The polynomial descriptions pertain to either analog or digital elements, as appropriate. 

6.6.52 IBN, ICNP, ICDP 
Format (315) 

1. Notes: 

ICNP = integer defining number of coefficients in the numerator polynomial 

ICDP = integer defining number of coefficients in the denominator polynomial 

6.6.6. 1 $ LISTING OF POLYNOMIAL COEFFICIENTS 
Format (FREE) 

6.6.6.2 IBN, (POLCON(I), 1=1, MAXPO) 

Format (15 , < MAXPO > (E10.4)) 

IBN, (POLCOD(I), 1=1, MAXPO) 

Format (15, <MAXPO > (E10.4)) 

1. Notes: 

1. The coefficients are to be listed in increasing order of polynomials. 

2. The numerator coefficients (POLCON) are placed in one row followed by the 
denominator (POLCOD) ones in the next row, for each block, one block at a time. 

3. Data to be prepared for each block, NB sets of data being the input. 

6.6.7. 1 $ GAIN INPUTS FOR EACH BLOCK 
Format (FREE) 

6.6.12 IBN, GAIN 

Format (5(15, E10.4)) 

1. Notes: 

Gains may alternatively be the input as multiplier of polynomial coefficients in the 

numerator. NB sets of data are the input. 

6.6.8. 1 $ SPECIFICATION FOR SYSTEM OUTPUTS, NYB = NY + NB NUMBER OF DATA 
Format (FREE) 


6.6.8.2 ISOl, IS02, . . ., ISONYB 
Format (1615) 

1. Description: This data is needed for closed-loop frequency response analysis only. 
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2. Notes: 


1. Plant output are numbered 1 through NY. 

2. Each block output is numbered as NY + IBN. 

3. Notes: 

ISOI = desired output from any sensor (corresponding row of C matrix for the plant) 
and any control element (augmented thereafter) 

6.6.9. 1 $ SPECIFICATION FOR SYSTEM INPUTS, NUV = NU + NV NUMBER OF DATA 
Format (FREE) 

6 . 6.92 ISI1, ISI2, . . ISINUV 
Format (1615) 

1. Notes: 

1. Plant input are numbered 1 through NU. 

2. Each block input is numbered as NU + IEXI. 

2 . Notes: 

ISII = plant input (corresponding column of B matrix for the plant) and external input 

6.6. 10.1 $ CONNECTION DETAILS FROM PLANT TO BLOCKS 
Format (FREE) 

6.6.10.2 NYTOV, NBTOU, NBTOK 
Format (315) 

1. Notes: 

NYTOV = number of connections from plant outputs to external inputs 

NBTOU = number of block outputs connected to plant inputs 

NBTOK = number of digital element outputs connected to analog element inputs 

6.6.10.3 IYTOV1, IYTOV2 
Format (215) 

1. Notes: 

IYTOV1 = row number of C matrix corresponding to output from plant to feedback 
control system 
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IYTOV2 = external input number which describes connection of plant output to control 
system 

2. Notes: 

Repeat NYTOV times, ISO to IEXI. 

6.6.10.4 IBTOU1, IBTOU2 
Format (215) 

1. Notes: 

IBTOU1 = block number to be connected to plant input 
IBTOU2 = column of B matrix to which block is connected 

2. Notes: 

Output NBTOU times, IBN to ISI. 

6.6.10.5 IBTOK1, IBTOK2 
Format (215) 

1. Notes: 

Output NBTOK times, IBN (ANALOG) to IBN (DIGITAL). 

6.6.1 1.1 $ FREQUENCY RANGE SPECIFICATION (Required if NLST * ()) 

Format (FREE) 

6.6.1 1.2 FREQI, FREQF, NFREQ 
Format (315) 

1. Notes: 

FREQI = initial frequency 

FREQF = final frequency 

NFREQ = number of frequencies within range, logarithmically spaced 

NFREQ = number of frequencies within range, logarithmically spaced 

2. Notes: 

Data to be repeated NLST times 
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6.6.12.1 $ LOOP DEFINITIONS 
Format (FREE) 

6.6.12.2 ILOOP, IPRINT 
Format (FREE) 

6.6.12.3 NBRAK1, NBRAK2 
Format (215) 

1. Notes: 


ILOOP 


= integer defining loop type 
= 0, for closed loop case 
= 1, for open loop case 


IPRINT = eigensolution print option for closed loop case 
= 0, prints eigenvalues only 
= 1, prints eigenvalues and vectors 


NBRAK1 = block having the output signal 


NBRAK2 


block having the input signal 


2. Notes: 


This data is needed for ILOOP = 1 case, only. 
Data of 6.6.12.3 to be repeated NDRESP times. 


(Required if ILOOP = 1 ) 
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7. SAMPLE PROBLEMS 

(STARS Integrated Aero-Structural-Control Systems Analysis) 


A simplified aircraft test model (ATM) is selected as a standard problem for the full spectrum of 
ASE analyses. In this section, the relevant data (fig. 22) for associated SOLIDS, AERO, and ASE 
modules are presented in detail. Each such data set is also followed by relevant output of results. The 
input data are prepared in accordance with procedures described in section 6. 

Three perfect rigid body modes (Y-translation, X-rotation roll, and Z-rotation yaw about center of 
gravity - <J> PR ) and two rigid control modes (aileron and rudder deflections - d> c ) are generated in this 
module along with eight elastic (<X> E ) and three usual rigid body modes (d> R ), of which the latter are ex- 
cluded from consideration as GENMASS data input. The perfect rigid body modes 0 PR are moved in 
the front through convert data input for subsequent ASE analysis (d> = d> PR + O e + ct> c ). 


7.1 ATM: Free Vibration Analysis (STARS-SOLIDS) 

These input data pertain to the free vibration analysis of the finite element model. The direct modal 
interpolation option is used for subsequent flutter and ASE analyses. 

The finite element model (fig. 23) of the symmetric half of the aircraft is utilized for the vibration 
analysis; only the typical antisymmetric case is presented here. Figure 24 shows a direct interpolation 
scheme for subsequent aeroelastic and aeroservoelastic analyses. 

STARS-SOLIDS input data: 


RERO TEST NODEL 

C 

C ANT ISVNNETRIC HALF NODEL 

C IINTP - 1, DIRECT INTERPOLATION OF flODRL DATA 

C 

C NCNTRL - 5, FIRST THREE TO GENERATE PERFECT RIGID BOOV NODES 
C V TRRNSLRT ION, ROLL AND VAU, PLUS AILERONS ANO RUDDER CONTROL 

C NODES. 

c iiminii/iuiuimiuuuun/ummiinnnum/nuujmiimi 


71, 

1 19, 

1, 1 , 

22, 5, 

0 , 

o , 


0 , 

0 




0 , 

0 , 

0 , o , 

0 , 0 , 

5, 

132, 


0 





1 , 

1 , 

0 , 0 , 

0 , 1 , 

0 , 

1 







2, 

0 , 

2, 0, 

1 









1 , 

15, 

0, 0 . 7500E*03, 0 .0000E+00 , 

0.0 







NODAL DATA 











1 

300.0000 

200.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

2 

312.5000 

200.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

3 

325.0000 

200.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

1 

337.5000 

200.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

5 

350.0000 

200.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

6 

287.5000 

150.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

7 

303.1250 

150.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

8 

318.7500 

150.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

9 

331.3750 

150.0000 

o.oooo 

0 

0 

0 

0 

0 

0 

0 

0 

10 

350.0000 

150.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

11 

335.3750 

119.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

12 

350.0000 

119.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

13 

275.0000 

100.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

11 

293.7500 

100.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

IS 

312.5000 

100.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

16 

331 .2500 

100.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

17 

332.2500 

100.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

18 

350.0000 

100.0000 

o.oooo 

0 

0 

0 

0 

0 

0 

0 

0 

19 

262.5000 

50.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

20 

281.3750 

50.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 

21 

306.2500 

50.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 
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22 

328.1250 

50.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 


23 

329.1250 

51.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 


21 

350.0000 

51.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 


25 

350.0000 

50.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 


26 

250.0000 

0.0000 

0.0000 

1 

0 

I 

0 

1 

0 

0 

0 


27 

275.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


28 

300.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


29 

325.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


30 

350.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


31 

200.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


32 

100.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


33 

150.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


31 

150.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


35 

100.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


36 

500.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


37 

50.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


38 

560.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


39 

0.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


10 

600.0000 

0.0000 

0.0000 

1 

0 

1 

0 

1 

0 

0 

0 


11 

520.0000 

0.0000 

50.0000 

1 

0 

1 

0 

1 

0 

0 

0 


12 

580.0000 

100.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 

0 


13 

555.0000 

0.0000 

50.0000 

1 

0 

1 

0 

1 

0 

0 

0 


11 

520.0000 

20.0000 

0.0000 

0 

0 

0 

0 

0 

0 

0 
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$ LINE ELEflENT BASIC PROPERTIES 


1 

1 .5000 

37.5000 

18.8000 

18.8000 

2 

0.5300 

3.8000 
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Figure 24. ATM antisymmetric case, direct-surface interpolation scheme. 
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STARS-SOLIDS analysis results: 


Table 15 depicts the results of the free vibration analysis. Figure 25 shows the eight elastic mode 
shapes, whereas the three perfect rigid body and the two control modes are shown in figure 26. In order 
to effect correct response from the controllers, the perfect rigid body and control modes need to be 
defined in the fashion shown in table 16. 

Table 15. AERO test model: Antisymmetric free vibration analysis results. 



Mode shape 

Eigenvalue 

Generalized 
mass, lb 

Mode shape 

SOLIDS AERO-ASE 

Hz 

rad/sec 

1 





— 

113.8 

Rigid body X-rotation 

2 

— 

— 

— 

2,384 

Rigid body Y-translation 

3 

— 

— 

— 

111.6 

Rigid body Z-rotation 

4 

1 

10.175 

63.931 

8.2 

Vertical fin first bending 

5 

2 

12.449 

78.217 

235.1 

Fuselage first bending 

6 

3 

14.632 

91.934 

44.62 

Wing first bending 

7 

4 

28.741 

180.584 

60.53 

Wing second bending 

8 

5 

29.810 

187.301 

204.3 

Fuselage second bending 

9 

6 

32.450 

203.890 

47.87 

Wing first torsion 

10 

7 

35.815 

225.030 

3.233 

Fin first torsion 

11 

8 

59.664 

374.880 

21.84 

Wing third bending 

12 

9 

— 

— 

2,534 

Rigid body Y-translation 

13 

10 

— 

— 

151,200 

Rigid body roll 

14 

11 

— 

— 

589,000 

Rigid body yaw at 275 in. 

15 

12 

— 

— 

128.60 

Flap deflection 

16 

13 

— 

... 

14.22 

Rudder deflection 


Table 16. Rigid body and control mode generation parameters. 


Motion 

Symmetric analysis 

Antisymmetric analysis 

X-translation 

1.0 in X 


Y-translation 


1.0 in Y 

Z-translation 

-1.0 in Z 


X-rotation 


-A in Z 

Y-rotation 

-Ain Z 


Z-rotation 


-Ain Y 

Flap 

-A in Z 


Aileron 


+A in Z 

Elevator 

-Ain Z 


Rudder 


-A in Y 


In the table, the term A is defined as 


dN~dA 
A ~ 12 
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where d^ and d A represent the coordinates of the node under consideration and the axis of rota 
tion, respectively. 


Mode 4, fin IB 



Mode 8, fuselage 2B 



Mode 5, fuselage 1 B 



Mode 9, wing IT 



Mode 6, wing IB 



Mode 7, wing 2B 


Mode 10, fin IT 





Mode 11, wing 3B 



Figure 25. ATM antisymmetric case, elastic (<I>e) mode shapes. 


(a) Rigid body mode, X-Y plane motion. 





(d) Control mode, rudder motion 



(e) Rigid body mode, X-rotation motion. 


Figure 26. ATM antisymmetric case, perfect rigid body (Opr) and control (<3>c) modes. 


7.2 ATM: Generalized mass analysis (STARS-AERO-GENMASS) 


This run is made by deleting the first three rigid body modes so that 

O — <I>£ + Oq 


STARS-AERO-GENMASS input data: 


$ AERO TEST tlODES, RHTISVfttlETRIC UERSION 
4 39 386.088 

t LRTERRLLV 110UING direct Interpolation output NODE NUfIBERS 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

M 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 


The calculated generalized mass is depicted in table 15. 
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7.3 ATM: Aeroelastic analysis 


(STARS-AERO) 


The input data used for eventual ASE response analysis are given in this section. These data also 
enable flutter and divergence analysis of the aircraft. For a k or p-k method of flutter solution, the 
number of reduced frequencies in the data is increased from 10 to 28, and rigid body and con- 
trol modes are eliminated. Figures 27 through 29 show the ATM relevant aerodynamic element 
arrangements. 

STARS-AERO input data - ASE analysis: 


AERO TEST NODEL - RHT ISVNNETRIC CRSE - NCNTRL - 5 
SET UP FOR RSE SOLUTION, 

DIRECT SURFACE INTERPOLATION, 

EIGHT ELASTIC NODES, PLUS FIUE ADDED RIGID BODV-CONTROL NODES; 

REUISED RIGID BODV V TRANSLATIONS, ROLL, VAU, PLUS AILERON AND RUDDER NODES. 
NACH NO.- 0.90 ALTITUDE: SEA LEUEL 

I 15 3 10 1 00000 

I 000000000 
I 000000000 


0 0 
1 

0 0 

0 0 

0 99 

0 0 


1 

38.89 

0.90 





11000.0 

1000.0 

100.0 

50.0 

10.0 

5.0 

0.667 

0.500 

0.25 




.10 

-.10 

1100. 

80.0 



1 .0 






77.78 

15000. 





-1 15 

1 1200 

0 0 

1 

90.0 



508.0 

580.0 

532.0 

580.0 

20.0 

80.0 

0.0 

0.0 

1 1 

0.0 



0.0 

0.3333 

0.6666 

1.0 



0.0 

0.3333 

0.6666 

1.0 






90.0 



532.0 

600.0 

510,0 

600.0 

60.0 

100.0 

0.0 

0.0 

2 5 

0.0 



0.0 

0.2353 

0,1705 

0.7059 

1 .0 


0.0 

1.0 


90.0 



580.0 

600.0 

580.0 

600.0 

20.0 

80.0 

0.0 

0.0 

1 2 

0.0 



0.0 

1.0 





0.0 

0.3333 

0.6666 

1 .0 



255.0 

350.0 

262.5 

350.0 

20.0 

50.0 

0.0 

0.0 

3 5 

0.0 



0.0 

0.25 

0.50 

0.75 

1 .0 


0.0 

0.5 

1.0 




262.5 

328.125 

275.0 

331 .25 

50.0 

100.0 

0.0 

0.0 

1 1 

0.0 



0.0 

0.3333 

0.6666 

1.0 



0.0 

0.31 

0.66 

1.0 



275.0 

331 .25 

287.5 

331.375 

100.0 

150.0 

0.0 

0.0 

1 1 

0.0 



0.0 

0.3333 

0.6666 

1 .0 



0.0 

0.31 

0.66 

1.0 



287.5 

350.0 

300.0 

350.0 

150.0 

200.0 

0.0 

0.0 

1 5 

0.0 



0.0 

0.25 

0.50 

0.75 

1 .0 


0.0 

0.31 

0.66 

1.0 



328. 125 

350 . 0 

331.25 

350.0 

50.0 

100.0 

0.0 

0.0 

1 2 

0.0 



0.0 

1.0 





0.0 

0.31 

0.66 

1.0 



331.25 

350.0 

331.375 

350.0 

100.0 

150.0 

0.0 

0.0 

1 2 

0.0 



0.0 

1.0 





0.0 

0.31 

0.66 

1.0 
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50.0 


520.0 

580.0 

0.0 

0.0 

0.0 

0.4167 

0.0 

0.5 

527.5 

580.0 

0.0 

0.0 

0.0 

0.4167 

0.0 

0.34 

580.0 

600.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.5 

580.0 

600.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.34 

- 5.0 

600.0 

- 20.0 

0.0 

0.0 

0.1074 

0.5868 

0.6777 

0.0 

1.0 

- 5.0 

600.0 

0.0 

20.0 

0.0 

0.1074 

0.5868 

0.6777 

0.0 

1.0 


527.5 580.0 

3 0.0 

1 .0 

t .0 

540.0 580.0 

3 0.0 

1 .0 

0.66 1.0 

580.0 600.0 

2 0.0 

1.0 

580.0 600.0 

2 0.0 

0.66 1.0 

- 5.0 600.0 

211 0.0 

0.2149 0.3223 

0.7769 0.8678 


- 5.0 600.0 

211 0.0 

0.2149 0.3223 

0.7769 0.8678 


20.0 

50.0 100.0 

20.0 50.0 

50.0 100.0 

0.0 20.0 

0.4298 0.5083 

1 .0000 

20.0 0.0 

0.4298 0.5083 

1.0000 


0.0 

0.0 

14 0 1 

- 15.000 

25.000 

65.000 

295.000 

335.000 

365.000 

605.000 

645.000 


0.0 

10.0 

20.0 

40.0 

40.0 

40.0 

20.0 

15.0 


1 0 

0 2 

16 1 


0,0 

76 95 


145.000 

205.000 

245.000 

425.000 

485.000 

565.000 

20.0 

20.0 

40.0 

30.0 

30.0 

20.0 


T 16 1 

4 0 11 


7 

502.0 


2.0 

542.0 

100.0 





2.0 

20.0 


40.0 

55.0 

70.0 

85.0 

100.0 

7 

542.0 


2.0 

560.0 

100.0 





2.0 

20.0 


40.0 

55.0 

70.0 

85.0 

100.0 

7 

578.0 


2.0 

578.0 

100.0 





2.0 

20.0 


40.0 

55.0 

70.0 

85.0 

100.0 

2 

598.0 


82.0 

598.0 

100.0 





82.0 

100.0 







580.0 

20.0 

580.0 

80.0 




2 

0 1 

1 







2 

582.0 


22.0 

582.0 

78.0 





22.0 

78.0 







2 

598.0 


22.0 

598.0 

78.0 





22.0 

78.0 







T 

44 1 








6 

0 1 

1 







11 

252.0 


2,0 

302.0 

200.0 





2.0 

25.0 


50.0 

70.0 

90.0 

110.0 

130.0 


150.0 

170.0 


185.0 

200.0 




It 

270.0 


2.0 

313.0 

200.0 





2.0 

25.0 


50.0 

70.0 

90.0 

110.0 

130.0 


150.0 

170.0 


185.0 

200.0 




1 1 

298.0 


2.0 

324.0 

200.0 





2.0 

25.0 


50.0 

70.0 

90.0 

110.0 

130.0 


150.0 

170.0 


185.0 

200.0 




11 

323.0 


2.0 

336.0 

200.0 





2.0 

25.0 


50.0 

70.0 

90.0 

110.0 

130.0 


150.0 170.0 185.0 200.0 

3 348.0 2.0 348.0 48.0 

2.0 25.0 48.0 

4 348.0 152.0 348.0 200.0 

152.0 170.0 185.0 200.0 

328.125 50.0 334.375 150.0 

2 0 11 


330.0 


52.0 

336.0 

148.0 


) 

70.0 


90.0 

110.0 130.0 

148.0 

348.0 


52.0 

348.0 

148.0 


) 

70.0 


90.0 

110.0 130.0 

148.0 
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T 15 I 

3 0 11 


6 

522.0 


20.0 

512.0 

100.0 




20.0 

10.0 


55.0 

70.0 

85.0 

100.0 

6 

550.0 


20.0 

562.0 

100.0 




20.0 

10.0 


55.0 

70.0 

85.0 

100.0 

6 

578.0 


20.0 

578.0 

100.0 




20.0 

10.0 


55.0 

70.0 

85.0 

100,0 

580 . 0 

20.0 


580.0 

100.0 



2 

0 1 

1 






6 

582.0 


20.0 

582.0 

100.0 




20.0 

10.0 


55.0 

70.0 

85.0 

100.0 

6 

598.0 


20.0 

598.0 

100.0 




20.0 

10.0 


55.0 

70.0 

85.0 

100.0 

12 

2 M 
0.0 

100.0 


150.0 

200.0 

250.0 

300.0 


350.0 

100.0 


150.0 

500.0 

560.0 

600.0 


0 
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Buttline, 

in. 


Waterline, 

in. 


Figure 27. ATM antisymmetric case, half aircraft aerodynamic boxes. 




0 100 200 300 400 500 600 


Fuselage station, in. 
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Figure 28. ATM antisymmetric case, aerodynamic panels. 




0 100 200 300 400 500 600 

Fuselage station, in. 

10161 


(a) Aerodynamic boxes. 


-15 645 



0 100 200 300 400 500 600 

Fuselage station, in. 

(b) Slender body definitions. 


Figure 29. ATM antisymmetric case. 


STARS-AERO input data - k-type flutter analysis: 

These data pertain only to changes required in the corresponding ASE analysis case and occur within 
the first 16 lines; such data are presented next. 


AERO TEST tlODEL - ANTISVflllETRtC CASE - I INTP - 1 

K-FLUTTER SOLUTION 

DIRECT SURFRCE INTERPOLATION 

STARS STRUCTURAL tlODEL, BVPRSS RIGID BODV NODES IN GENNASS 
NACH NO. - 0.90 ALTITUDE: SEA LEUEL 

in/n/miinuiinmuninnininuiiuimunninni 

110 3 26 1 0 0 0 0 0 

I 100000000 
1 000000000 
0000000 99 00 
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38 . 89 
.350 

0.90 
. 745 

.940 

1.491 

1 .615 

1 .698 

1.864 

2.000 

2.450 

2.750 

3.150 

3.270 

3.400 

3.850 

4.150 

4.551 

5.250 

7.000 

9.000 

11.110 

15.000 

19.000 

24.070 

50.000 

140.000 

315.774 

616.746 

1200.000 


STARS -AERO analysis results: 

Table 17 provides the results of flutter analysis by various procedures using direct interpolation of 
modal data. The flutter solution based on the ASE method utilizing state-space formulation employs a 
data file as in section 6.5. Figures 30 through 32 depict the pattern of root location as a function of 
velocity for the k, p-k, and ASE methods. In this connection it may be noted that only the elastic modes 
are considered in these analyses. In the ASE method, the real (a) and imaginary (b) parts of the 
eigenvalues, termed as damping and frequencies, of the state-space plant dynamics matrix (A) are 
plotted against the air speed. In the k and p-k methods, the damping term is expressed as g’ = 
2 aco ( i/co n 2 , where c% is the relevant natural frequency [C 04 = (a^ + p 2 )l/ 2 ] 


Table 17. AERO test model: An aeroelastic antisymmetric analysis using a direct interpolation 

for AERO paneling. 


Mode 

Instability 

number 

k-SOLN 

p-k 

ASE 

Velocity, 

keas 

Frequency, 

rad/sec 

Velocity, 

keas 

Frequency, 

rad/sec 

Velocity, Frequency, 
keas rad/scc 

Fuselage first bending 

FI 

626.9 

74.3 

632 

73.7 

597 79.8 

Wing first bending 

F2 

... 

... 

— 

... 

820 101.9 

Fin first bending 

D1 

795 

0 

... 

... 

780 0 


Analysis notes: 

1) F - Flutter point 

2) D - Divergence point 

3) Mach = 0.90 

4) Altitude = Sea level 
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10162 


Figure 30. STARS ATM-k flutter analysis — damping (g*), frequency (P), velocity (v) plot, 
antisymmetric case, using direct interpolation where g' = g x 200. 
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Velocity, V, keas 


° *2 / 10 ° Q 


Damping, 


SB 

siis^ 


Mode 7 


Mode 6 


Mode 5 


Mode 1 


Mode 4 


Mode 2 


Mode 3 


Damping, g' 


Velocity, V, keas 


Figure 31. STARS ATM-pk flutter analysis — damping (g*), frequency (p), velocity (v) plot 
antisymmetric case, using direct interpolation where g’ = g x 200. 


3 















Velocity, V, keas 



Figure 32. STARS ATM-ASE flutter analysis — damping (a), frequency (p), velocity (v) plot, 
antisymmetric case, using direct interpolation. 
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7.4 ATM: Aeroelastic analysis 


(STARS-ASE-CONVERT) 


These input data enable appropriate reordering of generalized matrices. Thus, the three perfect rigid 
body modes (O pr ) are placed in the front, followed by eight elastic modes (O e ) and two rigid control 
modes (O c ) for the ASE solution. For the flutter analysis, only the eight elastic modes (O e ) are used for 
the solution. 

STARS-ASE-CONVERT input data: 


$ CONUERT FILE FOR RSE SOLUTION 
13 

$ tIODRL SELECTION RNO ORDERING 
13, I 
11,2 
15,3 
1,1 

2.5 

3.6 

1.7 

5.8 

6.9 
7,10 

10,11 

16,12 

17,13 


$ CONUERT FILE FOR RSE FLUTTER AND DIUERGENCE SOLUTION 
8 

$ nODRL SELECTION AND ORDERING 
I , I 
2,2 
3,3 
1,1 

5.5 

6.6 
7,7 

10,8 
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7.5 ATM: Aeroservoelastic analysis 


(STARS-ASE-PADE) 


These input data effect curve fitting of unsteady aerodynamic forces employing Pade polynomials. 
The state- space matrices are also formed in this module. Version I of the input file pertains to the ASE 
flutter solution, whereas version II corresponds to the subsequent ASE frequency response and damping 
solution. 

STARS-ASE-PADE input data: 


S ATfl ASE FlUTTEB ANRLVSIS, 0.9 fIRCH AT SEA LEUEL - UEASION 
0, 8, 0, 0, 0, 10, 2, 1.0, 871.0, 3.2, 

$ TENSION COEFFICIENTS 
0.1 0.2 
i GENERALIZED MASS 


0. 2570E + 00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.7309E+01 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 1387E+01 

0.0 

0.0 

0.0 

0.0 

0.0 

0. 1882E+01 

0.0 

0.0 

0.0 

0.0 


0.635IE+01 

0.0 

0.0 

0.0 



0 . 1 188E+0 1 

0.0 

0.0 




0. I005E+00 
0 . 6790E+00 

Q.O 






$ GENERALIZED DROP I NG 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
0 . OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
J NATURAL FREQUENCIES 

0 . 6393 1 008E + 02 0 . 7821 6623E+02 0 . 9 1 931026E+O2 0 . 1 8058112E+03 
0. lfl?30l1IE + 03 0 . 20388996E+03 0 .22502991 E+03 0 . 3718801 OE+03 
$ UELOCITJES FOR FLUTTER AND DIUERGENCE ANALVSIS 
I .0 


I DATA 
0, S3 


100.0 

200.0 

300.0 

100.0 

500.0 

600.0 

700.0 

800.0 

900.0 

1000.0 

1 100.0 

1200.0 

1210.0 

1220.0 

1230.0 

1210.0 

1250.0 

1260.0 

1270.0 

1280.0 

1290.0 

1300.0 

! 100.0 

1500.0 

1600.0 

1700.0 

1800.0 

1900.0 

2000.0 

2050.0 

2100.0 

2150.0 

2200.0 

2250.0 

2300.0 

2350.0 

2100.0 

2150.0 

2500.0 

2550.0 

2600.0 

2650.0 

2700.0 

2750.0 

2800.0 


146 



0 , 


0 


2850.0 

2900.0 

2950.0 

3000.0 

3200.0 

3400.0 

3600.0 


$ fiTtl fiSE ANALYSIS, 0.9 NACH AT 40K FEET - UERSiOH II DATA 

3, 8, 2, 0, 2, 10, 2, 0.217, 871.0, 3.2, 

$ TENSION COEFFICIENTS 
0.4 0.2 

$ GENERALIZED NASS 


0 . 7878E+02 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 470 IE+01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.183IE+05 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 2570E+00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.7309E+01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0 . 1 387E+01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 




0. 1882E + 01 

0.0 

•0.0 

0.0 

0.0 

0.0 

0.0 





0.6351E+OI 

0.0 

0.0 

0.0 

0.0 

0.0 






0. 1188E+0I 

0.0 

0.0 

0.0 

0.0 







0. 1005E+00 

0.0 

0.0 

0.0 








0 . 6790E+00 

0.0 

0.0 









0, 3998E+01 

0.0 











0.411 9E + 00 

t GENERALIZED DANPING 

o.ooooooooe+oo O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 

$ Natural Frequencies (radians) 


0.0 


0.0 0.0 0.0 0 . 63931008E+02 

0.782I6623E+02 0 .91934026E+02 0 . 18058112E+03 0, I8730I4IE+03 
0 . 20388996E + 03 0 . 22502991 E+03 0 . 37488010E+03 0.0 

0.0 

$ PHI, THETA, PS I , US, US, US, PS, OS, RS, PHIO, THAD, PSID, HDOF 
0.0, 0.0, 0.0, 871.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, -3 

$ SENSOR DATR 
0 


300.00 

0.0 

50.0 




0.0 

0.0 

0.0 

1 .0 

0.0 

0.0 

300.00 

0.0 

50.0 




0.0 

0.0 

0.0 

0.0 

0.0 

1 .0 


STARS-ASE-PADE analysis results: 

The state-space matrices generated in this module by the Version I data file are utilized for the flutter 
solution, the results of which are given in table 17. Results derived through utilization of Version II data 
are used for subsequent ASE frequency response and damping analysis in the next section. 
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7.6 ATM: Aeroservoelastic analysis 


(STARS-ASE-FRESP) 


Input data presented herein pertain to the frequency response analysis of the ATM at Mach 0.9 and 
40,000 ft altitude. Thus, phase and gain margins for the equivalent open-loop case as well as damping 
and frequency values for the closed-loop configuration are generated from this module. Figure 33 shows 
the block diagram for the ATM lateral mode analog control system. 

STARS-ASE-FRESP input data: 

Open-loop case - 

S RTN ANTISVHNETRIC THREE RIGID, EIGHT ELASTIC, AND TUO CONTROL NODES 


18, 12, 

1, 

1, 58, 

0.0, 

0.0 

10 1 

10 

6 10 1 

2 1 1 


% BLOCK 

CONNECT IUITV 


1 

3 

0 

0 

0 

2 

10 

0 

0 

0 

3 

7 

0 

0 

0 

1 

6 

0 

0 

0 

5 

0 

0 

0 

1 

6 

0 

0 

0 

2 

7 

0 

0 

0 

3 

8 

0 

0 

0 

1 

9 

5 

0 

0 

0 

10 

8 

0 

0 

0 

$ TRANSFER 

FUNCTION DESCRIPTIONS 

1 

1 

2 



2 

1 

2 



3 

2 

2 



1 

2 

2 



5 

1 

3 



6 

1 

3 



7 

1 

1 



6 

1 

1 



9 

1 

3 




10 I 1 

$ LISTING OF P0LVN0N I AL COEFFICIENTS 

1 . 2Q00E+02 .0000E+00 .OOOOE+OO 
0 . 2000E+02 . I 000E+0I .0000E+00 

2 . 2000E+02 . 0000E + 00 .OOOOE*QO 
0 . 2000E+02 . I000E+0I .OOOOE+OO 

3 . 5000E+01 . 1000E + 00 .0000E+00 
0 . OOOOE+OO . 1000E+01 .0000E+00 

4 . 0000E+00 . 1 000E+01 ,O0O0E*OO 
0 , 1 OOOE+OO . 1 Q00E-*01 .0000E+00 

5 . 1 877E+05 . OOOOE + OO .0000E+00 
0 . 1 877E*05 . 1 930E+03 . 1000E*01 

6 . I 677E«-05 . OOOOE+OO .O0OOE*OO 
0 . 1877E+05 . 1 930E+03 . 1000E*01 

7 . 1 000E*01 . 0000E*00 .OOOOE+OO 
0 . 1 Q00E+01 .OOOOE+OO .OOOOE+OO 

8 . 1 000E+01 .OOOQE+OO .0000E+00 
0 . 1 OOOE+OI . 0000E+00 .0000E*00 

9 . lOOOE+OO . OOOOE+QO .OOQOE+OO 
0 . 1000E*02 . 1 100E + 02 .1000E+01 

10 . 1000E*0I .0000E«00 .0000E*00 

0 . 1000E+01 . 0000E+00 .0000E+00 
S GAIH INPUTS FOR EACH 8L0CK 

1 . 1000E+0I 2 . tOOOE*OI 3 .tOOOE+OI 1 .10Q0E+00 5 .I000E+0I 

6 . 1 OOOE+OI 7 . 1000E+01 8 .1000E*0! 9 .1000E+01 10 .I000E+01 

J SPECIFICATION FOR SVSTEN OUTPUTS 

7000000000000000 
0 0 0 0 0 0 

J SPECIFICATION FOR SVSTEN INPUTS 

78000000 
J CONNECTION DETAILS FRON PLANT TO BLOCKS 
2 2 0 

7 1 

8 2 

1 1 

2 2 

$ FREQUENCV RANGE SPECIFICATIONS 
0.1, .9, 200 

1.0, 9. 200 

10. 90. 200 


148 



too. 500. too 
J LOOP DEFINITIONS 
t 0 

4 8 

9 7 
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STARS-ASE-FRESP input data: 

Closed-loop case - 

i ATtl ANT I -SVMHETRI C-THE ROLL RNO VflU CLOSED LOOP CASE 
18, 12, 1, 1, 53, 0.0, 0.0, 3, 3 

12 1 10 6 10 0 0 I I 

$ BLOCK CONNECT IUITV 

1 3 0 0 0 0 0 

9 5 0 0 0 0 0 

It 0 0 0 0 3 9 

3 7 0 0 0 0 0 

5 0 0 0 10 0 

7 0 0 0 3 0 0 

2 10 0 0 0 0 0 

1 6 0 0 0 0 0 

12 0 0 0 0 10 1 

6 0 0 0 2 0 0 

8 0 0 0 1 0 0 

10 8 0 0 0 0 0 

$ TRANSFER FUNCTION DESCRIPTIONS 

1 1 2 

2 1 2 

3 2 2 

12 2 

5 1 3 

6 I 3 

7 I 1 

8 l 1 

9 1 3 

10 I I 

11 1 I 

12 I I 

J LISTING OF POLVHOn I RL COEFFICIENTS 

1 . 2000E+02 . OOOOE+OO .OOOOE+OO 
0 . 2000E+02 . IOOOE+OI .OOOOE+OO 

2 . 2000E+02 . OOOOE+OO .OOOOE+OO 
0 . 2000E+02 . 1QQOE+01 .OOOOE+OO 

3 . 5000E+01 . 1000E+00 .OOOOE+OO 

0 .OOOOE+OO . 1000E+01 .OOOOE+OO 

1 .OOOOE+OO . 1000E+Q 1 .OOOOE+OO 
0 , 1000E+00 . 1000E+0 1 .OOOOE+OO 

5 . 1877E+05 .OOOOE+OO .OOOOE+OO 
0 , I877E+05 , I 930E+03 .IOOOE+OI 

6 . 1877E+05 .OOOOE+OO .OOOOE+OO 
0 . 1877E+05 . 1930E+03 . IOOOE+OI 

7 .IOOOE+OI .OOOOE+OO .OOOOE+OO 
0 .IOOOE+OI .OOOOE+OO .OOOOE+OO 

8 .IOOOE+OI .OOOOE+OO .OOOOE+OO 
0 .IOOOE+OI .OOOOE+OO .OOOOE+OO 

9 . 1000E+00 .OOOOE+OO .OOOOE+OO 
0 . 1000E+02 . 1 1 00E+02 .IOOOE+OI 

10 . IOOOE+OI .OOOOE+OO .OOOOE+OO 
0 . IOOOE+OI .OOOOE+OO .OOOOE+OO 

11 . IOOOE+OI .OOOOE+OO .OOOOE+OO 
0 . IOOOE+OI .OOOOE+OO .OOOOE+OO 

12 .IOOOE+OI .OOOOE+OO .OOOOE+OO 

0 .IOOOE+OI .OOOOE+OO .OOOOE+OO 
J GRIN INPUTS FOR EACH BLOCK 

1 .IOOOE+OI 2 .IOOOE+OI 3 .IOOOE+OI 1 .1000E+00 5 .IOOOE+OI 

6 .IOOOE+OI 7 .IOOOE+OI 8 .IOOOE+OI 9 .IOOOE+OI 10 .IOOOE+OI 

11 - . I00E+01 12 - , 1 00E+01 

$ SPECIFICATION FOR SVSTEfl OUTPUTS 

7600000000000000 
00000000 
i SPECIFICATION FOR SVSTEN INPUTS 

78000000 
$ CONNECTION DETAILS FRON PLANT TO BLOCKS 

2 2 0 

7 I 

8 2 

1 1 

2 2 

S LOOP DEFINITIONS 
0 0 
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Filter 


Roll rate sensor 



Filter Yaw rate sensor 

10165 


Figure 33. ATM lateral mode analog control system. 


STARS-ASE-PADE analysis results: 

Figures 34 and 35 depict the lateral loop gains for the roll and yaw modes, respectively. The gain 
margins are tabulated in table 18. 


Table 18. ATM phase and gain margins. 


Mode 

Phase crossover, rad 

Gain margin, dB 

Roll 

19.88 

79.59 

Yaw 

2.80 

-4.62 


The closed loop damping and frequency plots are shown in figure 36. 
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Figure 36. Closed-loop equivalent velocity vs. damping and frequencies. 


Ames Research Center 
Dryden Flight Research Facility 
National Aeronautics and Space Administration 
Edwards, California, September 28, 1990 
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APPENDIX A — PREPROCESSOR MANUAL 


The preprocessor program PREPROC forms an integral part of the set of routines that compose 
the STARS program. It has been developed to automate generation of finite element models and 
corresponding data files. Instead of defining a complete structure by independently describing each 
node and element, the preprocessor allows the formation of such data automatically for a structure. 

The preprocessor minimizes the information required from the user and eliminates the time-consuming 
task of data editing, enhancing the efficiency of the STARS program. 

To run the preprocessor, the user may type the command GRUN followed by the command 
PREPROC, and immediately after this the program will prompt with a list of different terminals. 

The user may then choose the type of terminal to be used, namely E/S PS 390, Tektronix, and var- 
ious compatible terminals. Next, the user will be prompted with menu options in a progressive fashion. 
At any level of the menu, the user can exit by entering Control-Z. 

Only a brief description of the primary menu is given next since, because of the interactive nature of 
the program, the user is automatically exposed to more extensive details. 



PREPROC MENU 


MENU OPTIONS: 

0 STOP 

stop the program 

1 COMPUTER-AIDED DESIGN 

generate graphics objects 

2 PROPERTIES AND ANALYSIS 

specify STARS data 

3 READ 

read STARS data file 

4 WRITE 

write STARS data file 

5 DELETE 

delete the current structure 
1 COMPUTER-AIDED DESIGN 
DESIGN OPTIONS: 

0 QUIT 

quit this menu 

1 LINES 

generate line segments 

2 SURFACES 

generate surface segments 

3 SOLIDS 

generate solid segments 

4 SYNTHESIS 

generate surfaces from existing line segments 

5 REPRODUCE 

generate new segments using existing ones 

6 DRAW 

plot the current structure 

7 EDITOR 

modify existing data 
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1.1 LINES 


0 QUIT 

quit this menu 

1 STRAIGHT LINE 

generate straight line segment 

2 PARABOLIC 

generate parabolic line segment 

3 CIRCULAR CURVE 

generate circular line segment 

4 ELLIPTIC CURVE 

generate elliptical line segment 

1.2 SURFACES 

0 QUIT 

quit this menu 

1 SIMPLE SURFACE 

generate four node surface segment 

2 COMPLEX SURFACE 

generate nine node surface segment 

3 ELLIPTICAL SURFACE 

generate elliptical surface segment 

1.3 SOLIDS 

0 QUIT 

quit this menu 

1 8-POINT SOLID 

generate eight node segment 

2 ELLIPTICAL SOLID 

generate solid cylinder 

3 4-POINT SOLID 

generate four node segment 

4 6-POINT SOLID PRISM 

generate six node segment 



1.4 SYNTHESIS 


0 QUIT 

quit this menu 

1 ARC: LINE SEGMENT ~> SURFACE 

generate surface segments by moving a line segment along a curve 

2 GLIDE: 2 LINE SEGMENT -> SURFACE 

generate surface segments using two line segments 

1.5 REPRODUCE 

0 QUIT 

quit this menu 

1 COPY 

reproduce by method of direct copying 

2 MIRROR 

produce a mirror image 

3 ROTATE 

reproduce by rotating the original about an axis 

1.6 DRAW 

The preprocessor will draw the generated structure on a standard terminal with multiple options. 
2 PROPERTIES AND ANALYSIS 

This option enables automatic generation of a complete STARS data set for a structure in which 
the user is prompted for appropriate data. 
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APPENDIX B — POSTPROCESSOR MANUAL 


The POSTPLOT program is designed to provide graphic depiction of analysis results pertaining 
to the three major modules, namely SOLIDS, AEROS, and CONTROLS. This is effected by the main 
command GRUN, followed by the POSTPLOT command. The program runs on a variety of terminals 
such as E/S PS 390, Tektronix, and other PLOT 10/PHIGS -compatible machines. 



1.0 Basic Menu 


1.1 On-off switches 

1.2 Load Database 

1 . 3 Delete Database 

1.4 Miscellaneous 

1.5 Exit 

1.1 On-off switches 

1.1.1 Original structure 

1.1.2 Deformed structure 

1.1.3 Dynamic response 

1.1.4 Displacement as a function of time 

1.1.5 Stress as a function of time 

1.1.6 Node number 

1.1.7 Element number 

1.1.8 Element group 

1.1.9 Depth clipping 

1.2 Load Database 

1 .2. 1 Deformed or mode shape 

1 .2.2 Rendering deformed or mode shape 

1.2.3 Dynamic response 

1.2.4 Rendering stress 

1.2.5 Rendering deformation 

1.2.6 Displacement as a function of time 

1.2.7 Stress as a function of time 
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1.2.8 Node numbers 


1.2.9 Element numbers 

1.2.10 Numerical renumbering 

1 .2. 1 1 Aerodynamic paneling plots 

1.2.12 Interpolated mode shape for aerodynamic loads calculations 

1.2.13 Aerodynamic pressure distribution 

1.2.14 Frequency-damping-velocity plots, k, p-k, and ASE solutions 

1.2.15 Phase and gain plots as a function of frequency for analog and digital 
systems 

1.2.16 ASE damping and frequency plots as a function of velocity 

1.3 Delete Database 

Essentially any one of the loaded databases, given above. 

1.4 Miscellaneous 


A host of additional options. 



APPENDIX C — SYSTEMS DESCRIPTION 


The STARS computer program is set up using a main directory and many subdirectories. The setup 
described in this section uses the directory names employed on the DEC VAX computer system at 
NASA. The top-level directories are shown in figure 37. The [KGUPTA.STARS] is the main directory 
which contains the five major subdirectories named COMMANDS, SOURCES, OBJECTS, EXECU- 
TIONS, and TESTCASES. The COMMANDS subdirectory contains the command files which are used 
to guide the user in running the STARS program system. The SOURCES subdirectory contains the 
source elements for the program. It is further subdivided into the SOLIDS, AERODYNAMICS, CON- 
TROLS, and GRAPHICS subdirectories. The OBJECTS subdirectory contains the object elements re- 
quired for creating the execution elements. The object elements have been combined into various object 
libraries to ease the linking process. The EXECUTIONS subdirectory contains the execution elements 
used to run the program. The TESTCASES subdirectory contains a variety of representative example 
problems that facilitate the learning and debugging of the program. 
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Figure 37. STARS systems description, 




















APPENDIX D — NONLINEAR MULTIDISCIPLINARY SIMULATION 


Present effort is directed towards generating a number of additional analysis modules for effective 
simulation of nonlinear systems. These current developments include efficient fluids unstructured grid 
generation, and finite element CFD and heat transfer capabilities. Figure 38 provides details of the vari 
ous modules of the emerging code, which is currently being utilized (ref. 18) for simulation of highly 
nonlinear systems such as the national aero-space plane (NASP). 



Figure 38. Nonlinear multidisciplinary simulation (modular environment). 
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